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ABSTRACT 

A plane-stress analysis is presented for an edge-stiffened 
isotropic or orthotropic elastic rectangular plate subjected to 
prescribed loads and prescribed temperature distributions. 

Along one or more. of the edges the stiffeners are assumed to 
rigidly maintain a prescribed shape (e.g., straight), while the 
stiffeners along the remaining edges are assumed to have 
negligible flexural stiffness. All four stiffeners are assumed 
to be uniform and to possess finite axial stiffness. The plate 
edges are assumed to be integrally attached to the stiffeners 
along the originally straight centroidal axes of the stiffeners. 

This work represents a generalization of earlier studies in 
which only boundary loadings (rather than boundary displacements) 
were prescribed. 

The analysis is by means of Fourier series. Numerical 
results are presented for the specific cases of an isotropic 
square stiffened plate with one edge, two opposite edges, or all 
four edges held straight. In all three of these cases the plate 
was assumed to have a pillow-shaped temperature distribution. 

In the second case stretching forces parallel to the non-straight 
edges were also considered. 
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INTRODUCTION 


In reference 1 a plane-stress analysis, by means of Fourier series, 
was presented for an elastic rectangular plate with four edge stiffeners, 
subjected to any equilibrium system of boundary loads and any prescribed 
temperature distribution. The plate was assumed to be isotropic or ortho- 
tropic, with elastic constants independent of position, and, if ortho- 
tropic, with axes of elastic symmetry parallel to its edges. The 
stiffeners were assumed to have finite extensional stiffness but negligible 
bending stiffness, and the attachment between plate edge and adjoining 
stiffener was assumed to be along the originally straight centroidal axis 
of the stiffener. The stiffeners were assumed to be either uniform in 
cross section or tapered in such a way as to result in any prescribed 
variation of stiffener cross-sectional stress along the length of stif- 
feners. This structure is shown schematically in figure 1. 

In this earlier work, just described, the boundary conditions 
were entirely those of prescribed normal and shearing loadings along the 
outer periphery of the stiffeners. The purpose of the present paper is 
to generalize this previous work by assuming that along one or more edges 
the boundary condition of prescribed normal loading is replaced by one 
of prescribed normal displacement - i.e. along one or more of the edges 
the stiffener is assumed to be rigidly held in a prescribed shape. We 
can imagine that such a stiffener has been endowed with infinite flexural 
stiffness after being bent to its prescribed shape, while its axial 


1 



stiffness remains finite. Along such an edge, one may still prescribe a 
distribution of external normal forces along the stiffener, but only the 
resultant force T and moment M of this distribution are significant, since 
the stiffener shape has been assumed to be rigidly fixed. On the remain- 
ing edges the boundary conditions continue to be those of prescribed normal 
and shear loading applied to perfectly flexible stiffeners. Only the case of 
constant - area stiffeners is considered. 

By this generalization the range of applicability of the analysis is 
widened to include cases in which, for example, certain edges of the plate 
are forced to remain straight. As an illustration, one can cite an interior 
bay of an airplane shear web, bounded by spar caps along two opposite edges 
(the top and bottom) and upright along the other two edges. Continuity of 
the plate across the uprights suggests that when the adjacent bays are 
identical to the one under consideration and similarly heated or loaded 
the plate edges along the uprights will tend to remain straight. Thus 
such a bay would correspond to an edge-stiffened rectangular plate with 
two opposite edges kept straight. Similarly, an interior skin panel of a 
multi-spar multi-rib wing could correspond to the case of all four edges 
held straight. 

The most general results of the present analysis are in the form of 
equations for stresses in the case of shape, resultant external normal 
force, and resultant external moment prescribed along (a) one edge, 

(b) two opposite edges, and (c) all four edges. The use of these equa- 
tions is described in detail. In addition, numerical examples are given 
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corresponding to each of these cases. These numerical examples are for 
doubly symmetric square plates with one edge, two opposite edges, or all 
four edges held straight, and, for the most part, a prescribed "pillow- 
shaped" temperature distribution, i.e. a temperature distribution in which 

k 

the stiffener temperatures are constant at one value and the plate experi- 
ences a temperature rise, relating to the stiffeners, which varies as a 
half sine wave in both directions across the plate. For the case of two 
opposite edges held straight on overall tension loading is also considered 
which stretches the plate in the direction of the unconstrained edges, 
while the constrained edges remain straight and parallel to each other 
but move apart. 

For the sake of simplicity, the term "temperature distribution" has 
been and will often again be used in this paper. It should be understood, 
however, that what is meant is a distribution of temperature rise with 
respect to some datum temperature distribution at which the structure is 
assumed to be stress-free in the absence of applied loads. Usually, the 
datum temperature distribution is a uniform one, and then it becomes cor- 
rect to speak of thermal stresses due to a temperature distribution rather 
than to a temperature-rise distribution. 

The structure, loading, boundary conditions, and results are described 
in more detail in the following sections. The symbols are defined where 
they are first used, and the definitions are also compiled in appendix A for 
convenience of reference. Those details of the analysis not essential for 
the understanding and use of the results are given in appendixes B through E. 


DETAILED DESCRIPTION OF STRUCTURE 


Geometry and Coordinate System 

The combination of the plate and edge stiffeners is as shown 
schematically in figure 1. The plate has a length a and a width b. 
The position of any point in the plate is given by its coordinates x 
and y in a Cartesian coordinate system whose origin is at a comer 
of the plate and whose axes coincide with two adjacent edges, as shown 
in the figure. The cross-sectional areas of the stiffeners are denoted 
by A^, A . £ , A^ and A^ for the stiffeners located at x = 0, x = a, y = 0 
and y = b, respectively. In the analysis the stiffener axes are assumed 
to coincide with the plate edges, but in figure 1 the stiffener axes are 
shown slightly offset from the plate edges for clarity. 

Thermal Strains 

We assume that a datum temperature distribution exists for which 
the unloaded structure is free of internal stress. When the structure 
is in this datum condition the strains will be considered zero. For any 
other temperature distribution, those strains that would be produced if 
the thermal expansions of every infinitesimal element were permitted to 
occur without restraint from neighboring elements will be called the 
"thermal strains." The temperature distribution and coefficients of 
expansion are assumed to be known throughout the stiffeners and plate, 
hence the thermal strains are assumed to be everywhere known. 
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The notation for thermal strains is indicated in figure 2 and is 


as follows: e-^Xy), e 2 (y), e.j(x) and e^(x) are the thermal strains in 

the stiffeners located at x = 0 , x = a, y = 0 and y = b* respectively; 

e (x,y) and e (x,y) are the thermal strains in the plate in x and y 
x y 

directions, respectively. The strains are positive for elongation. 

Because the analysis allows for an orthotropic plate, e and e need 

x y 

not equal, but since the x and y axes are parallel to the directions of 


elastic symmetry, there is no thermal shear strain e 

xy 


Stress-Strain Relationships and Elastic Constants 

The notation for the internal forces is indicated in figure 3. 

P^(y) , P 2 (y), P 3 (x), P^(x) denote the cross-sectional tensions and 

e^y), e 2 (y), e 3 ( x )» £4 00 the total strains (thermal plus elastic) 

in the stiffeners located at x = 0 , x= a, y = 0 , y = b, respectively. 

The plate stress resultants (force per unit length) are represented by 

N (x,yD and N (x,y) for normal stress and N (x,y) for shear stress, 
x y xy 

as shown in figure 3. The corresponding total strains are symbolized 
by e (x,y), e (x ,y) and Y VTT ( x »y)- 

x y Ay 

In terms of this notation, the assumed stress-strain relations 
for the stiffeners are 

E i- e i + r-f- tt-1,2,3,4) (1) 

i i 

with Young's moduli and E 2 independent of y, E^ and E^ independent 
of x. The stress-strain relations of the plate are taken as: 
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e - e + C.N - C,N 
x x 1 x 3 y 


e - e + C-N 
y y 2 y 


C 3 N x 


( 2 ) 


y - C.N 
'xy 4 xy 


where the elastic constants C^, C 2 , and are independent of x . and y. 
If the plate is homogeneous . and isotropic, with thickness h. Young's 
modulus E, and Poisson's ratio v then 


C 1 = c 2 - (Eh) 

C 3 = v (Eh) -1 (3) 

C 4 = 2(l+v)(Eh)“ 1 


The assumption that the elastic constants are independent of spatial 
coordinates also implies the assumption that they are independent of 
temperature. 


Boundary Conditions 

Figure 1 shows the boundary conditions purely of prescribed 
loading which were employed In reference 1. This loading consists of 
end forces P^(0), P^(b) etc. applied to the centroids of the end cross 
sections of the stiffeners, and distributed external shear flows 
q^(y), ^ 2 ^* q^(x), and distributed external normal forces 
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N^(y), N 2 (y) , N^Cx), N^(x) applied to the stiffeners. The distributed 
loadings have the dimensions of force per unit length. The shear flows 
are assumed to be acting along the centroidal axes of the stiffeners, 
although for charity they are shown displaced from these axes in 
figure 1 . Because the edge stiffeners have zero bending stiffness, the, 
normal force is transmitted directly through them into the edge of the 
plate* 

In the present paper the boundary condition of prescribed normal- 
force distribution is replaced by that of prescribed and rigidly main- 
tained stiffener shape and prescribed external normal-force resultants 
along one edge, the opposite edges, or all four edges. (Along the re- 
maining edges two stiffeners are still assumed to be perfectly flexible.) 
These cases are shown schematically in figure 4 , in which the cross-hatched 
stiffeners are those whose shapes are rigidly prescribed. T^ through T^ 
are the prescribed external normal force resultants , and through 
the prescribed external moment resultants, acting on these stiffeners. 
Letting u(x,y) and v(x,y) denote , the displacement components in the x and 
y directions, respectively, we can describe the boundary condition of 
prescribed shape mathematically as follows for case (a) of figure 4 : 




x-0 


is a prescribed function of y. 


Similarly, the shape boundary conditions for case (b) are: 


g2 u g2 u 

(t— 7) and (77) are prescribed functions of y. 


'3y' 


x-0 


k 3y z/ 

J x=a 


And those for case (c) are: 
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( i% 


and 


x«0 


,d £ Uv 


are prescribed functions of y 


x=a 


, 1 2 


3 2 v, 


(S=t) and (f-r) 


3x 


y**0 


3x^ 


y=b 


are prescribed functions of x. 


Note that the prescribed shape of an edge is described by means of its 
curvatures. If an edge is held straight the prescribed curvatures are 
zero. 

The boundary condition of prescribed shape along an edge can be 
achieved, hypothetically, by forcibly bending the edge stiffener (in 
its attached condition) to the prescribed curvatures, and then endowing 
it with infinite flexural stiffness while its axial stiffness remains 
finite. If the now flexurally-rigid stiffener is free of all external 
forces, the values of T and M shown in figure 4 for that stiffener became 
zero; the system of internal normal forces acting between that stiffener 
and rest of the structure will then constitute a self-equilibrating 
system, and a mathematical statement to that effect must be included in 
the analysis as a boundary condition. On the other hand » if the stiffener 
is not completely free of external force then the internal normal forces 
along the edge of the stiffener plus the external T and M applied to 
that stiffener must together be self-equilibrating. 

Of course, the total system of external forces must at the outset 
constitute an equilibrium system. 
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Comer Conditions 


The stiffeners will be assumed to be, hinged where they meet at 
the comers. This assumption is really superfluous in the case of figure 
4a or 4b, for in these cases one of the two stiffeners meeting at every 
comer is a perfectly flexible one, incapable of developing any bending 
moment. The assumption is . significant only in the case of 4c, where 
at each comer two rigid stiffeners meet. 

If the loading and temperature distribution for the case of 
figure 4c are such that there is no tendency for the stiffeners to undergo 
relative rotation at the comers (as was true in the numerical examples 
considered in this paper), then, of course, the plate stresses and stif- 
fener tensions are the same for rigid joints as for hinged joints. 
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FOURIER SERIES EXPANSIONS FOR PRESCRIBED BOUNDARY LOADS, 
PRESCRIBED BOUNDARY CURVATURES, AND PRESCRIBED 
THERMAL STRAINS 


The results of the present analysis, to be discussed shortly, con- * 
slst of formulas for the stiffener and plate stresses In terms of the given 
loading) the given boundary curvatures, and the known thermal strains. How- 
ever, these prescribed quantities do not appear explicitly in these formulas; 
it is rather the Fourier coefficients of these quantities which are required. 
In anticipation of this requirement it is assumed that the given quantities 
can be expanded as Fourier series, with known coefficients, in the forms given 
below. 


Prescribed Boundary Loadings 

When the prescribed normal-force distributions N^(y), ^(y), N^(x) , 
N^(x) exist, that is when their existence is not pre-empted by boundary con- 
ditions of prescribed shape, then they will be assumed known in the form of 
the following series: 


V y) ‘ l, K sl " S b £ 

n=l 

N 2 (y) - l 

n~l 

No (x) = J B-'" sin 2221 

3 m-1 m 

N. (x) ** l B ,m sin 222 

4 m a 

m=l 


(0 < y < b) 


(0 < y < b) 


(0 < x < a) 


(0 < x < a) 


( 4 ) 
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m. 


where 


■5 


K-z U 


(y) sin 


mn 


dy^ etc. 


(5) 


And the shear flows will be assumed to be given in the following form: 


N 


q x (y) 

- i 



n=»0 


N 


q 2 (y) 

- Z 

<5 


n**0 


M 


q 3 00 

- Z 

Q... 


m»0 

m 


M 


V*) 

- Z 

m*0 

<c 


cos 


tnrv 


cos 


BI2. 


cos 


rnirx 


cos 


mux 


( 6 ) 


where 


2-6 


% 


n(J 


q x (y) 


cos 


nir\ 


dy 


etc. 


(7) 


and 6 is Kronecker’s delta. 

nO 

When any of the loadings N 1 (y), N 2 (y), N 3 (x), or N^Cx) exists, it 

is transmitted, through the assumedly perfectly flexible stiffener, into the 

plate. Thus the corresponding Fourier coefficients B', B". B ,,, , or B"" 

n n m ?n 

describe not only the externally applied normal loading but also the internal 

distributed tension N (0,y) , N (a,y), N (x,0) , or N (x.b), respectively, 

* x y y 

acting mutually between the stiffener and the edge of the plate. 
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Finite upper limits M and N are shown for the summation indexes in 
these (and subsequent) series in expectation of the fact that, for practi- 
cal computational reasons, it will normally be necessary to use truncated 
rather than infinite series. 


Prescribed Boundary Curvatures 


When curvatures are prescribed, they will be assumed to be given 
by one or more of the following series, depending on which edges have the 
prescribed curvatures: 


where 


<&> 


( 2 % 


^x 2 ' 


''Sx 2 ' 


K’ = 
n 


N 


— 

y 

K f sin 

(0 

< 

y 

< 

b) 

(8) 

x=Q 

L 

n=l 

n 

D 








N 









ss 

y 

K" sin 

(0 

< 

y 

< 

b) 

(9) 

x=a 

L 

n=l 

n 

D 








M 


mirx 







= 

y 

K"’ sin 

(0 

< 

X 

< 

a) 

(10) 

y=0 

L 

m=l 

m 

a 





M 

y 

K"" sin 

mitx 

(0 

< 

X 

< 

a) 

(ID 

y=b ” 

L 

m=l 

m 

a 






b 










M 


n sin 
x=0 

n-rry 

b dy 

, etc. 



(12) 


0 


If an edge is held straight, then all the coefficients in the Fourier 
series expansion for the curvatures along that edge will, of course, be 
zero. 
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Prescribed Thermal Strains 


If there are any discontinuities in thermal strain between the 

stiffeners and the edges of the plate, these will be represented by the 

following Fourier series with known coefficients: 

N 
V 

L 

n=l 


^(y) - e y (° »y) = _I n sin ^ 


(0 < y < b) 


N 


(y) - e (a,y) = j T" sin (0 c y < b) 


n=l 


M 


n 


e^(x) - e_(x,o) = J T" ' sin 


m=l 


m 


mux 

a 


(0 < x < a) 


(13) 


M 

e. (x) - e (x,b) = £ T"" sin (0 < x < a) 

4 x m a 

m-1 


where 


2 r 

t ' = r 

n b 


[e (y) - e (0 ,y) ] sin ^ , et 
j -L y d 


c. 


0 


(14) 


The quantity 3 2 e^/3x 2 + 3 2 e x /3y 2 is assumed known throughout the 
plate and representable by the following Fourier series in the open 
region 0<x<a,0<y<b: 


3 2 e 3 2 e 




3x^ 3y‘ 


M N 

r „ „ mnx . nuy 

, ) T sin ■ sm 

L , K mn a b 

m=l n=l 


where 


b a 

r r 


v 

0 


3 2 e 3 2 e 


T mn " lb I i 


3x" 


_x n . mux , nuy , , 

+ — — T“*J sm sm v dxdy 

9y a b J 


(15) 


(16) 
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Integration by parts in the above equation gives the following 

alternate form which permits T to be evaluated from the first partial 

derivatives of e and e instead of the second partial derivatives : 
x y 

b a g e 

_ mrr f f V „ m?rx mry , , 

T nm " ‘ T 7b j \ ST COS T Sln "b d * d y 
0 0 

r f 3e 

x mry . mTrx , , 

g — cos Tjp sin dydx 


n-fT 4 
b ab 


(17) 


0 0 


Equation (17) may be used for piecewise continuous e or e with finite 

y x 

discontinuities provided that 3e /3x and 3e /3y are regarded to be infi- 

y * 

nite in the manner of the Dirac delta function at points of discontinuity. 
If e^ and e^ are continuous in the closed region 0_<x^a, 0 _< y _< b , 
further integration by parts gives 


T — miT 4 

mn a ab 


[e (a,y) cos mir - e (o,y)] sin — dy 

y y b J 


b a 


/ 1| ' 2 4__ 
^a ' ab 


, . iuttx . nTry , , 

e (x,y) sm sm ■ ,r - dxdy 

y a b 3 


0 0 


mr 4 
b ab 


[e (x,b) cos ntr - e (x,o)] sin m , T ^ dx 
x x a 


S ; ab 


a b 
f f 


t \ • mry • m-rrx , , 

e (x,y) sm “H- sm dydx 

x b a 


0 0 


(18) 


Finally, the known quantities (3e^/3x) x _ 0> etc. are assumed to be 


expandable in the following Fourier series: 
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where 


3e 

(— 1 ) 
v 3x 'x=0 


9e 

( — 2.) 

3x x=a 


3e 

( — “) 

Sy ; y=0 

3e 

(— ) 

«y y=b 


N 


tltry 


I v’ sin 
S n b 

n=l 


N 


mry 


l V" sin ~ 
n b 

n=l 


M 


m=l 

M 


y V"’ sin 

m a 


y v m 

m 


sin 


mirx 


m=l 


(0 < y < b) 


(0 < y < b) 


(0 < x < a 


(0 < x < a) 


( 19 ) 


( 20 ) 


( 21 ) 


( 22 ) 


V 1 

n 



3e 

^3x ^x=0 



0 


etc. 


( 23 ) 
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GENERAL RESULTS 


The analyses are carried out in appendixes B through E. The results 
are of two kinds: (a) general results valid for any set of geometric, 

material and loading parameters, and (b) numerical results for specific 
cases. The former 'are presented in this section, the latter in the next 
section. 

The general results of analyses are in the form of Fourier series 
for computing the plate stresses and stiffener tensions. In the subsec- 
tions below the pertinent Fourier series for each of the various stresses 
will be pointed out. In the final subsection (Evaluation of Series 
Coefficients) the procedures for evaluating the Fourier coefficients in 
these series will be summarized. Enough detail will be given to enable 
the reader to make calculations based on the general results without 
having to study the derivations in the appendixes B through E. 

Series for the Interior Plate Stresses and Interior Stiffener Tensions 

The stress resultants N , N and N at interior points of the 

x * y xy r 

plate are given by the double Fourier series in equations (B19) , (B22) 
and (B25) . The stiffener tensions P^(y), etc. are given by the single 
Fourier series in equations (B16) for all but the end sections of the 
stiffeners. 

Series for Stresses along the Plate Edges 

The shear stress resultant N for points along the plate edges can 
be obtained from the same Fourier series, namely equation (B25) , as for 
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interior points. However the normal stresses along the plate edges 
require special single Fourier series, which will now be given. Equa- 
tions (B20) and (B21) can be used for evaluating along the edges 
y = 0 and y = b; equations (B23) and (B24) similarly give along the 
edges x = 0 and x = a. The stresses N x along the edges x = 0 and x = a 
are given by the series for N^(y) and ^(y) , respectively, in equations (4); 
and similarly along the edges y = 0 and y = b are given by N^(x) and 

N^(x) series, respectively, in equations (4). 

The series which are referred to above for the normal stresses 
N x and along the plate edges are valid at all points except the 
corners. Special formulas to be used for the comer values of N x and 
will be given subsequently. 

Series for the End Tensions of the Stiffeners 

Where two perfectly flexible stiffeners meet, as at the point 
(a,o) of figure 4a, the tension at the end of each is merely equal to 
the externally applied load, such as P^(a) or p2(0). Where a perfectly 
rigid and perfectly flexible stiffener meet, as at the point (0,b) of 
figure 4a, the tension at the end of the rigid stiffener is again equal 
to the applied load, such as P^b); but the tension at the end of the 
flexible stiffener, which is provided by the rigid stiffener, is an 
unknown rather than a prescribed quantity. Similarly, where two rigid 
stiffeners meet, as at the comer (a,b) of figure 4c, each stiffener pro- 
vides an unknown tension to the end of the other stiffener. The tension 
which one stiffener exerts on the end cross-section of the adjacent stlf- 
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fener is experienced as a shear force on the end section of the first 
stiffener. This mutual action between the stiffeners at the corner (a,b) 
of figure 4c is illustrated in the sketch below, in which the stiffeners 
are shown disjointed for the sake of clarity. 


P 2 (b) 



Pi* (a) 


Pi* (a) 


The equations for computing the unknown end tensions are as follows : 

i) Equations (C14) and (C15) for the case of figure 4a. 

ii) Equations (D13) through (D16) for the case of figure 4b. 

iii) Equations (Ell) through (E13) for the case of figure 4c. 

Formulas for the Plate Normal Stresses at the Comers 

None of the series for the plate stress resultants N and N thus 

x y 

far presented are valid at the corners of the plate. However, once the 

stiffener end tensions have been determined, these can be used to obtain 

the comer values of N and N . The procedure is as follows for the corner 

x y r 

x = 0, y = 0: From equations (1) and (2) the stiffener strains at this 

comer are 
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while the plate strains are 


e x = e 3 (0) + p 3 (0)/A 3 E 3 
e y = e^O) + P 1 (0)/A 1 E 1 


£ x = e x (0,0) + “ C 3 N y (0,0) 


e = e (0,0) + C„N (0,0) - C„N (0,0) 
y y 2 y * 3 x * 

Equating the stiffener strains and the corresponding plate strains (by 
virtue of the assumed continuity between stiffeners and plate) and 
solving the resulting equations simultaneously for N x (0,0) and N y (0,0) , 
one obtains 


N x (0,0) 


C 2 [e 3 (0)-e x (0,0)] + (^[e^O-e (0,0) ] + C 2 P 3 (0)/A 3 E 3 + (0) /A^ 


c i c 2 - c : 


(24) 


N (0,0) 


C 1 [e 1 (0)-e v (0,0)] + C 3 [e 3 (0)-e x (0,0) ] + 0^(0) /A^ + C 3 P 3 (0)/A 3 E 3 

C 1 C 2 “ C 3 


In a similar manner, one obtains the following formulas for calculating the 
plate normal stresses at the comers (0,b) , (a,0) , and (a,b) : 


N x (0,b) = 


N (0.,b) = 


C 2 Ce 4 (0)-e x (0,b)3 4- (^[e^-e (0,b)] + 0^(0) /A^ 4- C^OO/A^ 


C 1 C 2 “ C 3 


(25) 


C 1 [e 1 (b)-e (0,b)] + C 3 [e 4 (0)-e x (0,b)] + C^OO/A^ + C 3 P A (0)/A4E 4 


C 1 C 2 - C 3 
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N x (a,0) = 


C 2 [e 3 (a)-e x (a,0)] + C 3 [e 2 (0)-e^a,0) ] + C 2 P 3 (a)/A 3 E 3 + C 3 P 2 (0)/A 2 E 2 


C 1 C 2 - C 3 


(26) 


Ny(a»0) - 


C 1 [e 2 (0)->e v (a,0)l + C 3 [e 3 (a)-e x (a,0) ] + C-^W/A^ + C 3 P 3 (a)/A 3 E 3 

~ C 1 C 2 " C 3 


C 2 [e 4 (a)-e x (a,b)] + C 3 [e 2 (b)-e (a,b)] + C 2 P 4 (a)/A 4 E 4 + C 3 P 2 (b)/A 2 E 2 

N * (a ’ b) 5 ^ 7 ^ 


(27) 


N y (a,b) = 


C 1 [e 2 (b)-e (a,b)] + C 3 ie 4 (a)-e x (a,b) ] + C 1 P 2 (b)/A 2 E 2 + C 3 P 4 (a)/A 4 E 4 


C,C„ - Cf 


12 3 

When the plate is isotropic and there is no discontinuity of thermal strains 
between stiffeners and plate at the corners, then the above formulas are 
reduced to 

P (0) P (0) 

v°.°> - ispr » 


Fh P l (0) P 3 (0) 

y°>°> = v vT 1 


(28) 


with similar expressions for the other three comers. 


Evaluation of Series Coefficients 


In order to use the series referred to in the above sections for 

numerical calculation of stresses, one must first evaluate the coefficients 

c T , c", g' , g" , B’ etc. (if these coefficients are unknown), c , g , 1 , 

n n m ’ m n mn °mn’ J mn* 

S n* S n’ S m' ’ anc * S to" a PP earan 8 ^- n th em * The procedure for evaluating these 
coefficients will now be outlined in detail. It will be seen that the first 
four groups of coefficients, namely c^, c^, g^, g^, are the key to all the 
others „ 
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i) The case of figure 4a. - The c'. c" . g' . and g" are defined bv 

a — n n m 

the system of equations (C29) to (C32) and can be determined by solving 
these 2N + 2M equations simultaneously for the c^, c||, g^, and noting 
the definitions in equations (B66) . As an alternative, equations (C46) 
may be solved simultaneously for the and g^; the c^ and c” are then 
obtained directly from equations (C44) and (C45) . The alternative is 
preferable because it requires the solution of only 2M simultaneous equa- 
tions, regardless of how large a value is selected for N. With the c^, 
c" , g' , and g" known, equations (C9) will give the values of B*. 

If the structure, loading and thermal strains are symmetrical about 

the centerline y = b/2, considerable simplification results. The c’ , c", 

n n 

g’ , and g" are then defined by equations (C50) through (C52) . Only the M 
m m 

equations (C51) need to be solved simultaneously. They give the g^, after 
which the g^, c^, and c” are obtained directly from equations (C50) and 
(C53) . Equations (C54) and (C55) may then be used to determine the B^. 

In this case the size of N again does not influence the number of equa- 
tions that have to be solved simultaneously. 

Once the c^, c”, g^, and B^ known, equations (B34) , (B35) , and (B57) 
through (B61) will furnish the remaining coefficients directly, for either 
the symmetrical or unsymmetrical case. 

ii) The case of figure 4b. - The c T . c" . e’. e" are defined bv the 

° n* n °m °m J 

system of equations (D28) through (D31) and can be determined by solving these 

2N + 2M equations simultaneously for the c’, c" , g' , and g". As an alternative, 

n n °nr 6 m * 

equations (D49) may be solved simultaneously for g' and g"; the c’ and c" are 

mm n n 
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then obtained directly from equations (D42) and (D43) . The alternative 

requires the solution of only 2M simultaneous equations , regardless how large 

a value is selected for N. With the c’, c" , g', g" known, equations (D9) wil 

n n* °m m 

furnish the values of B* and B". 

n n 

If the structure, loading and thermal strains are symmetrical about 

both centerlines, considerable simplification results. The g' and e" are 

then defined by equations (D52) and the simultaneous system of (M+l)/2 

equations (D57). The latter may be solved simultaneously for the g^, after 

which the c^, the cj|, B^, Bj| are obtained directly from equations (D56) , 

(D53), and (D52) . The size of N again does not influence the number of 

equations that have to be solved simultaneously. 

With the c’, c" , g' , g" , B 1 , and B" known, equations (B34) , (B35) , and 

(B57) through (B61) will furnish the remaining coefficients directly, for 

either the symmetrical or unsymmetrical case. 

iii) The case of figure 4c. - The B', B" . B"J B"" . c’ . c" . g'. and 
a n’ n’ m* m* n* n* & m’ 

g^ are defined by the system of equations (E23) through (E30) and can be 

determined by solving these 4N 4M equations simultaneously for B^, B”, 

B" * , B"", c r , c", g', g" and noting the definitions in equations (B66) . 

For the case of a square plate (b=a) with structure, loading and 

thermal strains symmetrical about both centerlines (x = a/2 and y = b/2) and 

diagonals, considerable simplification results. The B’, B" , B'", B n ", c f , c" 

n n m m n n 

g^, and g^ are then defined by equations (E31) , (E32) , and (E33) , and the 
M + 1 equations (E34) and (E35) . Only equations (E34) and (E35) have to be 
solved simultaneously. 
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With these known, equations (B34) , (B35) , and (B57) through (B61) 
will furnish the remaining coefficients directly, for either the symmetri- 
cal or unsymmetrical case. 

Limiting case of large stiffener area . - In the appendixes various 
limiting cases are considered in which some or all of the stiffener cross- 
sectional areas are allowed to approach infinity by comparison with the 
plate cross-sectional area. The results of these limiting cases will now 
be described. It will be seen that the calculation procedures for these 
limiting cases are much simpler than for the general case. Except for two 
of the limiting cases considered, it is no longer necessary to solve simul- 
taneous equations in order to determine the Fourier coefficients in the series 
expansions for the stresses. 

(a) The case of figure 4a: For this case all four stiffener 

cross-sectional areas were assumed to approach Infinity simultaneously while 

maintaining constant ratios with each other. The results are contained in 

equations (C61) to (C64) , which give explicit expressions for c^, g^, g^ 

and c', correct to terms of the first degree in 1/ (a^E- .. A- E- ) . More accurate 
n lx 1 1 

results, correct to terms of the second degree in this quantity, are repre- 
sented by equations (C66) , (C67) , (C68) , and (C70) . 

(b) The case of figure 4b: For the configuration of figure 4b 

attention was restricted to the case in which structure, loading, and thermal 
strains are symmetrical about both centerlines (x = a/2, y = b/2) , and five 
different types of limiting case were considered. These five cases are 
described below together with the results obtained form them. 
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(1) =: °°, A^ = A^ remaining finite. The quantities c^, c^, g^, 

and g” are defined explicitly by equations (D52) , (D59) , and (D60) . 

(2) A 3 = A^ 00 , followed by A^ = A^ -*• 00 . The results for this iterated 

limiting process are represented by equations (D52) , (D62) , and (D64) 

(3) A^ = 00 » A 3 = A^ remaining finite. Equations (D52) , (D 66 ) , and 

(D 68 ) result. In this case a system of simultaneous equations (D 66 ) has 
to be solved. 

(4) A^ = A^ -*■ 00 , followed by A 3 = A^ -*■ 00 . Equations (D52) , (D69), 
and (D70) are the result of this limiting process. 

(5) A^(= A^) and A 3 (= A^) approaching infinity simultaneously while the 
ratio of A^ to A 3 remains constant. The results are given by equations 
(D52), (D72) and (D73) . 

In the above-cited results, c' , c" , g' and g" are given correct to terms 

n n m m 

of the first degree in 1/ (a^E^A^E^) . Conditions (2) , (4) , and (5) , 
which are physically identical, also turn out be mathematically identical, 
as an examination of the cited equations will show. 

(c) The case of figure 4c: Here attention was restricted to the 
square case (b = a), with symmetry about both centerlines and diagonals, 
and a particular loading consisting only of prescribed thermal strains. 

The stiffener areas, all equal, were then allowed to approach infinity. 

In view of the highly specialized assumptions no general results will be 
cited for this case. The reader may refer to appendix E for a detailed 
description of the procedure. In this limiting case it is also necessary 
to solve simultaneous equations. 
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NUMERICAL RESULTS 


The foregoing analytical results were used to obtain numerical 
stress data for twelve illustrative problems, ten of them thermal -stress 
problems involving non-uniform temperature distribution without any 
applied loads, and the other two prescribed-force problems corresponding 
t° T]_ = T 2 = T 5 * 0 in figure 4b, with all other loads vanishing and 
temperature uniform. 

In all twelve problems the plate was square (b - a) and isotropic, 
with Young's modulus E, Poisson's ratio v, and thickness h. In each 
problem the two x-wise stiffeners were assumed to be identical to each 
other (i.e. = A^) and the y-wise stiffeners were similarly assumed 

to be identical (i.e. A^ = A^) . Furthermore, in all but one of the 
problems the x-wise stiffeners were taken to be identical to the 
y-wise stiffeners; in presenting the results for those cases the symbol 
A will designate the common value of A^ , A^ , A^ and A^, and the symbol 
X the common value of the area-ratio parameters A^ and X 2 » defined as 
follows: A^ = 4ah/ir 2 A^, A 2 = 4bh/7r 2 A2« The stiffeners were also assumed 
to have the same Young's modulus as the plate. 

In the thermal stress problems, the stiffeners were considered to 
be at a uniform temperature, while the plate was assumed to have a 
pillow-shaped temperature rise, relative to the stiffeners, of the form 
0 sinCux/a) sin(iry/b); thus 6 denotes the temperature rise of the plate 
center relative to the edges. The symbol a will denote the coefficient of 
thermal expansion of the plate material. 
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The table below summarizes the twelve problems for which calculations 
were made. The first column indicates for each problem which of the cases 
is being considered, that is, whether one edge, two edges, or all four 
edges are held straight. The second column indicates the loading condi- 
tion. In this column "PSTD" stands for the pillow-shaped temperature dis- 
tribution described above and connotes a purely thermal-stress problem, 
with no applied forces. The notation T^ = = T which appears for two 

of the problems associated with figure 4b denotes a pair of stretching loads 
applied in the x-direction to a structure of uniform temperature. The third 
column shows the value of Poisson's ratio used in the calculations; except 
for one calculation Poisson's ratio was taken as 0.3. The fourth and fifth 
columns give the values of the area-ratio parameters, defined as follows: 

X^ = 4ah/Tr 2 A^, X 2 = 4bh/7T 2 A^. The value of zero for these parameters refers 
to the limiting case in which the stiffener cross-sectional area approaches 
Infinity by comparison with the plate cross-sectional area. Column 6 
gives the equations used to compute the basic unknowns leading to the 
stresses. Column 7 tells in which figures of the present paper the 
results can be found. The results are in the form of dimensionless plots 
of plate stress and stiffener tension, represented by the solid curves of 
figures 6 through 17. The dashed curves which appear on some of the 
figures are selected results from reference 1 for the case in which all 
four stiffeners are perfectly flexible, included for the sake of comparison. 

Additional information of interest concerning the calculations is 
appended in columns 8, 9, and 10. Columns 8 and 9 give the upper limits 
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(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

No. of 
edges held 

straight 



Loading 

V 

X 1 

X 2 

Equations 

employed 

Results 

M 

N 

IBM 7074 
computing 
time (min.) 

One (fig. 4a) 

PSTD* 

.3 

2 

2 

(C50)- (C55) 

Fig. 6 

30 

59 

21 

VI 

II 

II 

1 

1 



11 

7 

It 

It 

21 

11 

VI 

1 

0 


(C61)-(C64) 

8 

79 

79 

3.5 

II 

II 

B 

0 

0 

If 

9 

tl 

II 

3.5 

Two (fig.4b) 

PSTD* 

.3 

2 

2 

(D52), (D53), 
(D56) , (D57) 

1 

10 

59 

59 

11 

If 

II 

II 

1 

1 

If 

11 

1 

it 

■ 

11 

II 

II 

II 

1 

0 

(D52), (D59), 
(D60) 

■ 

■ 

■ 

2.5 

II 

It 

II 


0 

(D52 ) , (D62), 
(D64) 

■ 

79 

79 

2.5 

ft 

- - i 

T =T =T 
1 2 

II 

I 

■ 

(D52), (D53), 
(D56), (D57) 

m 

59 

59 

10 

II 

It 

If 

0 

0 

(D52), (D62), 
(D64) 

15 

H 

II 

2.5 

Four (fig. 4c) 

PSTD* 

tl 

1 

1 

(E34) , (E35) 

16 

II 

II 

9 

n 

II 

II 

0 

0 

(E52)_(E56) 

17 

IV 

II 

5.5 


* PSTD stands for 


"pillow-shaped temperature distribution 


n 
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of the summation indexes employed in the assumed Fourier series for stress 
function and stiffener tensions; M pertains to the x-direction and N to 
y-direction. Only the highest M and N values used for each problem are 
shown. In general, calculations were also made for smaller M and N combi- 
nations in order that the convergence of the calculations could be observed. 
Column 10 shows the IBM 7074 computing time that was required to obtain, 
for the given M and N, all the results plotted in the respective figure. 

The stresses were computed at x/a and y/b intervals of 0.1. Because of 
symmetry it was sufficient to make calculations for only one-half or 
one-fourth of the structure. 
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Discussion of Numerical Results for the Thermal-Stress Problems 


Figures 6 through 13 , 16 and 17 present the computed results for the 
thermal stresses due to the pillow-shaped temperature distribution. The 
primary effect of assuming some of the stiffeners to be held straight is 
the creation of a running normal stress between these stiffeners and the 
adjacent edge of the plate. These stresses, which would be zero if the 
stiffeners were perfectly flexible, are depicted by N x (0,y) in figures 6a, 
7a, 8a, 9a, 10a, 11a, 12a, 13a, 16a, and 17a, and by N^(x,0) in figures 
16b and 17b. The maximum running normal stress between a rigid stiffener 
and adjacent plate edge is seen to range from about 35% to about 80% of the 
normal stresses produced in the middle of the plate. 

An interesting discontinuity occurs in the mutual normal stress 
between a plate edge and an adjacent perfectly flexible stiffener where 
that stiffener meets a rigid stiffener, e.g. the stress N^(x,0) at the 
point x = 0 in figure 6b and the point x = 0 in figure 10b. This stress 
is zero everywhere along the edge of the plate except at the corner where 
the rigid and perfectly flexible stiffeners meet. At that point the plate 
stress must jump to a value consistent with the strains at the edges of 
the two meeting stiffeners. The strain at the end of the rigid stiffener 
is zero, but the strain at the end of the flexible stiffener is not zero. 

It is the latter strain, in conjunction with Poisson's ratio, that gives 
rise to the non-zero plate stress right at the corner. 
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Comparison of the dashed curves with the corresponding solid curves 
shows that a further effect of holding one or more edges straight is to 
increase at least one of the normal stresses at the plate center. In 
figure 17a, for example, the normal stress at the plate center is seen to 
be increased almost 80% as a result rigidizing the stiffeners. In figure 
7a and 7b, it is seen that the effect of keeping the left side straight 
is to increase x-wise compressive stress and decrease the y-wise compres- 
sive stress at the plate center due to a pillow-shaped temperature rise. 

A similar result is shown in figures 11a and lib when both the left and right 
sides are forced to remain straight. 

On the other hand, comparison of dashed and solid curves for N 

xy 

shows that the maximum shear stress associated with the x and y directions 
is generally reduced by the existence of one or more straight edges. In 
the case of all four edges held straight, the reduction is seen to be quite 
drastic (figs. 16c and 17c). The corner shear stress, in particular, is 
reduced to zero for this case, as one should expect. The non-zero comer 
shear stresses in the other cases are a consequence of the zero bending 
stiffness for at least one of the two stiffeners meeting at every corner. 
Finite stiffener bending stiffness plus rigid joints at the corners would, 
in an actual situation, tend to reduce the corner shear stresses virtually 
to zero. 

In the case of perfectly flexible stiffeners (ref. 1), the stif- 
fener tensions approach zero at the ends of the stiffeners. In the present 
case, however, a finite end tension is produced in a perfectly flexible 
stiffener where it meets a rigid stiffener, or in each of two rigid stif- 
feners where they meet. This is shown in figures 6f, 7f, 8f, 9f, lOe, lie, 
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12e, 13e, 16d and 17d. However, the maximum tension in a stiffener still 
occurs away from the end. 

As is to be expected, the ratio of stiffener cross-sectional area 
to plate cross-sectional area has a noticeable influence on the plate normal 
stresses. Increasing this ratio (i.e. , diminishing X) tends to increase the 
plate normal stresses; this can be seen, for example, by comparing the central 
value of N in figures 6b and 8b . 

y 

Discussion of Numerical Results for the Applied-Force Problems 

In figures 14 and 15 are plotted the plate stresses and stiffener 
tensions produced by the pulling apart of two opposite stiffeners of infinite 
bending stiffness (fig. 4b with = T and all other loads and tempera- 

tures zero) . From figure 14a and e or 15a and e it is seen that this loading 
produces a nearly uniform tensile stress in the direction of the stretching 
(note that the scale in figure 14a does not start at zero). Deviations 
from uniformity are most pronounced near the ends (x = 0) . The plate stresses 
in the direction transverse to the stretching are seen to be quite small 
(figs. 14b and 15b) except near the ends, where the Poisson contraction of 
the plate is partially suppressed (fig. 14b) or completely suppressed 
(fig. 15b) by the end stiffener. 

For the case X -K), the plate cross-sectional area becomes negligible 
compared to the stiffener cross-sectional areas, and one should therefore 
expect that the applied load T would be transmitted from one side to the 
other entirely through x-wise stiffeners, causing P^(x)/T to be uniform at 
the value 0.5. Figure 15e shows that this predicted behavior was very nearly 
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achieved by the calculations, which led to values of 0.5 at x = 0 and 
.495 at x = 0.5a. The largest deviation from unif-orm stiffener tension 
occurred at x ~ 0.1a, where a value of .483 was obtained as a result of 
what appears to be a Gibbs phenomenon. For the tension P^(y) in the 
other two stiffener (fig. 15 d) , the calculation led to zero, which is 
value one should expect for this limiting case. 

The behavior of N in the neighborhood of the plate comer is 
interesting. It seems to become infinite for all values of X as the 
corner is approached (see figs. 14c and 15c). This singularity was not 
proved mathematically but the evidence for its existence is almost 
unmistakable from the computed values of N x ^(0,0) as the upper summation 
limits M and N in the assumed series are increased. This is shown in figure 
15f, where the computed dimensionless corner shear stress is plotted as 
a function of M on semi-log paper for the case M = N. The linearity of 
these graphs indicates rather convincingly that the comer shear stress 
would be infinite for the presumably exact solution (M = N ■> °°) . 

This behavior of the comer shear stress must be construed as a 
consequence of the assumption of perfect flexibility for one of the two 
stiffeners meeting at each corner. A similar behavior is obtained when all 
four stiffeners are perfectly flexible, for the case of a step-like tempera- 
ture discontinuity between stiffeners and plate or tensions applied to the 
stiffener ends, as is noted in references 1, 2 and 3. In an actual situation, 
finite flexural stiffness in the stiffeners which were here assumed to have 
none would tend to eliminate the shear stress singularity at the plate corners. 
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APPLICABILITY OF RESULTS 


There are many plane thermal-stress analyses in the literature for 
rectangular plates with free edges. In aircraft and spacecraft applica- 
tions, however, the plate elements that one is concerned with are almost 
always attached to stiffening members (e.g. , spar caps, rib caps, or 
shear-web uprights) along their edges. It was with such plate elements 
in mind that the work in references 1, 2, and 3 was done. In these 
references the stiffeners were assumed to have finite axial stiffness but 
negligible flexural stiffness and the boundary conditions were entirely 
those of prescribed load. Thus, while references 1, 2 and 3 represent 
a step toward more realistic detailed analysis of rectangular-plate 
plane-stress problems, they are not directly applicable to problems in 
which the boundary conditions are those of prescribed shape along one or 
more of the edges. The most obvious examples of such prescribed-shape 
boundary conditions arise when one considers plates that do not end at the 
stiffeners but are continuous across the stiffeners, forming a one-dimensional 
array of panels, as in a shear web, or a two-dimensional array, as in the 
cover of a multi-web multi-spar wing or the skin of a multi-ring multi- 
stringer fuselage. In these cases it may be reasonable to assume that an 
edge which is common to two adjacent panels is forced to remain straight. 

It is with such applications in mind that the present analysis and calcu- 
lations were performed. The assumption of perfect flexibility for the 
stiffeners was retained along some edges; along the other edges this 
assumption was replaced by one of prescribed stiffener shapes plus pre- 
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scribed resultants applied externally to the stiffeners of prescribed 
shape. Although the prescribed shape that one would use in most appli- 
cations is that of straightness, the prescribed shape was left arbitrary 
in the present analysis for the sake of generality and because such 
generality entailed very little additional complication. 

Figure 18 illustrates some of the types of structure to which the 
present paper and references 1, 2 and 3 may be applicable. Part a of 
figure 18 represents a single-bay structure made up of non-coplanar flat 
rectangular plates with stiffeners at their junctions and ends. The 
non-coplanarity of the sides suggests that there is little restraint 
against deformation of the plate edges. The assumption that any panel, 
such as 1234, is a plate with four perfectly flexible edge stiffeners 
may therefore be valid, and consequently the method of references 1, 2 and 
3 may be applicable for the detailed stress analysis of such a panel. 

Part b of figure 18 depicts a two-bay structure. If the loading 
and temperature distributions are approximately symmetrical about the 
central bulkhead or ring, then a panel such as 1234 may be considered as 
a plate bounded by the perfectly flexible stiffeners (12, 23, and 34) and 
one stiffener (41) that is forced to remain straight. Then appendix C 
and figures 4a, 6, 7, 8, and 9 of the present paper may be applicable. 

Part c of figure 18 denotes a multi-bay wing. In the shear web 
of such a wing, an interior panel, such as 1234, may be approximately like 
a rectangular plate with two opposite edges held straight (14 and 23) 
and with perfectly flexible stiffeners along the other two edges (12 and 
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34). Appendix D and figures 4b, 10, 11, 12, 13 of the present paper may 
therefore be applicable to such a panel. 

Part d of figure 18 represents the cover of a multi-spar multi-rib 
wing or a plate with integral waffle-like stiffening. An interior panel, 
such as 1234, may be approximately a plate with four edge held straight, 
and appendix E and figures 4c, 16, 17 of the present paper may be appli- 
cable to its analysis. 

In the above cases, wherever a stiffener is common to two adjacent 
plate, such as stiffener 12 in figure 18a or stiffener 14 in 18b, its 
cross sectional area should be divided in a reasonable way between the two 
plates which it serves. For example, if in figures 18a, b and d the adjacent 
plates are nearly identical in geometry, loading, and temperature distri- 
bution, then it is reasonable to divide each shared stiffener equally between 
the two plates which it serves. The same can be said about stiffeners 14 
and 23 in figure 18c if the adjacent panels of the shear are similar in 
geometry, heating, and loading. However, the two plates which share 
stiffener 34 in figure 18c (one in the shear web and one in the top cover) 
could be quite different. In such a case it might not be easy to decide 
on the effective cross-sectional area to be used for stiffener 34 when 
analyzing the plate 1234. When such ambiguities result in a significant 
uncertainty in the computed results, then an analysis of the present kind, 
based on isolating one plate and its surrounding stiffeners as a unit, 
might not be applicable. In the unlikely event that the cover plates are 
missing in figure 18c, or have negligible cross-sectional area, then the 
ambiguity referred to above would, of course, not exist. 


35 


The above remarks concern the applicability of the present analysis 
in general. Some brief comments are, perhaps, also in order on the appli- 
cability of the computed thermal-stress data plotted in figures 6 to 13, 

16, and 17. The geometry and temperature distributions to be expected in 
actual design situations are, of course, manifold. It is therefore not 
suggested that the geometries and temperature distributions assumed in the 
thermal-stress calculations will correspond exactly to a broad spectrum 
of particular design problems. Nevertheless, it is felt that the calcu- 
lations made and the results obtained serve two purposes. First, they attest 
to the feasibility of the calculation procedure, thereby encouraging its 
use for other specific geometries and temperature distributions which might 
correspond more closely to a particular application. Second, they provide 
qualitative information on the thermal plane stresses that may be encountered 
during aerodynamic (or other) heating as a result of a temperature lag between 
a thin plate and the more massive or better insulated stiffening elements 
to which it is attached. The pillow-shaped temperature distribution 
0 sin (tx/a) sin (iry/a) employed in the calculations is not too much dif- 
ferent from the experimental temperatures obtained in a square plate that 
was heated by lamps and cooled by water circulating through hollow edge 
stiffeners (fig. 11 of ref. 4), although the experimental distribution is 
somewhat flatter in the central region of the plate and might be better 
represented by the temperature distributions used in reference 3. In the 
latter reference the temperatures were assumed constant in a central plateau 
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region of the plate, zero at the stiffeners, and varying as one-quarter 
of a sine wave in the transition zone between the stiffeners and the 
central plateau. 

In the thermal-stress calculations (fig. 6 to 13, 16, and 17), 
the rigid stiffeners were assumed to be entirely free of external re- 
straint. That is, the external resultants, such as and in figure 
4b, were assumed to be zero along with all the other force loadings. 

Thus, although certain edges were forced to remain straight, these edges 
were free to translate. In an actual situation, there might exist some 
restraint against this translation due to the presence of surrounding 
structure. The effect of such restraint would be to introduce non-zero 
values of the external resultants T^, T£ , etc., with magnitudes depending 
upon the relative stiffness of the plate-stiffener combination and the 
surrounding structure. The stresses due to such externally developed 
restraining forces would have to be superimposed as corrections upon 
those shown in figures 6 to 13, 16 and 17. Figures 14 and 15 show the 
nature of these correction stresses for the case in which restraining 
forces T^ = T^ = T are developed in the structure of figure 4b. From 
figures 14a and 15a, it can be seen that one effect of preventing the 
separation of the end stiffeners (through the development of a negative 
value of T) is to superimpose a nearly uniform compressive stress in the 
x-direction. From figures 14c and 15c, it is seen that another effect 
is, possibly, the introduction of high plate shear stresses near the 
comers. However, as mentioned earlier, the shear stress singularity 
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at the plate corners would be eliminated in an actual structure by the 
presence of some flexural stiffness in the x-wise stiffeners. 


To illustrate the use of data of the type shown in figures 14 
and 15, let us consider the thermal-stress problem whose results are 
given in figure 11, but this time assume that constraining forces 
= T£ = T exist, which are of appropriate (negative) magnitude to 
prevent completely the overall thermal expansion of the structure in 
the x-direction. In order to determine T, it is necessary to determine 
the elongation A^ due to thermal stress from the data of figure 11, the 
elongation A^ due to T from the data of figure 14, and equate their sum 
to zero and solve the resulting equation for T. The two required elon- 
gations can be obtained by integration of the stiffener strains e^(x) 
or of the plate strains e x (x,y) along any line parallel to the x-axis. 
The stiffener strains are slightly more convenient in this case, because 
of the absence of the thermal contribution e^Cx); they give rise to the 
following expressions for A^ and A^: 


A i ' 


A 2 ' 


P 3 (x) 

~AE 


dx = 20Xota 


P 3 (x) „ Ta 

1 —** = 2 AE 


1/2 

' 

0 

1/2 

. 

0 


P Q (x) 

-.3 

QAAEa u V 


P 3 (x> x 

-T— d < 7 > 


Using numerical integration (trapezoidal rule) in conjunction with figures 
lie and 14e to evaluate the right-hand sides of these expressions, one 
obtains 
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A x = 0.2393 eActa 


A 2 = 0.2185 Ta/(AE) 

The condition A^ 4- A 2 = 0 then gives 

T = - .4439 eaahE 

The plate and stiffener stresses due to this value of T can now 
be obtained from figure 14. Superposition of these stresses and the 
ones in figure 11 will give the stresses produced by the pillow-shaped 
temperature rise when the overall expansion of the plate in the x-direc- 
tion is prevented. The results are shown by the dotted curve of figure 
19. The solid curves are the data of figure 11, repeated for comparison. 

It is seen from figures 19a and e that the normal stresses in the x-direc- 
tion are, of course, appreciably altered by the suppression of the overall 
thermal expansion in the x-direction. Figures 19b and d show only minor 
alterations in the normal stresses in the y-direction. Figure 19c shows 
that the shear stresses also are only slightly affected, except in the 
neighborhood of the comer (x and y < 0.1a), where the effect of the singu- 
larity in shear stress due to T is felt. Accepting this singularity as 
illusory for reasons already cited, it can be concluded that in this 
example the suppression of the x-wise expansion has a significant effect 
only on the x-wise normal stresses. 


CONCLUDING REMAEKS 


Many detailed plane-stress analyses are available in the literature 
for rectangular plates with unstiffened edges, subjected to boundary 
forces or non-uniform temperature distributions. Rectangular plates with 
stiffened edges have been treated to a lesser extent in the literature, 
mainly in "shear-lag" analyses. These are generally restricted to uni- 
form temperatures and unidirectional loading and are characterized by 
simplifying assumptions regarding the stresses or elastic constants. 

Rectangular plates in practice generally have stiffening members 
along their edges and may be subjected to non-uniform temperature rises 
and multidirectional loading. In addition it may be important to have 
some of the detailed information about the stresses which is lost in the 
usual type of shear-lag analysis. Therefore in reference 1 a unified 
analysis of the edge-stiffened rectangular plate was presented, including 
both non-uniform temperature distributions and prescribed boundary load- 
ings of a very general nature. The plate material was assumed to be 
homogeneous, linearly elastic, and orthotropic, but no further simpli- 
fying assumptions were introduced regarding the stress distributions 
or the elastic properties of the plate. The stiffeners were assumed to 
have finite extensional stiffness but negligible flexural stiffness. The 
boundary conditions considered were entirely those of prescribed running 
normal and shear loading along the outer periphery of the stiffeners and 
prescribed forces at the stiffener ends. 
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The present paper represents a generalization of reference 1. 

It permits the boundary condition of prescribed normal loading along 
some of the edges to be replaced by the boundary condition of pre- 
scribed shape of the stiffener axis. Generally, the prescribed shape 
that one would encounter in practice is that of straightness due to 
symmetry, continuity, or rigid fixation. However, for the sake of 
generality, and because it involved no great complication, the present 
analysis also permits prescribed edge shapes other than straight. Along 
the remaining edges the stiffeners have been assumed to have finite ex- 
tensional stiffness but negligible bending stiffness, as in reference 1. 
Through this generalization of the boundary conditions, certain problems 
can be handled which are not directly solvable by the method of reference 
1 alone. 

The present method and that of reference 1 are based on Fourier 
series - double series for the plate stress function, single series for 
the stiffener tensions and certain other quantities. In some cases many 
terms of the series may be needed for sufficient accuracy, and therefore 
the feasibility of these methods depends on the availability of high-speed 
computers. This feasibility has been tested in the present paper through 
twelve specific numerical examples, solved with the aid of an IBM 7074 
computer. The machine time required to obtain detailed stress surveys in 
one quadrant or in one half of the structure, depending on the type of 
symmetry, ranged from 2.5 minutes to 21 minutes in these examples. 
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The assumption of negligible bending stiffness for those stiffeners 
whose shape is not prescribed sometimes produces infinite plate shear 
stress at a corner where such a stiffener meets a stiffener of prescribed 
shape. This was noted in figures 14 and 15 of the present paper. A 
similar anomaly was observed in references 1, 2 and 3. Its elimination 
would require the incorporation of finite stiffener flexural stiffness 
into the analysis . 



APPENDIX A 


SYMBOLS 


Remarks: (i) The subscript 1, 2, 3 or k on a symbol for a stiffener- 

related quantity identifies the stiffener location as x = 0, x = a, y = 0, 
or y = b, respectively. (ii) The Fourier Coefficients of known quantities 
(loads, thermal strains), and combinations of such coefficients, are 
generally represented by capital letters, while the Fourier coefficients 
of initially unknown quantities (e.g., internal stresses) are denoted by 
small letters. (iii) The symbols for convenience represent the known 
quantities in illustrative problems are not listed here. 


a 

a 

mn 


« t it _ * it _ mi 

8. J 3. , ct y a. 

n rr m m 


plate dimension in x direction; see figure 1. 

Fourier coefficients in series expansion for the stress 
function F(x,y); see equations (B15), (B17) and (B50). 
Fourier coefficients in series expansions for F(0,y), 
F(a,y), F(x,0), F (x,b) respectively; see equations 
(B13) and (Bl4). 


A ij A 2 jA 3' A 4 

A 


b 


B ' , B",B 1 
n ’ n* " 


B"" 

m 


B 


c 

mn 


stiffener cross-sectional areas. 

common value of the above when all four stiffeners are 
identical and uniform. 

plate dimension in y direction; see figure 1. 

Fourier coefficients in series expansions for N-^N^N^N^ 

respectively; see equations (4) and (5). 

a/b 

Fourier coefficients in series expansion for N^(x, y); see 
equations (B22) and (B35)* 
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_ r „ it 

c , C 
n 7 n 

Fourier coefficients in series expansions for N (0,y) 
N (a,y) respectively; see equations (B23) and (B24). 

y 

c ' , c " 
n 7 n 

c' C 0 mr/b, c" C_ nn/b. 
n 2 ' 7 n 2 ' 

C 1' C 2 7C V C 4 

plate compliances defined by equations (3). 

C 

A 5 /A 1 , if A x = A 2 and = A^. 

d 

mn 

Fourier coefficients in series expansion for crF/dx , 
see equations (B4l) and (B46) . 

D 

(7 < 2) ) 2 - 

n 

n n 

D' 

n 

defined by equation (Dll) . 

D" 

n 

defined by equation (D12). 

/ 

e 

mn 

Fourier coefficients in series expansion for c^F/cbc^; 
see equations (B 26 ) and (B37)* 


^ 1 (y)^e 2 (y),e 5 (x),e J+ (x) 

stiffener thermal strains; see figure 2. 


e x (x,y),e y (x,y) 

plate thermal strains; see figure 2. 

(^(mjt/a) 1 * + (C^ - 20^) (mjt/a) 2 (nn/b ) 2 + C^nn/bjV 
value of above with m = 1 and n = 1. 

E^E 2 ,E^,Ej^ Young's modulus for stiffeners. 

E Young's modulus for stiffeners and isotropic plate when 

all have the same Young's modulus. 

F(x,y) stress function for plate; see equation (B4). 

g^ Fourier coefficients in series expansion for N^(x,y); 


®m 7 ®m 


see equations (B19) and (B34). 

Fourier coefficients in series expansions for N x (x,0) 
and N x (xjb) respectively; see equations (B20) and (B21) . 
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i bO 


m* ®m 


mn 


H 


mn 


mn 


'mn 


C L nut/a, g^ C L nut/ a 
thickness when plate is isotropic. 

3 3 

Fourier coefficients in series expansions for crF/djr ; 
see equation (b 42) and (b47). 

(l/E mn )(2/a)(nVb) 5 C 1 - (2/a)(b/mt) 

1, 2, 3, or 4. 

Fourier coefficients in series expansion c^F/By^; 
see equations (B27) and (B38). 

Fourier coefficients in series expansion for If ; see 
equation (B25); see equation (B 56 ) for value of j 


K 


mn 


K',K",K'",K"" 
n’ n' m ' m 


combinations of known Fourier coefficients, defined 
by equation (B 69 ) . 

Fourier coefficients in series expansions for prescribed 
boundary curvatures; see equations (8) to (12). 

3 2 

Fourier coefficients in series expansion for (iuF/8x 8y); 
see equations (B43) and (b 48). 
summation indexes (integers), 
upper limit on m, p, and i . 
loading resultants; see figures 4 . 
summation index (integer) . 
upper limit on n and q. 

N 1 (y),N 2 (y),N 3 (x),N 4 (x) 

external running tensions, force per unit length; see figure 1. 
plate stress-resultants, force per unit length; see figure 3« 
summation index (integer). 


mn 

m,n,p,q 

M 

m 1 ,m 2 ,m 5 ,m 4 

n 

N 


N ,N ,N 
x 7 y xy 


mn 


4/22 

Fourier coefficients in series expansion for 8 F/dx by ; 
see equations (B 28 ) sind (B39)* 
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P 1 (y),P 2 ( y ).i > 5 (x), p ^(x) 

stiffener cross-sectional tensions; see figure 3* 
P 1 (0'),P 1 (b),P 2 (0),P 2 (b),P 5 (0),P 3 ( a ),P 4 (0),P 4 (a) 

stiffener end loads; see figure 1. 

P, Q integers appearing: in sinusoidal temperature distribution 

expansion; see equation (C7l) . 

q summation index (integer). 

q 1 (y)^q 2 (y)>q 5 (x),q 4 (x) 

external shear-flow loadings; see figure 1. 

Q^, Q^, Q^", Fourier coefficients in series expansions for q^, q^q 

respectively; see equations (6) and (7)* 

R' ,R",R'",R ,m combinations of known Fourier coefficients, defined by 
n 7 n 7 m 7 m 7 J 

equations (B 68 ) and (B 69 ) . 


s',s",s' ,s Fourier coefficients in series expansions for the stiffener 
n 7 n m m 

cross sectional tensions; see equations (Bl6) and (b6o). 

S (D s (2) s (3) s 00 

n 7 n 7 m 7 m 

completely known loading terms; defined by equations (C17) 

through (C20). 

s (5) jS (6) (7) >s (8) 

n 7 n 7 m 7 m 

completely known loading terms; defined by equations (Dl8) 
through (D21). 

^n 7 ^n 7 ^m" 7 ^m" cc^pletely known loading terms; defined by equations (E15) 
through (El8) . 

known quantities; defined by equation (C21) . 
known quantities; defined by equation (D22) . 

known loading terms, defined by equations (C58b) and (C59a) • 


s d) 

mn 

s< 2 > 

mn 

s ( 5 ) ,s (4) 

m m 
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t',t",t , ",t" ,, Fourier coefficients in series expansions for the derivatives 
n’ n* m ' m ^ 

of the stiffener cross-sectional tensions; see equations 


T 

ran 


(B29) and (B33) • 


Fourier coefficients in series expansion for 
2 2 

d e x /c)y ; see equation (l5)> also equations 
and (l8) . 


S 2 e y/Sx 2 + 

(16), (17) 


T' ,T" ,T' ", 
n* n* m * 


iplltt 

m 


T i>V T y T i. 


T 

o 


u,v 

U 

n 


mp 

V’ 

mp 

V" 

mn 

V' 

n' n m 7 m 


w 


mn 


x,y 

y 

y',y' 


Y , ,Y",Y ,M ,Y"" 
n 7 n 7 m 7 m 


Fourier coefficients in series expansions for thermal- 
strain discontinuities between stiffeners and plate edges; 
see equations (13) and ( 14 ) . 
loading resultants; see figure 4. 

stiffener temperature for the sinusoidal plate temperature 
distribution. 

x and y components of displacements in plate, 
known quantities, defined by equation (C57a) • 
plate displacement component in y- direction, 
known quantities, defined by equation (C58a) . 
known quantities, defined by equation (D58a) . 
known quantities, defined by equation (D58 c). 

Fourier coefficients in series expansions for (dey/dx) x _Q 
etc.; see equations ( 19 ) through ( 23 ). 

3 2 

Fourier coefficients in series expansion for (d^F/chc^y ); 
see equations (B44) and (B49). 

Cartesian coordinates; see figure 1 
Cartesian coordinate; see figure 1 
dummy variables representing y 

known loading terms; see equations (C70a) , (c 66 a) , (c67a) , 
and (C68a) . 
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Z^,Z^ known quantities; defined by equations (D64a) and (D5&b). 

a coefficient of thermal expansion of plate and stiffeners 

in numerical examples. 

a 1 (n),a 2 (n),a 3 (m),a i+ (m) 

known quantities; defined by equations (b 67 ). 

P^(n),Pg(m) known quantities, defined by equations (b 67 ). 

known quantities, defined by equations (C22), (C23), (D23) 
and (C56d) . 


7 (D y(2) 

'n * ' n 

r , ,r , ,r"' 

m m m 

,(1) r (2) 


m 


m 


mp 7 mn 7 mn 

r (2) r (3) r (M 

) Tyiyi ^ 


mp 


mn 7 mn 


& (!) 6 ( 2 ) 

n 7 n 

8 ', 5 ", 6 '", 8 "" 

n 7 n m m 


known quantities, 
known quantities, 
known quantities, 
known quantities, 


given by equations (CIO) and (d 6 ) . 
given by equations (E 19 ) through (E2l) . 
given by equations (E25a) and /E25b). 
given by equation (C 36 ). 


known quantities, given by equation (D. 35 )- 
known loading terms, defined by equations (Cll) . 
known loading terms, defined by equations (D7) and j[D 8 ). 
known loading terms, defined by equations (E23a), (E2Ua), 


(E25c) and (E26a). 


8 . . 
ID 


A 

n 

A', A" A'" 
n 7 n 7 mn 


Kronecker's delta, unity when both subscripts are equal, 
zero otherwise. 

completely known loading terms, defined by equation (c64a). 
known quantities, defined by equations (C56a), (C56b) and 


(C56c) . 


n 7 pn 
e x (x,y),e (x, 

e i(y)> e g(y), e 


known quantities, defined by equations (D 6 lb) and (D63a) . 

y),7 xy (*,y) 

plate total strains; see equations ( 2 ). 

3 (x),e ]+ (x) 

stiffener total strains; see equations (l) 


i iil 
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^n 

known 

quantities 

^n 

known 

quantities 

t (l) t ( 2 ) 

known 

quantities 

4' > T l" 
'mn' 'mn 

known 

quantities 

(l) (2) 

r Wi ’ Vn 

known 

quantities 


0 ’ , 0 " , 0 '",©"" 

rap mp mp mp 


known quantities, given by equations (C40), (C42). 


©(i) (2) (3) ©(4) 

mp 7 mp mp mp 


0 ' , 9 " 
mp 7 mn 

V X 2 


4 


A' , A" 
mn’ mn 

p' ,p" 

mn’ mn 


v ' , v 
mn mn 


known quantities, given by equations (D 38 ) and (D4o). 
known quantities, defined by equations (D57a) and (D57b). 
area-ratio parameter used in numerical example, defined 
by equations (C 80 ) 

common value of 1 ^ and 4^ when they are equal, 
known quantities; defined by equations (E23b) and (E25d) . 
known quantities, defined by equations (C12) and (|E25e). 
known quantities, defined by equations (C13) and (E25f). 
Poisson's ratio for isotropic plate. 


t ' p" 

hrm’ b mn 

known 

quantities, 

defined by 

equations 

(C35). 

.U) ,( 2 ) 

b mn ’ 5 mn 

known 

quantities, 

defined by 

equations 

(D34). 

^pn 

known 

quantities, 

defined by 

equations 

(D 56 c) . 

_ n 

known 

quantities. 

defined by 

equations 

(C56e). 

— 1 — " 

- n’ ~n 

known 

quantities, 

defined by 

equations 

(C 56 f) and 

n 

mn 

known 

quantities. 

defined by 

equations 

(E22) . 

P', p" 
nr m 

known 

quantities. 

■defined by 

equations 

(C4 3 ) . 

('l) (2) 

> P t ' 

m m 

known 

quantities. 

defined by 

equations 

(D4l). 
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- ! t »t -Mil 

0 » 0 • 0 
mp 7 mrr mn 


n 


known quantities, 
known quantities. 


^n’^n’^a ’^n' kno " n quantities. 


0 (D, 0 <a) >0 O) >0 OO 

r n ’ r n r n r n 


n " n 

known quantities. 


defined by equations 
defined by equations 
defined by equations 


defined by equations 


(C37) • 
(D 56 a) . 

(C53). 


(D32). 


$ 


n 


n 


d>. 


( 1 ) 


n 


T mn 

■Jr ' 
nnn 


% 


mnp 


ijr ' 

mnp 


Cf’Cf nn - 0"cf"\ 
r n r n r rv n 

jjd)/*) . 0 (s) 0 (3) 

r n r n r n r n 


1 - E" 

n 

known quantities, defined by equations (C48) . 
known quantities, defined by equations (D5l)« 
known quantities, defined by equations (c47). 

known quantities, defined by equations (D50). 
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APPENDIX B 


ANALYSIS FOR THE CASE IN WHICH THE BOUNDARY CONDITIONS 
ARE ENTIRELY THOSE OF PRESCRIBED LOADING (ref. l) 


The case in which curvatures are prescribed along one or more edges 
can best be analyzed by making appropriate modifications in the basic 
analysis of reference 1, in which the boundary conditions are entirely 
those of prescribed loading. Therefore the analysis of reference 1 
for the case of constant- area stiffeners is summarized in this appendix, 
and portions of it will be used as needed in the subsequent appendices 
when various cases of prescribed boundary curvature are considered. 


Basic Equations 

With u(x, y) and v(x, y) denoting the x- and y-components, respectively, 
of infinitesimal displacement, the strain- displacement relations for the 
plate are 


Bu Bv Bu | Bv 

€ x ~ Bx £ y By ^xy - by bx 


(Bl) 


Equstions (Bl) imply the following compatibility condition on the strains 

(B2) 




bxby 


2SZ. 


b 2 e 


b 2 e 


JL - 


by 2 bx 2 


= 0 


The plate equilibrium equations, namely 


Bn BN 

x + _ o 

Bx By " ° 


Bn Bn 

x y + __y 

dx by 


(B3) 


imply the existence of a stress function F(x, y) such that 
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N 

x 


d 2 F 

Sy 2 



N = 

xy 


-&~F 

dxby 


(B4) 


Eliminating e , e and 7 in equation (B 2 ) by means of the stress- 
x y xy 

strain relations, equations (2), and then N^, and N through equations 
(BU), yields the following form of the compatibility equation: 


, c>V 
J 2 7 T 


+ (C), - 2C ) 


aV 

&V 1 + ax 2 


dV 


a 2 e 


= 0 (B 5 ) 


ax" - ^ ax^ay" x a y " a x " a y ‘ 

Considering infinitesimal lengths of the stiffeners as free bodies, 

and utilizing the third of equations (b 4 ) to express N at the plate 

xy 

edges in terms of F(x, y), one obtains the following equilibrium equations 
governing the longitudinal variations of the stiffener axial forces . 


dP L /dy - (a 2 F/a X ay) x=0 - q L (y) = 0 
dB 2 /dy + (a 2 F/a X ay) x=a + q^y) = 0 
dP^/dx - (a 2 F/a X a y ) y=0 - q^(x) = 0 
d^/dx + (a 2 F/a X ay) y=b + q^(x) - 0 


> (B6) 


The stiffeners and plate edges, being integrally attached, must have 
equal longitudinal strains along all four edges. Using the stress- strain 
relations (eqs. (l) and (2)), the stress function defined by equations 
(B 4 ), and the assumption of perfect stiffener flexibility (i.e., 

(N x )x_ 0 = N (y), etc.), these conditions lead to the following 
equations: 
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P 1 (y) b 2 ? 

+ E e i(y) - e(o,y)] - C F-|) + c N (y) = 0 

A 1 E 1 1 y 2 ax 2 x=0 5 1 


P p (y) 

Ag ~ Eg ' + t e 2 (y) - e y (a,y)] - C 2 (^-|) + c 5 N 2 (y) = 0 


■ b^F - 

k r 

c)x x=a 


> (BT) 


A^r + [e 3 (x) - e x (x ' 0)] - c i<r-f) + c 3 n j ( x> ■ 0 

3 3 oy y=0 


P (x) 2 

A- + [e 4 (*> - e x (x,b)] - C (f|) + 0 n 4 (x) = 0 

4 4 dy y=b 


Boundary Values of Stress Function F(x,y) 

The fact that the distributed boundary normal loadings N (y), 

Ng(y), N^(x), and N^(x) are transmitted directly to the plate means 

that the second derivative of F(x,y) in the direction of the edge is 

known. Therefore two integrations will give the variation of F along 

each edge in terms of the unknown corner values and the known N^, N , 

2 2 

N^, N^. As an illustration, along the edge x = 0, d F(0, y)/Sy =■ W 1 (y), 
and two integrations lead to 


b y' 


F(0,y) = F(0,0) + J [F(0,b) - F(0,0) - 


N 1 (y ,, )dy"dy l ] + 


0 0 


y y' 

+JJ N 1 (y")dy ,, dy- (B8) 

0 0 


(Here F has been assumed continuous along the perimeter of the plate, 
including the corners.) Replacing N^(y") by its series expansion, 
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equation (4), and carrying out the integrations indicated in equation 
(B8), one obtains 

1M 

2 


F(0,y) = F(0,0) + l [F(0,b) - F(0,0)] - £ sin ™ 

n=l 

(B9) 

Similarly, the variation of F(x, y) along the other plate edges is as 
follows: jj. 

F(a,y) - F(a,0) + *[F(a,b) - F(a,0)] B" (^) sin 


n=l 


(BIO) 


M 


F(x,0) = F(0,0) + |[F(a,0)-F(0,0)] B^"(^) sin 


mirx 


m=l 

M 


(Bll) 


F(x,b) = F(0,b) + ^[F(a,b)-F(0,b)] sin ^ 

m=1 (B12 ) 

For the later use it will be necessary to expand the boundary 

values of F in single Fourier series of the following form; 

N 

F(0,y) = ^ a n sin Cn*y/b ) (0 < y < b) ^ 

n=l 
N 


F(a,y) = ^ a" sin (nrty/b) (0 < y < b) 

n=l 

M 

F(x,0) = I a^” sin (mitx/a) (0 < x < a) 

m=l 

M 

F(x,b) = ^ a^" sin (mnx/a) (0 < x < a) 


> (B13) 


m=l 


Determining the coefficients in these series through the formula 


etc. 


b 

a n = ( 2 / b ) / F (°>y) sin ( nn y/ b ) 

o 

, with F(0,y), etc. replaced by the right-hand sides of equations 
(B9) to (B12) y one obtains 

2 

a ' = ~~ [F(0,0) - (-l) n F(0,b)] - (— ) B’ 
n nit \ / \ n:Jt / n . 

> (bi4) 


Series Assumptions for F(x,y) and P - ^y), P c (y), P^(x), and P ^ fx) 

The stress function F(x,y) will be assumed to be representable in 
the interior of the plate (i.e. in the open region 0 < x < a, 0<y<b) 
by the double Fourier series 
M N 

F(x,y) = II a^ sin (nmx/a) sin (nny/b) (B15) 

m=l n=l 

with as yet unknown coefficients. Equation (B15) is not valid at the 
plate edges; however there the values of F are already expressed in 
series form by equations (B1J) and (Bl4) . Similarly the unknown stiffener 
axial forces will be assumed in the form 


- (-l> n lHa,b)] - &) B” 

o 2 

a'" = [F(0,0) - (-l) m F(a,0 ) ] - (— ) B 1 " 

m mn ' ’ ' v mjr m 


a"" = — [f(o b ) - (-l) m F(a,b ) ] - (— ) B"" 

m mjt v v K ’ ‘ J m 
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P x (y) 

P 2 (y) 

p 5 (x) 

P 4 (x) 

At the end cross sections the stiffener forces are already known from 
the given loading (see fig. l) . 

The coefficients in the series in equations (B15) and (Bl6) are 
related to the left-hand sides through the usual formulas 


b a 


ab J J P(x,y) sin a sm dx dy 

(B17) 

0 0 


b 


i f P i^ y ^ sin ^ dy, etc. 

(B18) 


0 


N 


= ^ s^ sin (njty/b) (0 < y < b) 


n=l 

N 


E s" sin (nrty/b) 
n 

n=l 

M 

^ s^" sin (irntx/a) 

m=l 

M 

Z s"" sin (mitx/a) 
m ‘ 


(0 < y < b) 
(0 < x < a) 
(0 < x < a) 


> (B16) 


m=l 


Series for the Derivatives of F(x,y) an d P ^, P , 

The derivatives appearing in equations (b4) to (B7) will be assumed 
representable by the following series: 


M N 

= a 2 F/5y 2 sin ( m W a ) 

m=l n=l 


sin (rnty/b) 


(0 < x < a) 
(0 < y < b) 


(B19) 
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I 


( W xV=0 = (d 2 F/By 2 ) y=Q = X g m Sin ( mitx / a ) (0 < x < a) 

m=l 

M 

( N x ) y =b = ^ & 2 F/Sy 2 ) _b = X ®ta S±n ("Wa) (0 < x < a) 


M N 


N y = B 2 F/Bx‘ ! = IE sin (nntx/a) sin (nity/b) 


m=l n=l 


(0 < x < a) 


(0 < y < b) 


( N y ) x =0 = (d^/Bx 2 )^ = X C A Sin ( n5Ty / b ) (0 < y < b) 


^ N y^x=a = ^ 2F /^ x2 ) x =a = X C n sin ( niT y/ b ) (0 < y < b) 

n=l ( B24 ) 


M N 

= -B 2 F/BxBy = ^ cos (mirx/a) cos (nny/b) 

m=0 n=0 


(B25) 

M N (0 < x < a) 

cTF/Sx ‘-I i sin (mjtx/a) sin (mty/b) (B 26 ) 

m=l n=l (0 < y < b) 


MW (0 < x < a) 

4 4 ' V - * 

s F/By = 2^ ) i^ sin (imtx/a) sin (nity/b) (B27) 

m=l n=l (0 < y < b) 


ra=l n=l 


Ij. 2 2 V 1 V - * 

B F/Bx By = \ sin (mirx/a) sin (mty/b ) 


M N 


(0 < x < a) 


ra=l n=l 


(0 < y < b) 


(B28) 



N 


dP-j/dy = I V cos (nny/b) (° < y f: b ) 


n=0 

N 


dP 2 /dy = y~* cos (nny/b ) (0 < y < b) 


n=0 

M 


(B29) 


dPy'dx = ^ t^" cos (mnx/a) (0 < x < a) 


m=0 

M 


dP^/dx = I t^" cos (mnx/a) (0 < x < a) 


where 


m=0 


b a 


(d^F/dy 2 ) sin (mnx/a) sin (nny/b) dx dy, etc. 

(B30) 


0 0 


a 

gjj = | J (S 2 F/a y 2 ) y=0 sin (imtx/a) dx, etc. 


(B31) 


J ran 


( g - 6 mo» 2 - 5 n 0 > 

ab 


b a 


f J~ (S 2 F/Sx5y)cos(nmx/a)cos(nny/b)dxdy 

(B32) 


0 0 


2-6 


t' = 


nO 


n b 


D 

f ( dP x / dy ) 


cos (nny/b) dy, etc. 


(B33) 


The coefficients appearing in the series for the derivatives (eqs. 

(B19) to (B29)) are not independent of the coefficients in the series for 

the basic quantities (eqs. (B15) and (Bl6)). The former can be related 

to the latter by using integration-by-parts in the right-hand sides of 

equations (B30) to (B33)- For example, two partial integrations with 

Such a technique was employed for plate bending problems by A. E. Green 
(ref. 5), who ascribes its earlier use to S. Goldstein. 
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respect to y in equation (BJO) give 


g = Sjf g [ a .» - (_i) n a ,M ’] 

mn b b m v m 


f n jtx 
r J a 
b mn 


(B34) 


Similarly, 


c =22£ [a . 
mn a a n 


(-D m a^] 


(— ) a 

'a mn 


(B35) 


^ — [(-l)"^ 11 F(a,b) - (-l) m F(a,0) - (-1)^(0^) + F(0,0)J 


'mn 


+ fit f^nO [( . 1) n a „„ 
a b ' ' m 


m 


3 + ^ 


2 - & n 

nn m0 


b a 


[(-D^a” - 


a ' ] + 
n 


, mn nn 
+ — — a 

a b mn 


(B36) 


mn 2 r , 

e = Lc ' 

mn a a n 


3 0 U 

(_l) m c "] _ (™) 2 [a , _ ( . l) m ,,j + (M) a 

v ' n v a' a n ' n a' mn 


(B37) 


nn 2 
L mn b b 




(-i)XJ 


O 3 I 

b b m 


(-l) n a""] + (2g) a 
v / m ' b' mn 


(B38) 


= 3 T ^ [(-l) mfn F(a,b ) - (-l) m F(a,0) - (-l) n F(0,b) + F(0,0)] 


■ mn ab a b 


+ | (— ) ~ L(-l) n a"" 

b v a b ' ' m 


a J + ~(~r) — [(-l) m a" - a-J 
m a^b 7 a v n n 


2 2 
, /mn, mn, 

+ <■— > C-b> a mn 


(B39) 
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*a - ^ [(- 1 ) 11 V b ) - v °) ] + “ 


2-5 


t" = . 

n b 


no [(-l) n P 2 (b) - P 2 (0)]+^s^ 


tc- 1 )" V a) - V 0)] + H ?C 


> (BUo) 


2-5 


m 


f 9 - [(-D m P 4 (a) - P 4 (0)] + ^ s™ 


J 


Reference 1 did not require series expansions for cpF/dx"^, c)^F/dy^, 
d 5 F/dx 2 dy, and c^F/dxdy 2 ,, but the present paper will require such expansions 
in the subsequent appendices. These derivatives will therefore be 
assumed expandable in the following series: 

M N (0 < x < a) 


d^F/dx^ = / / d cos (mjtx/a) sin (mry/b) 

4 « • » mn 


m=0 n=l 


M N 


S^F/Sy^ = y / h sin (mrtx/a) cos (njty/b) 
i * » mn 


m-1 n=0 


(0 < y < b) 

(B4l) 

(0 < x < a) 

(0 < y < b) 

(B42) 


M N 


d^F/dx 2 dy =/ ) ^ sin (mnx/a) cos (mty/b) 

L * / ■ mn 


m=l n=0 


M N 


(0 < x < a) 

(0 < y < b) 
(B43) 


, 2 v — 1 r - " 1 , (0 < x < a) 

d^F/dxdy = > ) w cos (mjtx/a) sin (mty/b) 

n i “H (o < y < b) 

m=0 n=l 

(B44) 


where 


(2-6 ) p 

^ mo' 2 


b a 


mn 


^ J J (^F/Sx^) cos (mrcx/a) sin (mry/b) dxdy, etc. 

(B45) 


0 0 
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By using integration-by- parts in the right-hand sides of above equations, 
one obtains 


d = [(-l) m c" - c'] - -(— ) 2 [(-l) m a" 

mn a ' n n a a ' n 


a A ] 


/nut-> 

a 


mn 


Similarly, 


h 

mn 


U-D\ - a;] - |(f) 2 [(-i)V" 


a’"] 

m 


(B46) 


{.—■ rv a 
v b mn 


(B47) 


mn 


nut . 
a J mn 


2*~5 

2 no /nut \ r / s m+n , \ 

F(a,b) 


- (-!)“ F(a,0 ) - (-l) n F(0,b) + F (0,0 ) ] 


w = 
mn 


2-5 2 

no ,nut s r f , mi 

b ( ~V [( " 1) a m 


mt 


[ (— l) m a"— a ’ ] - (— ) a 

m a a b ' n n ^ a bmn 


(B48) 


b J mn 


(2-5 ) 

* m r~\ ' 


a f(-f) [ (- 1 ) F(a,b) - (-l) u F(a,0) - (-1)“ F(0,b) + F(0,0)] 


_ nut g nn ^ ^ n [(.qjV - a'] 

a b b ' 7 m m a b' n n 


nut / nit ^ 
a b a mn 

(b49) 

Through equations (B34) to (B4o) and (b 46) to (b 49) all the unknown 
coefficients in the derivative series are expressed in terms of the basic 
unknowns a : c’, c", g', g": s', s", s'", s"": F(0,0), F(0,b), F(a,0) 
and F(a,b ) 
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Considering first equation (B5), substituting into it the series 

expansions for equations (15) and (B27) to ('B29), eliminating e , i^, 

and p through (B37) to (B39), eliminating a', a", a’", and a"" through 
mn nnm m 

equations (Bl4), and then solving for a one obtains 


mn 2 

mnrc 


[(-l) m+n F(a,b) - (-l) m F(a,0) - (-l) n F(0,b) + F(0,0)] 


f T +-- [C ' - (-l) m c"] C 0 + | ~ [g 1 - (-l) n g"]C n 
3 mn a a n ' n 2 bb s m v ^m ] 


+ - — (— ) [B T - (-ljV'HC— ) C Q + (~) (C. - 2C )] 

a a nir n ' n v a 2 ' b ' 4 3 


+ f “ O K' - <-«“ c i + O - 2C 5 )U 

(B50) 


where 


E mn ' + (T* 2 + C l^ (WD 

Thus the unknown a have been obtained in terms of the smaller class 
mn 

of unknowns, c”, g^, g", through the compatibility equation. 

Turning now to the stiffener equilibrium equations (B 6), substituting 
the series from equations (B25), (B29) and (6), and utilizing equations 
(B^O), one obtains the relationships 


[(-l) n P x (b) - P i(0 )] + 25 s , . q. . £ 


(n=0,l,2,...,N) 


(B52) 
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2-X 


M 


no 


b [(-l) n p 2 (b) 


P 2 (0) ] + ^ s n + «n + [w'v- 0 


m=0 


2-5 


(n=0, 1,2, ...,N) 


N 


(B53) 


mo 


[(-!)“ p (a) - p ( 0 )] ♦ as s;" 


2-5 


n=0 

(m=0,l,2,...,M) 

N 


(B54) 


mo 


[(-l) m p 4 (a) - p 4 (0)] ♦ = + <C + Z t- 1 )" Jmn = 0 

(B55) 


n=0 

(m=0, 1,2, ... ,M) 


From equations (Bj6) and (Bl4), it is noted that 


J 00 = " ^(a,0) - F(0,b) + F(0,0)] (B?6) 


Jon h (B A - B n> 


j = ~ — (B ,M - B"") 
mo b nut m m 


for n ^ 0 
for m / 0 


(B57) 


(B58) 


Using these results in the n=0 and m=0 equations (B52) to (B55), one 
obtains four expansions for Jqq, of which three are redundant because 
the structure as a whole in equilibrium (see ref. 1). Selecting the 
first as non- redundant, 

M 

j = -q' + i (p (b) - P, (0) + V — (B"" - B'")) (B59) 

u 00 ^o b v 1' ' l v ' L-. mit v m m ' \ ssi 

m=l 

Finally, substituting into equations (B7) the series expansions 
from equations (4), (.13), (Bl6), (B20), (B21), (B23), and (B24), one 
obtains 


63 


s */ (A E n ) + T ' - C c ' + C,B ' = 0 
n' v 1 l y n 2 n 3 n 

s "/ (A 0 E ) + T" - C 0 c" + C,B" = 0 

n /v 22 ' n 2n 3 « 


s'"/(A,E,) + T'” - CL g* + C,B ,n = 0 

m 7 ' 3 3 hi 1 un m 

s""/(A ) E ) ) + T"" - Ck g" + CHB"" = 0 
m ' 4 4' m 1 nn 5 m 


) (b6o) 


Reduction in the Number of Simultaneous Equations 

Essentially the problem has now been reduced to the solution of 

equations (B 50 ), (B 52 ) to (B55) and (B60) for the unknowns a^^ s^, 

S n’ S m"' S m"* C n’ C n> ^ and Since equation (B 50 ) explicitly 

expresses a^ in terms of c^, c^, g^, and g^, the solution of equations 

(B52) to (B55), and (B60) for s', s", s'", s"", c', c", g\ and g" will 

' n’ n’ m ’ m ’ n* n' hn °m 

be sufficient. However, further reduction in the number of equations 
to be solved may still be attained by using equations (b6o) to express 
s', s", s^" and s^" in terms of c^, c^, g^ and g^. Equations (B52) to 
(B55), with the n=0 and m=0 equations excluded, can then be used to 
obtain a system of simultaneous equations in which the c^, c^, g^ and 
g^ are the only unknowns. This is accomplished by eliminating s^, s^, 
s'", and s"" with the aid of equations (b6o), and j by means of the 
following expression 
2 


mnir 


°mn 
(mjfo) 
( n /0 ) 


abE 


mn 


(T + — — 
mn a a 




+ = r - (§(— ) 3 [ B ’" - (- l ) 1 ^""]^ + - (^£) 3 [ B ' 
E b'a m N m 2 a b' r 
mn 


(-i)^]c 1 ) 

(b6i) 
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•which is obtained from equations (B 36 ), (Bl4), and (B50). 
form of equations (B52) to (B55) is: 




- K ' ( - 1)n K 


E 


(x 


m=l 


mn 


M 




-'=• Pl (n) + 5; <x>) - bj + § (£S) 2 £ (-D m 


E 


m=l 

N 


mn 


% 


; a 5 w - Z - 


c’ - (-l) m c" 

M » 7 VI 


n 


E 


(m= 


n=l 

N 


mn 


+ 5; “i>> = C - !<“> Z ( ‘ 1)n ~ 


c' - (-l) ra c" 


E 


n=l 


mn 


where 


C n = C n c 2 ( n V b Z c " = c n c 2 W b Z = g^C 1 (nut/a), i 


M 


“l< n > - A A + f Y. 


(mjt/a)‘ 


E 


m=l 
M 


mn 


a 2 (n) . A 2 E 2 + f Y 


(nut/a Y 


m=l 

N 


mn 


“ 3 W =V 5 + fZ ^ 

mn 


n=l 

N 


a l>> = A A + § £ 

_ mn 


M 


Pl(n) - f Z 


n=l 


(-l) m (mrt/a) c 


E 


m=l 

N 


mn 


M-)-f Z^™ 

_ t mn 


n=l 
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H 

L 


The resulting 

L=lj 2, . • • jN) 

(B6e) 

(B 65 ) 

‘1,2,... f M) 

(b64) 

- (m=l, 2 , . . . ,M) 

(B65) 

m = Eri C lW*) 
(B 66 ) 


(B 67 ) 



and R', R", R'”, R"" are the following combinations of known thermal 
n’ n' m ’ m e 

and loading quantities: 

R' = O' + r[P.(0) - (-l)V(b)] + — (B'-B") + A E ^(C B ' + T’) - V K 
n hi b L I s ' s ’ I s ann'' n n' 1 1 b v 3 n n' / i 

m=l 

M 

R" = -Q" + f[P o (0) - (-l) n P 0 (b) ] - —(B'-B") + A E — (C,B" + T" )+ Y (-1)“k 
n hi b 2 v ' 2 amt ' n n ' 2 2 b ^ 3 n n v ] 


M 


mn 


m=l 

If 


R"’ = q." + %p (o)-(-l) m P z (a)]+ r-— (B ' "-B 11 ") + A,E, — (C,B"'+T"’) - Y K 
m hi a 3 3 bmir x m m ■ 3 3 a v 3 m m ' / . mn 


n=l 


N 


R n " = _q"'< + — [ P (0)-(-l) 1D P l (a)]- r~(B”'-B"")+A,E ) — (C,B""+T"")+ Y (-l). 1 ^ 
m hi a 4 V v 4 bmn m m 44 a 3 m m v mn 


bnur' m ^m ^"4*4 a^3 m m ‘ L* 

n=l 

(B 68) 


with 

K = ) 5 C 0 [B’" - (-1)V"3 - -(^) 5 G [B 1 - (-l) m B"]) 

E abmn b> a 2 m x/ m a b I n ' ' n 


ran 


mn 


(B69) 


Equations (b 62) to (B 65 ) can be solved simultaneously for c^, c^, 
g^, and g^. With these known, equations (b 60) will furnish the values 
of s£, s^", s^", and equations (B57) to (B39) and (B6l) the values 
of the .1 . Equations (Bl6) will then give the stiffener stresses, 

equations (B19) to (B25) the plate stresses. 



APPENDIX C 


ANALYSIS FOR THE CASE OF ONE STIFFENER 
WITH PRESCRIBED DISPLACEMENT CONDITIONS 


In this appendix the case of figure 4a is considered. In this 

case the edge x = 0 of the plate is assumed to be forced into a 

prescribed shape by means of an attached rigid stiffener (shown shaded 

in fig. 4a) which also has this prescribed shape. The shape is defined 
2 2 

by (8 u/dy ) , which is assumed to be a given function of y and 

expandable in the form of a Fourier series, equation ( 8 ), with known 

Fourier coefficients (n = 1, 2, . N). 

By virtue of the new conditions along the edge x = 0, certain 

quantities which were considered to be known or given in the previous 

appendix are now unknown. These are (a) the N Fourier coefficients 

B^ which, through the first of equations (4), described the running 

tension N x ( 0 ,y) acting mutually between the stiffener and the edge of 

the plate, and (b) the tension P^(0) and P^(0) existing at the 

ends of the horizontal stiffeners. Because the stiffener along x = 0 

is now rigid, these N+2 quantities, needed in equations (B 68 ) and (B 69 ), 

are no longer known from the given external loading, but must be regarded 

as additional unknowns along with the c^, c”, g^, g^. 

What makes the present case still solvable in the face of this 

increase in the number of unknowns is the fact that there are exactly 

N+2 new conditions which must be imposed on the problem. N of these new 

2 2 

conditions state that the curvature o u/dy along the edge x =0 of the 
plate, produced by the stresses in the plate, must, when expanded in a 
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Fourier series, be consistent with the known Fourier coefficients 
(n=l,2, . . . ,N) associated with the prescribed shape of the edge. 

The remaining two of the new conditions express the fact that the 
rigid stiffener at x=0 must be in equilibrium under the action of the 
unknown P,.(o), P^(0), N^(0,y) and the known T^ and M ^ (see fig. 5)* 

It should be noted that P^(0) and P^(0) are no longer applied by an 
agent that is external to the entire structure, but are applied by 
the rigid stiffener along x=0. They are forces that now act mutually 
between the rigid vertical stiffener at x=0 and the two horizontal 
stiffeners. 

The explicit mathematical formulation of these N+2 new conditions 
and their incorporation into the analysis of the previous appendix will 
constitute the bulk of the present appendix. 


Formulation of Boundary Condition of Prescribed Curvat ure 


Dirf^rentiating the last of the strain-displacement equations, (Bl), 


with respect to y, one obtains 


>2 de 

= *y _ _z 

dy 2 dy dx 


Eliminating the strains in terms of the stresses by means of equations 
(2), and then the stresses in terms of the stress function through equations 
(B4), this becomes 


= (c 3 - c k ) - b ^ F 

dy ^ dxdy‘ 


2 C2 sx 5 


S J F _ 


Sx 


(Cl) 


2 2 

Thus the curvatures 5 u/c)y of the edge x=0 of the plate are 
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n2 • 

By x=0 


- (ty^x 


■%> - c a (4) 

BxBy x=0 Bx^ x=0 


Be 

Os-X) 

Mx ^x=0 

(C2) 


The terms on the right-hand side of this equation can be expressed in 
series form with the aid of equations (19), (B4l) and (b 44). The 
result is 


>2 

£f> - 

By x=0 


N M 


n=l 


m=0 


mn 


M 

I 

m=0 


mn 


- V'] 
n 


sin 


nity 

b 


(C3) 


Comparing equations (C3) and (8), one obtains the following N equations 
representing the condition of prescribed curvature along the edge x=0: 


M M 



m=0 m=0 


(n=l,2,...,N) 

(c4) 


The unknowns j and d^ in this system of equations can be expressed 

in terms of the basic unknowns B', c', c", g'.e". To accomplish this 

n n n m m 

it is first observed that equations (B57) and (B6l) can both be 
represented by the following single equation, in which any undefined 
quantities are to be considered zero: 


mn 


mnn: 

abE 

mn 


-{T 


mn 


2 nut r , 

Lc’ 

a a n 


(-D‘ 








+ =i- {§ (— ) 5 [B'" - |-1)V']C 

E l b v a m ' ' m 

mn 


o + (~) 3 [B'-(-l) m B"]C 1 } 

2 a v b' n ' ' n 1 


(C5) 


Furthermore, from equations (B46), (Bl4), (B50), and (B51), one obtains 


d = - [c" 
on a n 


c ' ] 
n 


(C6) 
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2 2 
/DUt \ /IlJt > 


V-nn - ^ V ' f ♦ (0, - *,)<=> Pf> ] 


(m^O) 

(n/0) 


- §(=*> <“> c* - (-D n «;] c x 


- 1 ^ (*f ) 2 te A - (- D m b »] Bl 

+ 1 *5 “ + (“) 2 < c 4 - 2C j> ] 

(C7) 

and these two equations can be represented by the following single 
equation in which, once again, undefined quantities are to be regarded 


as zero: 


E d = (— ) T 
mn mn v a mn 


2-5 2 

mo ><nrt, 


<“f £<>A - (- i) m CJHC^SS) 8 ♦ (=) 2 (<yac )] 


(Vo) 


+ f C2f) t*i - (-U B 


- f (=) <“> [B n - <- 1 > m b : 1C 1 


♦ 1 = k ’ - ( - i>n b ; ,,] « £ 5) °i + (“) < c 4 - 2c j )] 

(C8) 

Substituting from equations (C5) and (C8) into equation (c4), and 


separating B' terms, one obtains 
n 2 M 


2 2-6 


*a - -k - ^ c i 1 1; ' °8 < “ ) K ~^ E) 


m=0 


M 


— mn 


m=0 


2-5 


- Clj )(22) -c 2 (E) ][§ SStg^-DX)* -f(S)(-l)Vj 


-c s (f)[o 1 (^, 2 -o 3 (S|) 2 ]c|E(B;"-(-i)V”)- 2 ' 5mo “ 
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This expression can be solved for each unknown in terms of the 
corresponding V^, c^, c^ and all the g^ and g^, with the following 
result: 

& 


M 


M 


B A - 701 + JIT Y. “An^A'^ 15 ” ®A ] + JIT 

teg) 


'n m=l 


'n m=0 


.CD 


where 7 ' ,6 ' and v' are known quantities and are defined as follows: 

/ n ’ n mn mn 

4 l) - ^ C lZf t«V<V^) 2 - C 2 (2S) 2 ] ^S2 (CIO) 


m=0 


mn 


M 


2 2-5 


6 —K’-V'+V -^-{[<C -C,)(^) -C 0 (5S) ][T (— )+C (^) ( ■ m? )(-l) m B" j 

n n n / . E v 3 4 ' b ' 2^ a mn v a l' b' v a v n 

1 mn 


m=0 


-CgLc^ag) 2 -^^) 2 ] ss | = [b;-(-i)V«d 


m 


,1 _ r( , G _ c n/HUt _ n (EE\ t HU ^ 
mn E L ^ 3 °4^ b j °2^ a^ J ^ * 

mn 

2 


b b 


1 „ _ 2 2-5 

. -L r _ / Hit \ _ /imr \ 1 nut mo 

V mn = E~ [C 1 <~b> - C ^l) J T ~ 

mn 


(Cll) 
(G12 ) 

(C13) 


Thus the unknown have, in effect, through equation (C 9 ), been replaced 
by an equal number of known K^. If the edge x=0 is forced to remain 
straight, then the KJ are all zero. 


Formulation of Boundary Condition of Equilibrium 
The normal forces acting on the rigid stiffener (see fig. 5) must 
be self-equilibrating. Therefore 


D 

P 3 (0) + P 4 (0) + J N x (0,y) dy = T 1 


and 
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r T b 

P^(0}b + j y N x (0,y) dy = -g- + 


Substituting 


N 


N 


: (0,y) = £ B n Sin (B ? } 


n=l 


and solving for P^(0) and P^(0), one obtains 

N m 

r-— i -u T M 

pJo) = -> — B* + 

3 nn n 2 b 


n=l 


N 


T„ M. 


P k (0) = Y (-l) n — B' + -^ + 
4 Z_ i ran 2 


n=l 


(C14) 


(C15) 


Thus, in effect, the unknown P^(0) and P^(o) have been expressed in terms 


of the known T^ and M . 


Separating B' Terms in R', R" , R"', R"" 

— c 52 — n n 2 — n 2 - — m- 2 — m- 

Equations (Cl4), (C15), and (B 69 ) can be used to eliminate P^.(0), 
P^(0).> an< 3- from equations (b 68 ). If the terms are then written 


•from 

the rest, 

equations (b 68) become 

R' 

+ 

1 — 1 
v — f 

CO 

II 

7 ’B' 

n 

n 

'n n 

R" 

= s< 2 > - 

7 "B' 

n 

n 

'n n 



N 

R’" 

II 

+ 

Y B' H 

m 

m 

n mn 



n=l 


- S W - 

N 

R"" 

Y (-D n B' H 
Z_, n mn 

m 

m 



n=l 


(Cl6) 
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where and are completely known loading terms defined 

n ’ n ’ m m 


n ' n ' m 
as follows 


M 


gC 1 ) e q. + |[P (0)-(-l) n P-,(h)] - — B" + A E.. 22 T' - Y 
n ta b l v ' 1' ann n 1 1 b n mn 


m=l 

(C 17 ) 

i 2) - -<fi + I [y°> - *~sK + A 2 E 2 H + T n> 

M 


£ (- 1 )" s£> (CIS) 


m=l 


s i 3) - V -t - 1 )” 1 S P 3 (a ) 


, V a, • (B'" - B M ") + A X E, — (C,B ' " + T'") 

3 V bmn v m m ) 3 a 3 m m 

N 


-I 


(i) + h . f!i 

ab 


S'*' + 
mn a 


(C 19 ) 


n=l 


S (4) = _q"" _ ^(-i) m P ( a ) - r~ — (B 1 11 - B" M ) + A, Ei — (C B" " + T"") 
m rn a' ' 4' 7 bmn v m m 4 4 a 3 m m 


mn 


N 


T- 2M, 


+ V (-l) n S^ + — + — g (C20) 

/ v mn a ab 


n=l 


with 



s« 

mn 

1 

“ E 

j-mnn 

ab 


mn 


H , 
mn 7 

y' and 
1 n 

n 

?n 


m 


a' b' ~1 n 
(C21) 


equations: 


M 




(C22) 


m=l 
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M 


y n arur 


I (-D m | O o 1 

_ mn 


m=l 


and 


H 


1 2 /nit 


2 b 


mn E 


mn 


- (~) C - — — 
a b' 1 a nit 


(C25) 


(C24) 


Revision of Equations (B62) to (b6^) 

Substituting from equations (Cl6) into equations (b 62) to (B 65 ), one obtains 


c’ a n (n) 
n 


-c" Pi(n) 


0 M ” . / \ n — tf 

= s* 1 ) + r 'B' - f( 2 £) V 

n 'n n b' b' Z_ E 

, mn 

m=l 


(n=l,2, (C 25 ) 


o M “ 1 / ijH “tt 

■’C fS.(n) + c" a„(n) > S (2) - r "B' + S(22) y (-1)” SuLlI 5s 

n ' n 2 X ' n 'n n b' b' 7 E 

n mn 

m=l 


(n=l,2,...,N) (C26) 


N 


N 


g’ a,(m) - g" P 0 (m) = + V B’H “ -(— ) V 

m3 tii 2 V m nmn a' a' ^ 


c' - (-l) m c* 


n=l 


n=l 


mn 


(®=1,2,...,M) (C27) 


N 


N 


-g’ p (m) + g" a, (m) = - Y (" 1 ) n B ' H + "(— ) Y “ 

ti K 2 V ' & m 4 V / m v n mn a v a' L— 


c; - C-D m c 


n 


n=l 


n=l 


mn 


(ro=l,2, . . .,M) (C28) 

If equation (C9) is now used to eliminate the B^, one obtains the follow- 
ing systemof simultaneous equations in which c^, g^ and g^ are the 


only unknowns: 
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M 


- KK - sa - 1 ^ K - <--v n Q 

(n = 1, 2, 


m=l 


N) 


M 


-W + 5 nC ■ ^ - I s ™ - t' 1 )" V 

“ =1 Cn = 1, 2, . 

M N M 

E g 1 [o' (m)6 - ^ ' H — 7TT ^ ' 1 - ^ ' g" T" 

6 p 3 V ' mp Z_. mn (1) pn Z_ & p mp 


N) 


p=l 


N 


n=l 'n 
N 


V c' r"' - V* c" 
/_ n mn Z_, n 


p=l 
N 


n=l 


r"" = q(5) + V -X H 
mn m 2—> (l) mn 

n=l 'n 


(m=l,2, 3, . ..,M) 


M 


M 


N 


- \ g' o' + y g" [a, (m)5 - H y ■ v p ' 3 

& p mp Z_. p mp 4_ mn (l) pn 

p=l p=l n=l / n 

N 


N 


N 


- y c' a"' + Y" c" a 
x— n mn Z_. n 


"" = S 
mn m 


(M 


t — » u 

- 1 


H 


mn 


n=l 


n=l 


n=l 


n 


(m=l,2,5> • • ■ ,M) 


where 


M 

0 n - a l (n) -jfr £ 

' — r\ 


n m=0 
7' M 


K ‘ Pl<“> -TIT I ( - 1)B V in 


n m=0 


7 „ M 

C “ p i (n) " TI) I. 


mn 


n m=0 


(C29) 


(C50) 


(C3i) 


(C32) 


(C33) 
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nr - v»> - -rfy I (- 1 )' 


n m=0 


f. n.(l) , 7 n 8 n 

c n - S n -TlT 


5" = S (2) - ^2 
n n (1) 


yl 2 

*. = _. n . u . _ 2 ,n*v _1_ 

5 mn (l) mn b ' b ; E 

mn 

n 


^ a • .2 il±£ 


t” = a - -a. JC. 

s mn ( 1 ) mn b ' b ' E 


r mp “ P 2 (m)5 mp " E ( " 1)nH mn 'JIJ 11 , 


P» it _ 2 / inn % 1 _ 1 p. 1 , 

mn ~ a ^ ~ 7 (1) H ®n Z_. V 

7 n p=0 


(-D"|(5f) ji- 

mn 
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N 


j' = 8 B (m) - V (-l) n H -TTT^' 

mp mp K 2 V x mn (1) pn 


n=l 


'n 


o’" = f-i) n r- (— ) — - 

mn > ; v v. 


a ' a' E H mn (l) 


mn 


M 

I’i 


n p=0 


pn] 


(C37) 


M 


o”" = (-l) n [(-l) m - (— ) - H 

' 1 L ' a a E mn 

mn mn 


try Z ( ' 1)P v pn J 


T 17 

7 n p=0 


Some uncoupling of unknowns in equations (C3l) and. (C32) can be 
effected by adding and subtracting equations corresponding to the same 
value of m. By adding, one obtains 


M 


M 


N 


N 


Z g p % + Z K 9 mp + Z ’Cn - Z =n C - Pi 

P~ 1 p~ 1 1,3. • . n— 1,3*.* /p ^ q > 

(n odd) (n odd) ' 


By subtracting, one obtains 


M 


M 


Y g’e"' - Y i" G "" + Y 

L_> °p mp °p mp i — * 


N N 

c* n’ - Y 


P=1 


P=1 


where 


n T \nn Z— C n T Wi 

n=2,4. . . n=2,4. . . 

(n even) (n even) 


N 


!f _fl M 


= P, 


m 

(C39) 


0 ' = fa,(m) - P_(m)l 5 - 2 Y H h ' 

mp L 3 ' ' 2' ^ rap c— mn (1) pn 


n=l, 3**. 7 n 
(n odd) 

N 


(C40) 


9" = [ a, (m) - f3 0 (m)l S> - 2 Y ^ /■ \ P 1 

mp L 4' ' K 2 V '1 mp mn £1) pn 


n=l, 3* • • 

(n odd) 
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M 


, _ 4 /inrtv 1 2 V" , 

\n a ' E (l) H mn 2- V pn 


mn 7 


n 


p=0 


M 


<n “ <-*> m I O t ' III H *n 1 1- 1 )" *p 

v n p=0 


mn 


pn 


(C4l) 


N 


©"' = [a,(m) + £ (m)J 5 - 2 V 

mp 3 2 ' mp Z_. 


JJ ^ I 

mp Z_. mn (l) H pn 


n=2, 4 . . . 
(n even) 


N 


mp 


[ «kW * % ■ 2 1 V -rry^p 

n=2, 4 . . . ^n 

(n even) 


pn 


(C42) 


p , . s (j) + a w + 2 y _!a 

H m m m { 1 

n=l, 3 . . . ^n 
(n odd) 


T H 
) mn 


N 


- S (3) . s « + 2 V 

m m L - . 


n 


H 


n=2,4. . . 7 n 
(n even) 


(IT mn 


(C43) 


Equation (C 38 ) involves only the odd subscript and and equation 
(C39) only the even- subscript and c^. These equations may replace 
equations (C3l) and (C32). 

Reduction in the Number of Simultaneous Equations 
Equations (C29) and (C30), written for the same value of n, can be 


solved for c' and c" in terms of all the g' and g". The result is 
n n & m & m 
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M 


!' « ~ ft’ 0"" + t"0" + V (0"V - 0'V ) [g 1 - (-l) n g"j} 

n $ n lb n mi vy n s pn Mi s pn ; Le p v ; ®p J 


P=1 

M 


(C44) 


3" = i- 

n $ 


n 


f W + % + I>A'V - W c ^> '- <-u n S j) 

P" 1 (C45) 


where 


$ = 0'0"" - 0"0 ' " . 
n r n r n r n r n 

Utilizing equations (C44) and (C45) to eliminate the c^ and c^j 
in equations (C38) and (C39)> and combining like terms, one obtains the 
following simultaneous equations involving only the g^ and g^ as unknowns: 


M N M N N 

Z g’ {9' + y »if i + y g" (©" + y t ? = p* - y 

D p mp 4_ mnp °p mp Z_ mnp p m / 

P=1 n=l, 3- • . P=1 


* 


mn 


(n odd) 


n=l, 3* • • 

(n odd) 


n= 1, 3 • • 

(n odd) 

(C 46) 


M N M N N 

y < (e + y + j - y < (e"" + y * j = P " - y 

» p mp t— • mnp / ^ p mp * mnp m / . 


mn 


p=i 


n=2, 4 . . . 
(n even) 


p=l 


n=2 , 4 . . . 
(n even) 


n=2, 4 . . . 
(n even) 


where 


4 ' 4 

i = rSS(0 ,M, |* - 0"i" ) - (0«»|« - 0-|" ) 

mnp $ xr n pn r n pn r n pn r n pn 


(m = 1, 2, . . . , M) (C46) 

(C4T) 


4 ' 4” 

* = -^(t’0"" + t"0" ) - (t’0’" + t"0 1 ) 

mn 0 V5 n^n *n r n' <J> ' a n r n 3 n^n 

n n 


(C48) 


The advantage of this reduction is evident: whereas the original simultaneous 

equations system, equations (C29) to (C32) requires the solution of 2N+2M 
simultaneous equations, the reduced system, equations (c46), contains only 
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2M simultaneous equations. Thus N may he taken arbitrarily large with- 
out increasing the number of simultaneous equations that have to be 
solved. 


Procedure for Use of Equations 

The procedure for using the foregoing analysis will now be summarized: 
Equations (C46) are first solved for the and g^. With these known, 

equations (c44) and (C4-5) give directly the c^ and c^, and equation (C9) 
the B^. Equations (b6o) then give the s^, s^, s^", s^", and equations 
(B57) to (B49) and (B6l) the j . Finally, equations (Bl6) and (B19) to 
(B25) give the stiffener and plate stresses. 


Special Case: Symmetry About y - b/2 

When the structure and loading are symmetrical about the line 

y = b/2 considerable simplification of the foregoing equations is possible. 

The symmetry implies that = A^, e^(x) = e^(x), P-^(O) = P^(b), P^O) - 

P 2 (b), P 5 (0) = P 4 (0), P 5 (a) = P 4 (a), N^x) = N 4 (x), q ? (x) = -q 4 (x), T 1 

and M n = 0. It also implies that B'" - B"", O'" = - Q"", T'" = T"", and 

a^(m) = a 4 (m). In addition, as a result of the symmetry the following 

quantities all vanish for n even: K', T , B", O’, Q", T', T", and V'. 

^ n' mn' n' tt hr n’ n' n 

Consequently the following quantities vanish for n even: 


8< lj , S^ 2) (see eqs. (C17) and (Cl8)); B n (eq. (Cll)); (eqs. (C34)), 


^mn ^ eq * • 


and the following equalities hold: 


S<*> 



(eqs. (C19) and (C20)) 


0 ' = 0 " 
mp mp 


(eqs. (C40)) 
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(eqs . (C42)) 


0'" = 0”" 
mp mp 


d" = 0 
m 


(the second equation of eqs. (C43)) 


Therefore, equations (C46) become 
M N N 


Z (i p + ' g ; )Ce m P + Z 


Z 


if. 


Tlf ) = o' 

J p w p ' ' mp i — » mnp nm l—> ' mn 

p=l n=l, 3... n=l,3**» 

(n odd) (n odd) 


N 


N 


>(C49) 


V (g 1 - g")(e ,M + y if ) = o 
Z_ i p P mp £— > mnp' 


(m = 1,2, . . .,M) 


P=1 


n=2,4. . 
(n odd) 


Whence 


and 


M 


_ I tl 

g p “ S P 


N 


(p = 1, 2, . M) 


Z ®P (e mp + Z 


* 


mnp 


p=l 


n=l, 3, ■ • • 

(n odd) 


N 

> - - Z 

n=l,3, ... 

(m = 1, 2, ..., M) 


(C50) 


(C51) 


From equation (c44) and (C45), together with (C50) and the earlier 
consequences of symmetry, there results 


c 1 = 0 
n 


c" = 0 
n 


for n even, 


(C52) 


M 


and, for n odd, 

1 fC'0"" + C"0" + 2 Y (0""| ’ - 0"£" )i'} 

rtn b n^n vy n pn Vpn ;6 p 


c ' = 


n 


n 


P=1 

M 


(053) 


c" = i- 
n $ 


n 


{£'0'" + £"0' + 2 V (0 '"I ' - 0'|" ) i 1 ] 

n*n ^n’n Z_. v> n pn Vpn ; 6 p 

p=l 
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Similarly, from equation iC9) 


B' = 0 for n even, 
n 


(C5*0 


and, for n odd. 



fc' - (-D 


m 


n 



(C55) 


The procedure for the symmetrical case can now be summarized as 

follows: Solve equations (C5l) simultaneously for the g', then use 

P 

equations (C53) and j[C55) to compute the odd-subscript c^, c^ and B^. 
With these known, equations JiB6o) will furnish the values of s^, s”, 
s^", s^" (n odd), and equations (B37) and (B6l) the values of the j 
(n odd). Equations (Bl6) and (B19) to (B25), with the only odd values 
of n included, will then give the stiffener and plate stresses. 


Limiting Case of Large Stiffener Areas 
The case in which the stiffener cross-sectional areas are large 
compared with the plate cross-sectional area is of practical and 
theoretical interest. In order to study this case, let it be assumed 
that A^E^, A^^, A^E,, and wm- approach infinity while maintaining 

constant ratios with each other. Then equations (C29) to (C32) can be 
simplified through the following steps: First, divide equation (C29) 

by A 1 E 1 , (C30) by A^, (C3l) by and ( c 32) by A^E^ and rearrange 

them to obtain 


c Ei-: 

n 


A' 


n + A,E 


A" 

-] + c"[~' - 

1“1 n n A 1 E 1 A 1 E 1 


M 




mn 

,(D 


nxr r 

b 


m=l 


A'” 

+ -™L] 
A 1 E 1 


(C56) 
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I'” u 

•c' — — + c" [l + — — ] 
C n A 2 E 2 C n L1 A 2 E 2 J 


c: 


M 


E 0 “ A 0 E 0 7 ^ ^ ^mn 


2 2 2 2 


m=l 


(C57) 


M 


M 


7 s' [&_. + a^ t t " ^ - a 7 ' 7 g r, r m 

L_* p mp A E,, A-,E, ^ p m 

p=l J J J J p=l 


mp 


A 3 E 3 


N N 

Z c' p 1 " - ~i — >* c" r"" = 

n mn A E n mn 

3 3 


n=l 


n=l 


S'" 

m 

A 3 E 3 


(C58) 


A 4 E 4 


M M 

V £' a- + V g" [6 + ~S£_] 

S p mp 2^ § P m P \\ 


P=1 


P=1 


N 


N 


six E C + Z =n C = or (C59) 


4 4 


n=l 


n=l 


4 4 


where 


M 


M 


n 


p -yH «T 

, = 2 V (mn/a) n ST 

* * L * 7J) Z_ 


, mn 7 

m= 1 n m=0 


mn 


(G 56 a) 


M 


V = p.(n) - -4=-r V (-I)® V 
n (1) L- mn 


n m=0 


(C56b) 


u' 2 

mn __ 2 ,rnt^ 1 

"" 'K V v-, ' tr 


A'" - —r 7 

mn (1) / n b x b' E 
'n 


mn 


(G 56 c) 


with 


M 

y"" _ J 3 + V"" ^ ('BZl') Q 

^n anm ) E a ' b' 

L ~t mn 


m=l 


(C56d) 


M 


'n Tl) 

'n m=0 


Z , nn 
V mn b 


(C56e) 
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(C56f) 


M 


e A = 7(T) 


y (-i) m »■ si c, 

L_, mn b 3 


n m=0 
M 


M 


vIZ I <-« 

m=l ™ 7 n m=0 


mn 


N 


N 


2 y (njr/b _) _y 1 

b L E_ Z- mn 7(1) pn 


V =6 

mp mp - - r 

n=l n=l n 


N 


S .(3) = s (3) + y “n 


H 


mn 


s<*> - s< 4 > 

m m 


n=l 

E 8 

Z ( - 1} 7?) H mn 

n=l / n 


(C5Ta) 

(C 58 a) 

(C58b) 

(C59a) 


Examining the coefficients of the unknowns c7, c”, and g^ in equations 

(C56) through (C59), it is observed that some of these coefficients are 

2 

of the order of 1, while the others are of the order of l/(a E ^A-jE^). 
Retaining only terms of the order of 1 in these coefficients, one reduces 
equations (C 56 ) to (C59) to the following system: 


c ’ 
n 


[1 - E n ] - 


J 7 

A i E i 


M 


C _ 

n n 


Z [i A - ( ' 1>n 


m=l 


M- 

mn 



b ~3 

(c6o) 


c = 
n 


JL 

A 2 E 2 



g 


II 

m 


5 ( 5 ) 

m 

A 3 E 3 

3(4) 

m 

Vu 


(C6l) 

(C 62 ) 

( 063 ) 
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Using equations (C6l) through (C 63 ) to simplify equation (C 60 ), the 
latter becomes 



(c64a) 


Thus a solution >(eqs. (C6l) through (c64)) is obtained which gives 
the unknowns c^, c^, g^, g^ explicitly without the necessity of solving 
simultaneous equations. This solution can be seen to be correct to terms 
of the first degree in l/(a^E^^A^E^) . With c^, etc., known, the 
procedure for computing stresses is the same as described earlier for 
the general case. 

A solution correct to terms of the second degree in l/ (a^E^A^E^) 
can be obtained by the following procedure: First substitute from equations 


(C 62 ) to (C64) into (C57) to obtain 


c = 
n 


1 


U 

[1 + — £-] 
W 


t" 

(^s_ + 1 


A 


n 


A 2 E 2 


• r-i 


A 2 E 2 A 1 E 1 n 


M 


A 2 E 2 


I < 

m=l 


g(5) 

m 

Vj 


(- 1 ) 


n m 

A 4 E 4 


] 5" ) 
mn 


(C 65 ) 


Expanding 


U 




A 2 E 2 


in a power series of the form 


1 + c 


1 - e + 


- 2 - 
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and neglecting the terms which are powers of l/ (a^E^A^E^,) higher than 
second in equations (C 65 ) gives 


fit y" 

c" = —2- + 5- 

n A 2 E 2 


< a 2 V‘ 


(C66) 


where 


M 


Y" 

n 


'Vn + A 2 E 2 


C - a 2 e 2 I [ ; 


m=l 


g(5) 

m 

A 3 E 3 


S (4) 

(-l) n _®_] 

A ^ A 4 E 4 J 6 inn 
(C66a) 


By further application of this technique, using the first order approxi- 
mations given by equations (C6l) to (C64), the second order approximations 
for and g^ (from eqs. (C 58 ) and (C59 ) ) can be expressed by following 


equations: 

g v y ) y T M 

- , _ m + m 

^ " a j E 3 (a 3 e 3 ) 2 

S< 4) Y" n 

-it m + m 


^ A 4 E 4 ' (a 4 e 4 ) 2 


(c6t) 


(C68) 


where 


Y 1 " 

m 


M N 


-V S^+ 


II 


p=l n-1 
(p^m) 


c(3). 


1 


x -TTT^' + A,E, / 

mm m L-, P mn (1) pn 3 ; L. 


M g(4) 


r" 


3 3 4-- A 4 E 4 mp 
P=1 


N 


N 


' A 3 E 3 


7 -a 

L, A E 


r"' + a,e 


n=l 


1 1 mn 


3 3 


C" 

y b n p 1TII 

A 2 E 2 mn 


n=l 


( 067 a) 
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<c68a) 


Y"" = -V + A, E. 

m mm m 4 4 


M q( 3 ) m n 

7 !L + V V s<*> H i 

Z_ A E L-> P inn (1) ^pn 

2 2 p=l n=l ' n 

(p^m) 


P =1 

N 


N 


r 


+ a a a a Z 


mn 


n=l n=l 

Substituting from equations (c 66 ), (C67), and (c68) into equation (C56), 
one obtains 


[1 


£ ' 

— i f _!a_ 

, + _^L] A i E i 

-n A lEl J 


M 


I « 


m 


[-JL + 

A 2 E 2 


Y" 

n 


Y i " 

m 


ra=l ‘* 3*3 < A 3 E /' 


A" 

-][H' -2-] 

(a 2 e 2 ) 2 n A i E i 


s^ 4) 

- (-l) n ( m 


m 




)] 


Expanding 


p ' 

mn nn 


A"' 

r mu iiji mn -l I 

L ~(l) b c 3 AE^' 
n 


(C69) 


l/[l - H + — 2-1 
7 L n 


in a power series of the form 

1 11 


1-K+e 1 - K 1 •+ c' 1 - K 


1 (1- C + c 2 


e’ 5 + 


where 


G’ = 


1 - K 


z 

and neglecting powers of 1/ (a^E^A^E^) higher than second in equation 
(C69) gives 

A n Y' 

n (070) 


C A - a A 


1 - E n (A lEl ) 2 
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where 


M 


Y' = 
n 


- _ 1 - A'C - 

1 - — n 3 n 
n 


- — =-=— A' A n E n 
1 - - n 11 


m=l 


V 5 1 lj a A j 


3(3) 

m 



M 


?<3) 


i< 4 > 


M 


+ A^E 


y [— — c-i) n — — 3 A'" + (a , 2 y t 

1 k^ ^ 1 Ai.E,, J mn ^ 1 l' L* 

3 3 


Y> *• 

m 


m=l 


\ E k 


m=l 


-(-1) 


Y"" 
n m 


KV 


1 Pin nn n 

2 J “til b C 5 

n 


1 - _ 


C" 

A.E A' — -§- 

n 1 1 n A 2 E 2 


E ' + 
n 


C" 

A i E i k^r 2 


A" 

n 




Y" 

n 


(a 2 e 2 ) 


2 


(CTOa) 


Equations (C66), (C 67 ), (c68) and \Q10) constitute a solution correct 
to the second degree in l/(a^E ^A E^, in which the necessity of solving 
simultaneous equations is once again obviated. 


Illustrative Thermal-Stress Problem 
In order to illustrate the details involved in the application of 
the foregoing analytical results, a particular example will be considered 
which has the following characteristics: 

a) Edge x=0 kept straight; therefore the in equation ( 8 ) are all 

zero . 

b) Plate isotropic, therefore elastic constants are given by equations 

(3) • 

c) Plate and stiffeners have the same Young’s modulus E. 
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d) A x = A 2 , -A 3 - A v 

e) No force loading. 


f) Stiffener temperature constant at the value T q . 

g) Plate temperatures T(x, y) symmetrical about both centerlines 
(x = a/2, y = b/2) and varying sinusoidally in accordance with the 
following equation: 

(0 < x < a) 

T(x,y) = T + © sin sin 

(0 < y < b) 

(C 71 ) 

where 9 is a constant, representing the temperature rise of the plate 
center relative to the stiffeners, and P and Q are odd integers. 

h) Plate and stiffeners have the same coefficient of expansion a. 
These are the only specializations to be presented later, the problem 
was further specialized to the case of a square plate (b=a), with all 
stiffeners identical (A^ = A^ = A^ = A^), and subjected to a "pillow- 
shaped" temperature distribution (P = Q = 1) . 


Reduction of general equations to special case . - From the given 
temperature distribution one obtains the following plate and stiffener 
thermal strains: 


e = e = a[T + 9 sin sin 
x y o a b 


e. = OT 
l o 


i = 1, 2, 3, 4 


(072) 

(975) 


Therefore 


S 2 e S 2 e 22 

x , y _ ^ r / Pft \ , / Qn N n . Pjix . Q.ny 

— 2 + — s = ~ ae d-;) - (V 1 V 

ay ox 


(s-^) „ - ae 2 s sin 

ax x =0 a b 


( 97*0 
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and 


e 1 (y) - e y (0,y) = aT Q - Q0? o = 0 
e 2 (y) - e (a,y) = 0 

(C75) 


e 5 (x) - e x (x,0) = 0 
e^(x) - e x (x,b) = 0 

Equations (C75) reflect the absence of any temperature discontinuity 
between stiffeners and plate in this example. Substituting from equations 
(CjU) and (C75) into the right-hand sides of equations (l4), (l6), and 
<23), one obtains 


T' = T" = T"' = T"" = 0 (C7 6) 

n n m m ' 


T 

mn 


” 6 mP 6 nQ 



(P 2 + q, 2 b 2 ) 


(C77) 


and 


Pit 

V' = a© — 6 
n a nQ 


where 

B = a/b 


(C78) 

(C79) 


Due to the absence of prescribed forces, the following quantities are 
all zero: 

P 1 (0),P 1 (b),P 2 (0),P 2 (b),P 3 (a),P J+ .(a) (fig. 4a) 

T 1' M 1 ( fig ’ ^ 

B n ,B m" ,B m" (see eqs ’ (4)) 

(see eqS - (6)) 

.(it should be noted that P (0), P (0) and B’ do not necessarily vanish.) 

3 4 n 

It will be convenient to introduce additional dimensionless parameters 
and C, defined as follows: 
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(c8o) 


= 4ah/(jr 2 A^) = 

^2 = 4bh/(A 3 ) = 4bh/(jT 2 A 4 ) 

C = V A i 

and to note, from equations (B51) and (3), that 
E = (rtV^EhKm 2 + n 2 B 2 ) 


(C8l) 


Because in this example the structure and loading are symmetrical 

about y = b/2, the simplified system of equations, namely equations 

(C 50 ) through (C55), will be used for the determination of the g', 

P 

c^, c^ and B^. The quantities needed in order to use these equations 
will now be evaluated. 

Substituting from equations ( 3 ) , (C77), (C 78 ), and (C8l) into 
equations (CIO) to (C13)> one obtains 


M 


7^ - -(n 2 B 2 /aEh) ^ 


i=0 


[i 2 + (2+v)n 2 B 2 ](2-S. 0 ) 

r • 2 , 2^2 -,2 

Li + n B J 


(CIO') 


8 = 0)0 
n 


Pn 

a 


[P 2 + (2+v)Q 2 B 2 3 
[P 2 + q 2 b 2 ] 


*3 


8 


nQ 


P 


mn 


2nB 2 [m 2 + (2+v)n 2 B 2 ] 

r 2 . 2 2 -,2 

jt Lm + n B J 


v ' 
mn 


nB[n 2 B 2 - vm 2 ](2 - 8 ) 

mo 

rr[m 2 + n 2 B 2 ] 2 


(C11-) 


(C12 ' ) 


(C13’) 


Substituting into equations (C17) to (C24) from equations (C77) and 
(C8l), one obtains 
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s (l) . 5 „ 0SEh - 

n nQ r P 2 2-, 

* [P + Q B ] 




= 6 aeEh — ^ 
m "* (P 2 + Q, 2 B 2 ) 

S< 4 > = 6 p 09Eh 
" "* (P 2 + Q. 2 B 2 ) 


3 3 M 

rV y 

* Ei S 


7' = 1 + v ^nB + 2n^B' 

n - njrB ^ ^ [,2 + 


y" = _i_ + 2 **^ y 7 (-1) 1 

n mtB n Z- r . 2 2„2 -,2 

i=1 U + n B J 


M 


2n 5 B 5 


Hmn = *[m 2 + n 2 B 2 ] 2 " n * B 


(C1T’) 
(Cl8') 
(C19' 1 ) 
(C2 0 ' ) 

( C22 ' ) 

(C23 1 ) 
(C24 ' ) 


Substituting from equations (B67), (C10 ! ), (C13 1 ), (C22') and (023') 
into equations (C33)^ one obtains 


0 ' = A,E 0 T 
r n 1 r n 

0" = A,E 0" 
r n 1 >n 

0 ' " — A E 0'" 
~m J 


(C33’) 


where 0 T 
r n 


0"" = A E 0"" 

through 0^" are numerical constants 


defined as 


follows: 
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(C35'c) 



M 


-JI 


(?"» _ T + 

2 tl [±2 + n2B2]2 


M 


V 1 + ^ I c .2 1 

Li + n B J 


i=l 


M 

4n 2 B 2 V [i2+C2+V)n2B2j (2 ' V 5 


(-l) 1 [n 2 B 2 -vi 2 ](2-6. 0 ) 
[i 2 + n 2 B 2 ] 2 


i=0 


r .2 . 2 2-, 2 

Li + n B J 


(C33’d) 

Substituting from equations (CIO') (Cll 1 ) (C17 ' ), (Cl8 ' ) , (C22'), and 
(C 23 ') into equations (C3^), one obtains 


£' = Q0Eh 6 . 
m 3 n nQ 


= oGEh £" 5 « 
n *n nQ 


(C3^ ' ) 


where ^ and £" are known quantities given as follows 


M 

P{1 + + 2Q y ^ -i) 


T r = PQB 
n [p 2 + q 2 b 2 ] 


r~^[ l +Q B j [P c + Q B ] 


il [i 2 + (2+v)Q 2 B 2 ](2 - 5 ) 

Q;B d ) — 


i=0 


[i 2 + q 2 b 2 ] 2 


(C3W) 


M 


PCI - pqV X - 1} 

^ [i 2 +Q 2 B 2 J 2 [P 2 + Q 2 B 2 ] 


PQP 


^n [P 2 + Q 2 B 2 J 


i=l 


5, [i 2 + (2+v) Q 2 B 2 j(2 - 5 ) 

Q B l - 


i=0 


[i 2 + q 2 b 2 ] 2 


(c3^'b) 

Substituting from equations (C8l), (C12’), (C22 1 ) and (C23 1 ) into equations 
(C35), one obtains 
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(C35‘) 


5 ’ 

= A..E 

6 ' 

pn 

1 

pn 


n 

> 

bd 

F" 

pn 

1 

pn 

F" 

are known 

pn 




6 pn 


M 

[p 2 +(2+v)n 2 B 2 ]f^. + 4vn 2 B 2 +2k nVV — 1 - ~ - 

iTi [i + nB] \ n b5 


2 2 b f 2 + 2 _2,2 f C±2 + (2-v)n 2 B 2 ](2 - 8 ) 

2nBLp+nB] > — 


i=0 


r .2 . 2 P-,2 

Lx + n B J 


2[p 2 + n 2 B 2 ] 2 


(C35'a) 


M 

2, ,„...v_2„2 n -> V 


[ p + ( 2+ v ) n B ] \ fl + 2n B 


(-1)1 


6 " 

pn 


i=l 


r .2 2 2 ,2 

Li + n B J 


0 2^ r 2^ 2 2 .2 V 

2n B [p +n B J > 

i=0 


M [i 2 +(2+v)n 2 B 2 ](2-& i0 ) 


r .2 2 2 ,2 

[x + n B J 


(-1) P n 2 B 5 
2[p 2 + n 2 B 2 ] 2 


(C35'b) 


Substituting from (c8l) into the third and last of equations ( B6j ), 
one obtains 

x N 

a (m) = A E + A E -g B ^ “ 


and 


^ 5 5 2 ^[m 2 + n 2 B 2 J 2 

n=l 


e (m) - A E \ B 4 £ — 


n-1 


r C. , ^_2-|2 

Lm + n B J 


therefore 


N 


0=3 ( m ) “ P 2 (m) = A,E{1 + ^B^ ^ 


n=l,3, 


r 2 . 2 2 1 2 ' 
[m + n B J 


(C82) 


Substituting from equations (c82), (CIO 1 ), (C12 1 ) (C24') into the first 

of equations (c4o), one obtains 


95 


0' = A E [1 + X B 

mp j? 2 


N 

E - 

n=l,3-.. 

(n odd) 


n 


[in + n B ] 


2„2,2 1 6 mp 


H 


£ t— n a - aa - -b> ^ 

*.1,5... [ ” +nB 1 nB [p 2 tn 2 B g. jg V [1 +(2+v)n B ](2- 8 . 0 ) 


i=0 


r .2 . 2 _ 2-,2 
[l + n B J 


(C40* ) 

Substituting from equations (C8l), (CIO'), (C13')> (C24') into equations 

(C^l)^ one obtains 

tT =A.eV 
*mn 1 'mn 

* A 1 E 

where 


(Chi ' ) 


2 x [ bM _-A-l 

V , 1 [n a +n a B a 1 2 ^ 


M r 2 2 .2-1 /_ „ > 

^ [n B -vi ](2-& i0 ) 


'mn r 2 2 2 -.2 ' M r.2^,,^ x 2„2 n x Z_. r .2 ^ 2^2-12 

[m +n B J Lx t (2+v )n B ] (2-6. ) [i + n B } 

nB X 


•E 

i=0 


r . 2 , 2^2-12 

Li + n B J 


(C4l'a) 


(-D m \m e 


K [ 

1 [m 2 + n 2 B 2 ] 2 


nB J 


M 

I 


C-l) 1 [n 2 B 2 -vi 2 ](2-6 i0 ) 


— if _ _ 

™ [m 2 +n 2 B 2 ] 2 y, [i 2 +(2+v)n 2 B 2 ] (2-6^) ^ [i 2 + n 2 B 2 ] 2 

“ [i 2 +n 2 B 2 ] 2 (C4l'b) 


Substituting from equations (CIO'), (Cll'), (C19')> (C20'), (C24 1 ) into 
the first of equations (C4j5), one obtains 
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) r Jp" + (2 + v)Q^ Bgj_ i r Q V i 

[P 2 + Q. 2 B 2 ] UW ‘ ^ 

p' = O 0 Eh {& — i^ p - “ M ”2 p"^ } 

^ [P 2 +Q 2 B 2 ] 00 £ [i 2 +(2+v)Q 2 B 2 ](2-6_) 


n3t,3 


^ 2 Z 


iO' 


i=0 


[i 2 + Q 2 B 2 ] 2 


(C43*) 

Substituting from equations (C33’Z (C35'Z and. (C4l‘) into equation 
(C47) one obtains 


* 


A E[ q ’ (C0 ,m r - W'i" ) - q" ((#'" | t" )J 

1 'mn' r n pn r n pn' ■mn' r n pn r n pn 


-HI 


£n_ 


HL 


mnp 


C($' $"') 

r n r n r n r n 


(C47’) 

Substituting from equations (C33'Z (CH'), (C4l') into equation (C48), 
one obtains 




mn 


C(9Eh [q> (CC + £" <J") - q" (Ct, ' (?'" + £" ®')]B _ 
‘mn n r n n r n 'mn' n r n n r n nQ, 

C($' 3"" - $'") 

' r n r n r n r n 7 


(c48 ' ) 

Substituting into equation (C5l) from equations (C40'), (c43')> (c 47‘) 
and (c48'), one finally obtains the following system of simultaneous 
equations for the unknown g^ s 


W 


tZ 

n=l,3. . . 


q* ) - q" (C3"'|' - ^‘1" ) 

*mn r n mn r n mn 'mn' r n mn ^n mn 


N 


C($' 3"" - 3" S'") 

v^n *n ‘ 


- c C 1 + V Z 


n 


a=l,3*- (m 2 +n 2 B 2 ) 2 


N 


4 

■2\gB X! 


Z ^ 


, e.. 2_2 v 2 

n=l,5..! m +n B > 


-b! 

n B 


[m 2 + (2+v)n 2 B 2 j 


, 2, 2,2,2 V [i 2 *(^v)n 2 B 2 ](2-6 ) 

[m +n B ] 5 o o o Mm) 


i=0 


r .2 . 2 2 1 2 

Li + n B J 
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L 




^JtQB 

4[p 2 +q, 2 b 2 ] 



[ [p"+(2+v)Q 2 B 2 ] _ lf g^B 

A L o no J Loo 


3.5 


(P 2 + ft 2 B 2 ) 


(m 2 +Q 2 B 2 ) 2 


Ji [i 2 +(2+v)Q 2 B 2 ](2 - B ) 

Q B ) - 


i=0 


[i 2 + q 2 b 2 ] 2 



v [ ^ c W + : ^q. ( c W + W 


100 ( ¥q" " K n) 


M N 


E<I 

p=l n=l, 3- • • 




_£5L 


m. 


J£L 


C - 3"3"') 
vr n r n r n”n 


N 


■2\-)CB 


T 

n=l,3- • 


n 


r 2 2 2 ->2 2 U 

[m +n B J n B 


1 n [p 2 +(2+v)n 2 B 2 ] 
6 I, J x* 


M [i 2 +(2+v)n 2 B 2 ](2-6. 0 ) 


[p 2 + n 2 B 2 ] 2 £ 


i=0 


r .2 . 2 2 - i 2 

U + n B J 


(1 - B ) G(p) E 

v mp' m 


(G51') 


(m = 1 2, . . M) 


where 


g (p) = ^ 


O0Eh 


Procedure for numerical solution . - It will be observed that equations 
(C51') are in a suitable form to be solved by the Gauss-Seidel iterative 
procedure (ref. 6), since the equation for any particular value of m 
has G(m) as its only unknown on the left-hand side and all the G's except 
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G(m) as unknowns on the right-hand side. Solution by the Gauss-Seidel 
procedure involves the initial assumption that all of the G(p) on the 
right-hand side of the m=l equation are equal to zero. This permits the 
m=l equation to be solved for an approximate value of G(l). Substituting 
this approximate value, together with G(3) = G(4) = ... = 0 into the 
m =2 equation permits that equation to be solved for an approximate value 
of G(2). Continuing in such a fashion it is possible to obtain a set 
of approximate values for G(l) through G(M) . This set is called the 
first- iteration solution to the system of equations. A second- iteration 
solution is obtained in the same manner as the first except that the 
initial values of G(2), g( 3) , etc. are those given by the first-iteration 
solution. Third and higher iterations can be performed in a similar way. 
As one generates more sets of solutions to the system of equations 
there should appear a general observable trend whereby each individual 
G(p) tends to approach a certain value. The calculation is stopped 
when changes in all the G(p) values appear to be negligible (less than 
.000001 in the present calculations) from one iteration to the next. 

With the G(p) known, equations (C33') (C3 4'), (C35 ' ) , and (C53) 

yield 


M 


(*/4)& „[>., £ + X£"$"] + 2y' [C$'"V - ] ^77 G(p) 

v ' ' nQ L l b n^n ' 2 b n>n J Z_. *n b pn y n s pn BC 


C ' (n) = 




n(M"" - W'") 
^'rrn r n r n ' 


M 


C"(n) = 


(*/4)6 ]+2y [C^’V - ] I 77 G(p) 

x ' nQ l s n r n ‘ 2 b n r n L-, r n pn >n b pn BC ' 

Eni 


n(^’^ m ' - «'") 
v,r n r n r n r n ' 


(C53’) 
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where 


c' 

n 


C’(n) 


Q©Eh 


C" (n) = 


n 


aGEh 


With and known, equations (CIO 1 ) to (C13 ’ ) , and (G55) 

yield: 


B ' (n) = 


p p p 

. r f[P 2 + (2 + v) 0 2 B 2 ] ,1 ,7 j g(m) 

W” i ri , A.2, ' 1 f A; [rftnW 


[P‘ + 8"B g ] 


m=l 


2 Ci 2+ (2+v)Q 2 B 2 ](2-6 ) 

Q*3 2_ ~ 


M [i 2 +(2+v)n 2 B 2 ](2-6 i0 ) 


i=0 


[i 2 + Q 2 B 2 ] 2 


1=0 ti 8 ♦ n 


2 B 2 ) 2 


M 


^ [n g B 2 - vm 2 ](g-6 mo ) 

111=0 r 2 , 2 „ 2 1 2 
Lm+n B J 


[C'(n) - (-l) m C" (n) ] 


M [i 2 +(2+v)n 2 B 2 ](2-6. 0 ) 


(C55 ' ) 


I 

i=0 


r . 2 2 2 -,2 

[i +n B J 


where 


B ' (n) 


B' 

n 

aGEh 


Note that 

G( m ) = g'/(o©Eh) 

m 

C ' (n) = c^/ (aGEh) 

C"(n) = C^/ (aGEh) 

With as known, from equations (Cl4) and (C15).> one obtains the 
following tensions at the left ends of the x-wise stiffeners where they 
join the rigid vertical stiffener: 


100 


I 



N 

P 3 (o) = P 4 (o) = -o@a 3 e \ 2 Y, ( J C15 ' ) 

n=l,3-.. 


With the coefficients c^, c^ g^, g^ and B^ as known quantities, equations 
(B60) will furnish the values of s^, s^, s^", s^ u , and equations (B5T) 
to (B59) and (B 6 l) the values of the j . Equations (Bl 6 ) will then 
give the stiffener force, equations (B19) to (B25) the plate stresses. 

One obtains the following equations for these quantities. 

N 

P 1 (y) = 0A x Ea Y Cc'(n) " vB '( n )l sin (^) (0 < y < b) 

n=l,3... (C 83 ) 

N 

P 2 (y) = eA L m Y c "(n) sin (o < y < b) 

n=1 > 5 ’*’ (C84) 


M 

P 3 (x) = P 4 U) = ©A^Eo: Y G ( m ) sin C 2 ^) 

m=l 


(0 < x < a) 

(C 85 ) 


M N 

N = 09Eh Y Y G(m,n) sin (252£) s i n (22^) (0 < x < a) 

m=l n=l, 3 . • • (0 < y < b) 

(C86) 

M 

(N x } y-0 = (N x } y=b = ^ Z G(m) Sin ( ° < X < a) 

m=1 (C 87 ) 


M N 

N y = O0Eh Y Y sin sin 

m=l n=l, 3- . ♦ 


(0 < x < a) 
(0 < y < b) 

(C88) 
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(0 < y < b) 


N 

( N y)x= 0 = Q ® Eh lL sin ^^b^ 

n = l , 3 ' • • (089) 

N 

(N y ) x=a = o©Eh Y, G "( n ) sin C 2 ^) (0 < y < b) (C90) 

n=l,3- • • 


M N N 

= -09Eh{£ £ J(m,n) cos (^fO cos (^) + £ cos <?&)) 

m=l n=l, 3* • • n=l, 3- ■ • 

(0 < x < a.) 

(0 < y < b) (C91) 

where G(m,n), C(m,n) and j(m,n) are known quantities and are defined 
as follows 

Q(m,n) = o ' l- a g (m[C‘(n) - (-l) m c"(n)] + 2nB 2 G(m) ) 

it [hi + nB' J 


2m[m a + 2iA g ]Bi'(ni O^ 2 _ 5 

+ A 2 ] 2 ' ( p 2 -^ 2 ) " p 


C(m,n) 


2mB * (n) , m 2 [m 2 +2n 2 B 2 ] P 2 - 5 5 

n^it 1 [m 2 -^ 2 ] 2 ' " (P 2 + Q 2 B 2 ) ^ nQ 


2m 


ittm^nS 2 ] 2 


{m[G ' (n) - (-1) C"(n)j + 2nJT G(m) } 


j(m,n) = 


2nB 


r 2 2 ^ 2,2 

itLm +n B J 


fn^B'Cn) - m 2 [C ' (n) - (-l) m C ,, (n)] - 2mnB 2 G(m) ) 


PQB 


& _ 8 


(p 2 + q 2 b 2 ) ^ nQ 


(C83) 


102 



And from the first of equations (4), one obtains the following 
running tension between the rigid stiffener and the plate edge at x=0: 


Iff N 

U^O,y) = I B^ sin = 0©Eh ^ B'(n) sin (^p) 

n=l,3» • • n=l,3... 

(0 < y < b) (C84) 

Numerical results for \ 4 0 and X^ 4 - 0. - The numerical procedure 

^ — 

and equations described above were applied to the special case of a_ 
square plate (B = l), with all stiffener areas equal (A^ = A^ = A^ = A^), 
a pillow- shaped temperature distribution (p = Q = l)> and Poisson's 
ratio v equals to 0.3* The assumption that B = 1 and all stiffener 
areas are equal implies that \ = X^ (see eqs. (C7l)), and the common 
symbol X mil therefore be used for both X and X^. 

The results obtained for the stiffener tensions and plate stress 
are shorn in dimensionless form in figure 6 for X = 2.0, and figure 7 
for X =1.0. In general, stresses were computed at x/a and y/a interval 
of 0.1. 


Limiting case of large stiffener areas . - For the case in which the 
stiffener cross-sectional areas are very large compared to the plate 
cross-sectional area, equations (C60) to (c63) may be employed as approxi 
mations which become more and more accurate as the ratios of stiffener 
to plate cross-sectional areas approach infinity. The quantities 
E n and ZP needed in these equations are defined by equations (C 56 e) 
and (C 56 f ) which can be reduced to the following expressions if use is 
made of equations (3), (CIO'), and (C13 1 ): 
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-V 


n 


y - ™ 2 )(2 ' »_J 

m=0 


mo J 


(m 2 + n^ 2 ) 2 


M [m 2 +(2+v)n 2 B 2 ](2 - 8 ) 

' ' mo' 


I 

m=0 


(m 2 + n 2 !* 2 ) 2 


n v 


y ( 1)B (nV - vm 2 )(2-6„„) 
m=0 


(m 2 + n^ 2 ) 2 


mo' 


[m 2 + (2+v)n 2 B 2 ](2 - S^J 

m=0 


(m 2 + n ^ 2 ) 2 


mo 


(C56’e) 


(C56 1 f) 


Substituting from equations (CIO'), (Cll*), (C12'), (C19'), (C20'), 
.(C24 1 ) and (CjV) into equations (C6l) and (C62), one obtains 


c" = ctQ ~ c" 8 
n A-^ n nQ 

(c6i' ) 


5 

(C62 ' ) 


where 


+ 

n (p 2 +q 2 b 2 ) 


M 

Ci + f -r ^Va )f ?[p2 a ( 2 1 ) 2 ^ 21 - p ) 

, [m 2 + Q 2 B 2 ] 2 TP 2 + 0 2 B 2 1 


m =l 
M 

m=0 


[P + Q B ] 


QB 


Im 2 +(2+v)Q 2 B 2 ] (2-8 ) 


mo 


[m 2 + Q 2 B 2 ] 2 


(C6l'a) 


PQB 


2P(- 


[P 2 +(2+v)Q 2 B 2 ] _ 


[P^ + Q 2 B 2 ] 


l)f- 


QB 


(p^qS 2 ) ^ 


[m 2 +Q 2 B 2 ] 2 


} 


A [i 2 + (2+v) Q^ 2 ] (2-5. 0 ) 

2- r.2 2 2 2 

1=0 U QB J (C62'a) 
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With c”,g^ and known, (eq. (c6o)) can be reduced to the following 
expressions if use is made of equations (3), (CIO'), (Cl 2'), (C34 1 ), 
(C6l>), (C62 1 ) , (C56'e), CC56’f) 


v, C * 

_ h n 

C — Q0 -7 — T — 

n A, 1 - r. 

1 n 



c ' 


n 



where 


c 


n 


vPQ n 2 { 


M 

&Z 


m=0 


[P 2 + (2+v)Q 2 B 2 ] , 

\P 2 + o 2 b 2 ] 

[m 2 +(2+v)Q 2 B 2 j(2-& mo ) 


(C60’ ) 


(C 60 'a) 


M 


2„2 ! M ,2, 


g m [m c +(2+v)n^B cl ] ^ ^ ^ [m^+(2+v)n 2 B 2 ] (2-6 mo ) 


4vB 


c* = S'6 - c" H'B _ + „ 

n n nQ, n n nQ G 


r 2 2 2,2 r _2 . __ 2„2 -2 

[m +n B J 


m=l 


m=0 


[m + nB c j‘ 


(C60’b) 


with 


M 


{1+2qV^ 


1 1 f [ p 2 +(2+v)Q g B 2 J 

2, J>2 n 2 n r^.^2! 


1} P 


PQB 


-[m^+QfB^P [ P + Q B j 

m=l 


n (p 2 +q 2 b 2 ) 


qV 


M 

Z 

m=0 


[m 2 +(2+v)Q 2 B 2 ](2 - 5 mo ) 
[m 2 + Q 2 B 2 ] 2 


(C60'c) 


Substituting the above expressions for c^, c^, g^ and intp equation (C55), 
one obtains 


B' = Ct0Eh B' 6 * + O0Eh VB* 
n n nQ In 


(C55" ) 


where 
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L 


;• T 

n 4_i 


M [Q 2 B 2 - vm 2 ](2-5 mo ) 


B' = 
n 


m=0 


[n^+Q^ 2 ] 2 


Jt [m 2 +(2+v)Q 2 B 2 ](2-B_) 
(l-H n )Q* £ " 


mo 


m=0 


[m 2 + Q^B 2 ] 2 


P*f tp2 ^ V) ff ] - 1} 

[p 2 + q 2 b 2 ] 

PP A [m 2 +(2+v)Q 2 B 2 ](2-& ) 

Q B ) 


m=0 


r 2 , _2_2-i2 

Lm + Q B ] 
(C55"a) 


and 




M [m 2 +(2+v)n 2 B 2 ](2-8 mo ) 


nB 


I 

m=0 


r 2 2 t,2 -i2 

Lm - n B j 


r ILZL 
1 C 


M 

I 

m=l 


[m 2 +(2+v)n 2 B 2_1 

r 2 . 2^2i 2 

Lm + n B J 


c* 

n n V" 

4 2- 


M r 2 2 2 -i / _ - % 

Ln B -vm Jr2-8 ) 

x mo ' 


r 2. 2„2 n 2 

n ~~ n Lm +n B j 
m=0 


=„ M / ,.m 2„2 2-,/^ c >. 

nc v—, (-1) [n B - vm J (2-5 ) 

R n V mo y , 

" -- ~ j_j_ 7 o r\ O i 


nQ 


m=0 


r 2 . 2 2 x2 

[m + n B j 


(C54"b) 


With the known, equations (Cl4) and (C15) yield the following tensions 
at the left ends of the x-wise stiffeners where they join the rigid 
vertical stiffener, when terms of degree higher than 1 in and are 


neglected: 


jrB ' 


P 3 (0) = P 4 (0) . -aeA 3 EX 2 6 nQ 5^2 


(015”) 


With the coefficients c ' , c^, g^ and g" known, the stiffener 
stresses and the plate stresses are given as follows: 

* 

P x (y) = Q0A 1 E\ 1 Y [ 4^B ±^t ' vB n ] sin (0 < y < b) 

(C85) 


n= 1, 3 • • • 
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P 2 (y) 


O0AiEX-i 


I 

n=l,3. 


JtC 


5 nQ sin < 


P 3 OO 


= P,,(x) 


09A,EX 0 

3 2 


M 

I 

m=l 


Bjt ^n . ,mirxv 

T5T sin ( ~ } 


N = O0Eh 
x 


T. Z S n<i km sln ^ sin 
m=l n=l, 3- ■ • 


M B = 

(N ) = (N ) . = a0Eh\ o V sin (25£) 

x y=0 x y=b 2 Z_ Um 'a ' 


m=l 


M N 


N = OGEh 
y 


7 7 5 n e sin (252E) sin (^rp^) 

Z_ Z— nQ mn x a ' v b ' 


nQ 

m=l n= 3 • • • 


N 


(N ) _ = -a©Eh > 

y x=0 

n=l, 3 ■ 


c ' 


nQ 1 - - nir 


? sin (») 


M N 


(N ) 
y x=a 


-CtGEh 


(7 7 cos cos 


/mitx 

'nQ, °mn CUB ^ a 
m= 1 n= 1, 3 . . . 


N 

-I 

n= 1, 3 • 


B 1 


<55 8 n4l cos > 1 


where 


0 < y < b) 

(C 86 ) 

(O < x < a) 
(C87) 

(0 < x < a) 
(0 < y < b) 

(C 88 ) 

(0 < x < a) 

(C 89 ) 

(0 < x < a) 
(0 < y < b) 

(C90) 

(0 < y < b) 

(C91) 

(Sf) 

(0 < x < a) 
(0 < y < b) 

(C93) 
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1 1 ii ii ii a i !■■■■■ an in 


2mQp 2 c^ 


2 2 

- = _ Qja 6 

^nn (p 2 +Q 2 ) - \r_2.„2„2 n 2 


2m [m 2 + 2Q 2 B 2 ]B^ 


ir c ( l-^ n ) [m^+Q^B^ f jt[m 2 + Q 2 B 2 ] 2 


c - b i f^L- jL_ f m 2 tm 2 +2Q 2 B 2 ) , _ & ... 

11111 n Q 2 n B 2 [m 2 +Q 2 B 2 j 2 mP (P 2 + Q 2 B 2 ) 


2m^ c' 


n 


q/(l - H n )[m 2 + Q 2 B 2 ] 2 


(C94) 


PQB 


2m 2 Bc ' 


8 „ + 


2Q 5 B 5 B’ 
n 


11111 (P 2 +Q 2 B 2 ) mP jr 2 Im 2 +Q 2 B 2 j 2 (l-Z n ) n[m 2 + Q 2 B 2 ] 2 


With the B^ known, the first of equations (4) yields the following 
running tensions between the rigid stiffener and the plate edge at 
x = 0: 

N 

N x (0,y) = aGEh Y, 5 nQ sin ,(«) (0<y<b) 

n=l,3. • • 


Numerical results for limiting case of large stiffener areas . - The 
numerical results for = \^ = \ -> 0 for square plates (B=l) with all 
stiffeners identical, subjected to a pillow-shaped temperature distribution 
(p = Q = l) are presented in dimensionless form in figures 8 and 9* 

The former represents for v = 0.3, the latter for v = 0. 
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APPENDIX D 


ANALYSIS FOP THE CASE OF TWO OPPOSITE STIFFENERS 
WITH PRESCRIBED DISPLACEMENT CONDITIONS 


This appendix considers the case of figure 4b, in which the stiffeners 
at x = 0 and x = a are bent to prescribed shapes defined by known values 
of and in equations (8) and (9)* Correspondingly, the Fourier 
coefficients and B^, which describe the running tension between 
the stiffeners and the plate, are now unknowns. In addition the loading 
resultants T^, M^, and constitute four new knowns, supplanting 

P^(0), P^(a) and P^{a), which are now unknowns. Further intro- 
ductory remarks can be made for this case which are obvious generalizations 
of those in appendix C. 


Formulation of Boundary Condition of Prescribed Curvature 

2 2 

The boundary curvatures 8 u/dy of the edges x = 0 and x = a of 
the plate are (see appendix C) 


(% = (=3 - - C 2 (% - <J*> 

ay x=0 5 4 Sxcly x=0 ov? x=0 0X x=0 

(Jf) - (o, - - c 8 <p> - 

ay x=a oxoy x=a ox x=a x=a 


(Dl) 


The terms on the right-hand side of these equations can be expressed in 
series form with the aid of equations (l9)^(20), B4l) , and (b 44) „ The 
result is 
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N J4 M 

(^-£) = Y [(C.-C-) ~ ) j - C V d - V'] sin 

V 2' „ ' k y b 4— rmn 2 4— mn n 

ciy x=0 


nrfy 


n-1 

N 


m=0 

M 


m=0 


M 


* I «V<V “ E - ° 2 I (- 1 )“ d mn - W sin 

oy x=a 




n=l 


m=0 


m=0 


(D2) 


Comparing equations (D2) with (8) and (9), one obtains the following 2N 
equations representing the conditions of prescribed curvatures along the 
edges x=0 and x=a: 


M 


M 


K' 

n 


K" 

n 


*= < C 4 - C 3 > “ E J »m - C 2 I d , 


m=0 

M 


m=0 


- V' 

mn n 


M 


«\ - V T E ( - 1)m J mn - C 2 E (-1)"^ - V" 

m=0 m=0 

(n = 1, 2, M) 


(D3) 


The unknown and d in these equations can be expressed in terms of 
u im mn 

the basic unknowns c^, c^, g^, g” (see appendix C). 

Substituting from equations (C5) and (c8) into equations (D3), and 
separating B^ and terms, one obtains 

- -v b a(t> 2 C 1 E - o 2 (2|) 2 ](^) 

I—* mn 


m=0 


M 


2 2-6 


^(22) -C 4 )(ffi) - C a (=S) ] (-;=) 

_ ~ mn 


m=0 


M 


+ E E^~* « c 5 - C 4 )(S t ) - C 2 (! V ) ] 

^ mn 


m=0 
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-(~)C 0 Cc,(^) 2 - CJ — ) 2 ][| — (B"’-(-l) n B"")- ^5° M( c ._(_ 1 ) m c ")] 
b 2 1 b 3 a b a m v m a ^ ' r,/ 


b v n 


n 


+ !( s f)(^)0 1 I(0 5 -c,)(SS ) 2 - c 2 (SS) 2 ][g; . (-1)“ £]) 


(PM 


M 


K" = -V"-B'(^f C, Y (-l) m =£- [(C_-C, )(25 ) 2 - C 0 (22) 2 ](^-^) 
n n n x b' 1 £_. ' ' E v 3 b' 2 V a v a ' 


m=0 


mn 


M 


2 2-5 


c ilf - « E t“> K-^ 2 ) 

^ mn 


m=0 


M 


Z ^Y~ 't™ (=)[(o 3 -c 4 )(S|) - C 2 (=) ] 


m=0 




2-8 


mo nn> 


^(c--(-l) m c")] 
a b v n v n 


+ |(”) >0,1 (c 3 -c 4 ) (Sf ) _ -c a (=)' ] [ ei- «-D a <] ) 


These equations can be solved for each unknown B' and B" in terms of the 

n n 

corresponding c^, c^, and all the g^, g^. Rearranging equations 

(D4), one obtains 

M M 

B’7^ _ + y q’ [i:-(-l) n i"]+y V [c’ - (-l) m c"] 

n n n n n mn ti ' ' Tn L mn n x ' n 


m=l 


ro=0 


M 


M 


b- 7 C 2 )-b" 7 (i) = S (2) +V(-i)V LiM-DV] + y (-i) m v [c* - (-l) m C"] 

n n n n n / 1 ' mn B m' , & m ^ mn n n 


m=l 


m=0 


(D5) 

q 1 and V have been defined already in appendix C; where 
'n ’ mn mn J ^ ’ 
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7 < 2 ) 

"n 


2 “ / ,vm „ 2 2 2-6 

- (Sg) c il «VV O - = 2 (“) 1 (-T 32 ) 

1 mn 


n^O 


(D6) 


M 


i 1} -' K n- V n*X ] 

^ mn 


m=0 


■C^c^) -o (= ffi =g | = [B- - (-i)V”]) 


(D7) 


M m 2 2 

6 ( 2 ) = _k"-V" + V {T (— ) [ (C -G,)(^) - C_(— ) ] 

n n n £_, E mn v a ' 5 4'' 1 V 2 V a' 

mn 


m=0 


■ c 2 [Ci(5f) 2 - c (Sin £2 | si [b;- - (-u n b;*]) 


(D 8 ) 


Solving simultaneous equations (D5) for and B^, one obtains 


n 


B A = D 


D" 

nil n 

B n - D 

n 


where 


(D9) 


d - (^ 2) ) 2 - (r (l) ) 2 

n ' n ' w n ' 

M 


(DIO) 


M 


D’ = 7^6^ + 7 ^ Y (-1) I V [g* - C-l) n s” ] + W (-l)V [c'-(-l) m c M ] 

n n n n mn °m ' y °m n Z_ v mn n v n 


m=l 


m=0 


M 


M 


_ 7 ( 1 ) 6 ( 1 )_ 7 ( 1 )y’ Li^- (-lj^g” ]- Y v ' tc' - (-l) m c"] 
n n n L mn ti v n 4_. mn n v n 

(Dll) 


m=l 


m=0 


and 
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D ;=’'n 1)s n 2>+5 'n 1) £ (-lA^%(-l)^ : 


in=l 


M 


m=0 
M 


_ 7 (2) 6 (1)_ (2) y [£--(-l) n i"] - 7 (2) V V [3’-(-l) m c"] 
n n . n ^ mn v ' ^ n Z_ mn n > n 

(D12) 


m=l 


m=0 


Thus the unknown B^, have, in effect, through equations (D9), been 
supplanted by an equal number of known and K^. If the edges x=0 
and x=a are forced to remain straight, then the and are all zero. 


Formulation of Boundary Conditions of Equilibrium 
The normal forces acting on the rigid stiffener at x=0 must be self- 
equilibrating. Therefore equations (Cl4) and (C15) of appendix C apply 
also in the present case. They are; 


N 

P,(0) = - y — B’ + 
3 nit n 

n=l 

N 

P, (0) = y (-l) n B' 
4^ ' £_ nit n 

n=l 



(D13) 


(Dl4) 


Similarly, the normal forces acting on the rigid stiffener at x=a must 
be self- equilibrating, i.e. 

b 

? 5 (a) + P^(a) + J N x (a,y) dy = T 2 
0 

} T b 

P^(a) * b + J y N x (a,y) dy = — + M 2 

0 
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By substituting 

N 

N x (a,y) = B n sin 
n=l 


and solving for P^£a) and P^Xa), one obtains 


N 


T M_ 


P (a) = - V B" + -g - -g 
3 L nit n 2 b 


n=l 


N 


*4<«> " I h B n 


T M 
2 2 

2 b 


n=l 


(D15) 


(Dl6) 


Thus, through (D13) - (Dl6), the unknowns P^(0), P^(0), P^(a) and P^(a) 
have been expressed in terms of the knowns T^, M , T^ and M^. 


Separating B' and B" Terms in R ' , R", R'", R"" 

— £ - “ — n n n r m— *■ — in— 

Equations (D13), (Dl4), (D15), (Dl6), and (B69) can be used to 

eliminate P (0), P^(0), P (a), P^(a), andK^ from equations (b 68). If 

the B' and B" terms are then written separately from the rest, 
n n 

equations (b 68) become 


R' = S 


(5) 


7'B' 

n n 


"B- 

n n 


R" = 

s< 6 ) . 

- 7 m b ' + 

7>"B" 


n 

n 

n n 

n n 




N 

N 


R 

= s< 7) 

- y 3 

’ h ~y 

(-l) m B"H 

m 

m 

L , 

n mn L — ■ 

v n mn 



n=l 

n=l 




N 


N 

R"" 

m 

= s< 8 > 

m 

- 5 > 

ljHB'H + 

' n mn 

£(-i) m+n 



n=l 


n=l 


(BIT) 
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where are completely known quantities defined as 

n n m m 

follows: 

M 

= O' + |[P..(0) - (-l) n P-,(b)] + A E ~ T* - Y 
n ii b 1 ' x ' 1 ^ ' 1 1 b n / . mn 


(Dl8) 


X’A 

s n 6) - + |[p 2 fo)-(-i) n P 2 (»)] + Vs T T A + It- 1 )" s m ) 


(D 1 9) 


= O'" + t~ — (B 1 " - B"") + A,E — (C-B’" + T"’) 
m th bnut m m 3 3 a ' 5 m m 


-I 


T 2M T 2M 

S + — r - ('l) m t— h (D20) 

mn a ab v a ab 


£" - - B"") + A, E, — (C,B" " + T" M ) 

in bmir' m m 4 4 a 3 m m 


T, 2M 


T_ 2*L 


V (-l) n S^ 2 ^ + — + — ^ - (-l) m [— + — ~] (D21) 

v mn a ab ' a ab v 


( TT T mn - !<■-/ W' (-D n B;"]) 
mn 


(D22) 


H , 7 ', and 7 " have been defined already in appendix C: 7 "' are also 
nur n n n 

known quantities and are defined by the following equation: 


7'" = — + A E ^ C + V 
n ann 2 2 b 3 


1 2 /nrtv _ 

E a ^ b ; °1 
mn 


(D23) 


Revision of Equations (B 62 ) to (B 63 ) 


Substituting equations (D17) into equations \b 62) to (B 65 )., one 


obtains 


I 


c'CiL (n)-c"P n (n)=S^^+7 , B I - 7" B" - |(^) 2 Y ^4 — 
n r y nH_ v ' n nn nn b v b y 4_, E 

, mn 


m=l 

(n=l,2,3, . . .,N) 


(D24) 


-c'P 1 (n) + c"a Q (n) = s^-7"B' + 7’"B"+f (~) Y (-l) m ^4 — 
n K l v ' n 2 V 1 n n n n n b ' b' L-. v ' E 

n mn 

m=l 


(n=l,2,3...,N) 
N N 


n=l 


^3 (m) -^2 (B) - S m 7) + Z B iV-Z ( - 1)mB n^n - f(“) 2 I - 

n=l n=l 

(io=l,2 r 3,..., M) 

N N 

-ifzWefaW - S< 8) -£ (-1)^ + Y. ( - 1)mtn B n H ,»n 


(D25) 


H 5'-(-i) m c" 

' 1 n 


E 


mn 


(D26) 


n=l 


N 


n=l 
s m- 


0 2 ^ c'-(-l)V' 

£(-1)"-^ 2 o-l, 2,5,... ,M) 


n=l 


mn 


(D27) 


Equations (D2^) to (D27) can be used to obtain a system of 
simultaneous equations in which c^, c", g^, are the only unknown 
coefficients, by eliminating and B^ with the aid of equations (D9). 
The resulting system of simultaneous equations is: 


M 


^ - K^ 2) - ^ + I ^ - <- 1 > n (n=l,2,...,K) 

m=1 (D28) 

-"^P ) + l $i k) - ^n 2> Q (n=l,2, . . .,N) 

(D29) 


m=l 
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Z ®p f “5«VZ if V Cr n 2> (C- 1 ) P ^- 1 )”) - 7< 1) (l+(-l)“* P )]) -Z^mp' 

P=1 n=l n p=l 


N 


N 


y S' r< 3) - V S" r (4) 

/ n mn n mn 


n=l 


n=l 


N 


- S m 7) + Z ~!P C7 n 2)6 n 2) ' 

(m=l,2, ...,M) (D30) 


n=l 


-Z^ r Mp )+ Z «;KWvZ if^n t7^ ) ((-l) P + (-l) m )-7^ 1) (l + (-l) m+I ')]) 

~l 1 n XI 


P=1 P=1 


n=l 


N 


N 


-V (-i) n c' r^ 5) + V (-i) n c" r 
Z_ n mn Z_ n 1 


00 

mn 


n=l 


n=l 


N 


* S m 8) -Z ( ' 1)n ^ t7 n 2>6 n 2) ' 7^ 1) B' 1) -'(-l) m 7^ 1, 6' a ) + (-l) m 7( 2) 8( 1) ] 
n 

(m=l,2, . . . ,M) (D31) 


n=l 


where through 0^, £^ 2 \ ^ 2 \ and 

r n & ^n * ^n ’ b n J ^mn 7 s mn * mp ’ mn * mn 


are 


defined by the following equations: 

M 


M M 

= q: (n) - — [ 7 1 y" y^^ (_i) m v ' + i— v r 

^n u l w D L n 'n 7 n 7 n J Z- K } mn D L7 n 7 n Vn J Z- mn 
n _ ^ n 


m=0 

M 


m=0 

M 


0 2 = Pn(n) + — (-l) m v’ - — [ 7 ' y( 2 )-y"y(^ ] 

*n P i V ; D 7 n'n 7 n 7 n J Z- ' ' V mn D L7 n 7 n 7 n 7 n J l—. 
n n 


mn 


n=0 


m-0 

(D32) 
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0' 3) - P x (n) - jr^nf ’’Cf V in + 


mn 


m-0 
M 


m=0 

M 


= C£ (n) + 77- [ 7 " 7 (^_ 7 i «y(2) ■jV"' _ }^.[y"y(2) y, tt ~(l) -iV"* 

*n 2 V 1 D n L Vn f n 'n J Z_ ^ L) V mn D^Vn 7 n 7 n J Z- 


mn 


m=0 


m=0 


.CD 


( 2 ) 


(D32) 


t^L S (5) S — [y, 7 (l).y. 7 (2), + JL_ r> ,J2) 7 „Jl)-| 

~ b n D L7 n 7 n Vn J+ D L Vn Vn J 
n n 

B (l) r (2) 

+ — — [ 7 " 7 ( 1 )_ 7 ,n,y(2)-| n — [- ,. (2) _ , M (l) -J 

b n n D 'n'n 'n 7 n J D L7 n'n 7 n 'n •* 

n n 

(D33) 


mVi 2 V,L „ V^V^ V?V' 

d “~ ( - 1) D 
n n 


t (j-) = r -M m ' n, n ’ mn 'n'n ^mn , 'n'n r nm , 'n'n r mn 2,mt^ 1 

hen v ' D 


n 


- -(— ) 
b\ h ' 


D b ' b ' E 

n mn 


K ".y(2) M t -yl’-vC 1 ),,! 

^ - (-V D 


7 " 7 '"'n' 7 " 7 ^ 7 li 1 7' T y»»»y(2) , ? , m 

m n n mn 'n'n__fmn , ,m 'n n mn 'n 'n ^mn 2, mu (-1) 

n D _ ^‘ X; D D ’ b^ V E 

n n n n mn 


(D3*0 


r mp )=e 2 (m) VZ ( ' 1)n ^ “in ^n 2)[ (- 1 ) P+ <- 1 ) m] - f’ [l + (- 1 > m+P ” 

n 


n=l 


2 H M M 

- !<”> if - if f^ 2) X [(-D p + (-n m ] v' n .f ) y . v p n ] 


p=0 


p=0 


2 H M M 

= (- 1 ) 1 ” !<“) sr- - r 2 ’X >£[(-i) p + (-i) m ] v- ) 

mn n ^ F ^ 


p=0 


p=0 


(D35) 

Some uncoupling of unknowns in equations (D30) and (D3l) can be 
effected by adding and subtracting these two equations corresponding to 
the same value of m. By adding, one obtains 
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Sy s ^ tiag 
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pU) = s(^)+s( 8 )+2 y [^ 2 ) B ( 2 )_ 7 ( 1 )5( 1 )_^_i) m 7 ( 1 ) 6 ( 2 ) + (_i) m y( 2 ) & ( 1 ) ] 

m m m D n n n n ^ ■ n n ' ' c n n 


, 2 n 

n= 1, 3 • • • 


If 


(2) = S (7) -S (8) +2y r7 (2) 6 (2) -/ i;) & (l) -f-l) m 7 (l) 5 (2) +f-l) m 7 (2) & (l) l 

m m m L- D L 'n n / n n ^ ' 'n n ^ ' 'n n J 

n=2,4.?. 

(D4l) 

Equation (D 36 ) involves only the odd- subscript and c^, and equation 
^D37) only the even- sub script c^ and c^. These equations may replace 
equations (D30) and (D3l). 


m 


Reduction in the Number of Simultaneous Equations 
Equations (D28) and (D29), written for the same value of n, can be 

K and 


solved for c* and c" in terms of all the g' and g". The result is 
n n 


c' = (0^ V 1 ' ) -0^ 2) l^ 2) ) [£' - (“l) n g"j 

n 0 s n r n ^n r n Z_ ''n pn r n pn Lto p v B p J 

n P=1 

(D42) 

M 

> . i rK (lU(3) +? (2U(l) t y'( 0 (J) E (l) ,(1) jnj,,], 

n O' ls n ^n s n 'n [_ V5 n 5 pn & pn ' & p v ' & p 


p=l 


(D43) 


where 0’ = 0^0^ 

n r n r n r n r n 


Utilizing equations (D42) and (D43) to eliminate the c^ and c^ in 
equations (D 36 ) and (D37) and combining like terms, one can obtain two 
sets of simultaneous equations involving only the g^ and g^ as unknowns: 


M 


N 


M 


N 


N 


y g , {e^+y' f ] +y g" {©^+y t ■ > = p 1 '- 1 ^ - y 

/_ & p 1 mp mnp Z_ P m P l— mnp m 4_ 

p=l n=l, 3 ... p=l n=l, 3- • • 


* 


(m=l, 2, . . . ,M) 


mn 

n= 1, 3 • • • 
(D49) 
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M 


N 


M 


N 


N 


X ~Y "Crm 5 = p m^ " Y ♦m 

l — ■ p mp L — ■ mnp / . p mp L mnp m / ■ mn 

p=l n=2,4... p=l n=2,4 ... n=2,4... 


Cm=l,2, . . .,M) 


where 


V = ^ {sa (‘»)|(i)_ 0 C2)| Ca)) . &0(3) ( d) . 0 <o-) (2) } 

mnp <|>^ vr n pn r n pn <I>^ r n pn ~n pn 

(D50) 


(D49) ) 


, . 3s5L( t <i)/3> + s (2)dU>) 

T mn <J> ' '^n r n r n ' 4>' V1> n r n 3 n r n 1 

n n 

(D51) 

Whereas the original simultaneous equations system, equations (D28) 
to (D?l), requires the solution of 2N+2M simultaneous equations, the 
reduced, system, equations (D49), contains only 2M simultaneous equations. 
Thus N may be taken arbitrarily large without increasing the number of 
simultaneous equations that have to be solved. 


Procedure for Use of Equations 
The procedure for using the foregoing analysis will now be 
summarized: Equations (D49) are first solved for the g^ and g^. With 
these known, equations (D42) and (D4-3) give directly the c' and c", 


n 


n 


and equations ^D9) the and B^. Equations (b 60) then give the s^, 
s£, s^", s^", an< ^ equations (B57) bo (B59) and {B6l) the Finally, 

equations (Bl6) and (B19) to (B25) give the stiffener and plate stresses. 


Special Case: Symmetry About x = a/2 and y = b/2 

When the structure and loading are symmetrical about both centerlines 
x = a/2 and y = b/2, then A ± = Ag, = A^, P.^0) = P ] _(b) = Pg(o) = Pg(b), 
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P 3 (0) = P^(a) = P 4 (0) = P|^(a)> q^(x) = -q^(x), q-^y) = -q^y), e 1 (y) 


l, e ? 

(x) = e^(x). 

n 

= “ 

■V"; K 
n' i 

B’ = 

B" = 0 

for 

n 

even 

n 

n 




5' = 

c" = 0 

for 

n 

even 

n 

n 




= 

<=° 

for 

m 

even 

B' = 

B" 

for 

n 

odd 

n 

n 




c ' = 

c" 

for 

n 

odd 

n 

n 




= 

K 

for 

m 

odd 


> (D52) 


) 

We can simplify the simultaneous equations (d 49) as follows. From 
equations (D5) and (D52), one obtains 


M M 

B’ = -7T\ 1 ~7 7TT [ B (l) +2 7 d' g* + 2(Y v ' ) c'] 
n (1) (2) n mn'hn v mn' n 

'n 'n m=l,3*** m=l, 3... 


(tt=l,3,-..,N) 

Equations (D24) to (D27) are replaced by 


(D53) 


M 


c'Ca-CnJ-p-fn)] = S^ 5 ^+B ' (7'-7" ) - £ (^) Y 

n 1 1> n n v n n b v b' E 

^ mn 

m=l, 3* • ■ 
(i*=l,3,...,N) 

N 2 N 

n=l, 3* • • • n=l, 3- • • 

(m=l, 3> • ■ • jM) 

Substituting from (D53) into (D54), one obtains 


(D5*0 


c' 

n 

E 

mn 


(D55) 
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(n=l,3, ...,N) (D 56 ) 


M 


*s Jr- [£ + T" g' I ] 

n Q> b n fe p pn 

' n 


P pn‘ 
P=l,3- • • 


where 


M 


2 ( 7 * — 7 " ) , 

“n - « x (n) - 03 .Cn) - f!) J 2) Y. ^ 


y - y 

'n 'n m=l, 3* • • 


mn 


,(D 


So - s n 5) + (r A - >£) ,( 1 ) - .( 2 ) 

2 ( 7 ; - 7") 


y \ / _ y \ 

/ n 'n 


n 'n' _ 4 ,rw, 1 

b ' b ' E 


; P n 7 (l) _ 7 (2) pn 


n n 


pn 


(D56a) 


(D56b) 


(D 56 c ) 


Substituting from equations (D53) and. (D 56 ) into equations (D55)> one 
obtains 


M 

Z i ' (0 ' 

& p l mp 


P=l,3- ■ • 
(p odd) 



0" } 
mn 


£ 26^ 

’ S m + L “mn TTTU5) 

n=l, 3- • • n n 



(m=l,3, • • - ,M) 


(D57) 


0 " 

mn 


where 


mp 


[a 3 (m)-0 2 (m)]6 mp 


N 

Z 


H 


mn 


(1)_ (2) ^pn 

n=l,3... / n 'n 


(D57a) 


0 ” 

mn 



1 

E 

mn 


4h 

mn \ , 

’JTTjJT V p n 

n / n p=l,3--- 


(D57b) 


The number of equations in the system (D57) is (M+l)/2, regardless 
of the value selected for N. 

Equations (D57) can be solved simultaneously for g^ (p odd). With 
gp's as known, equations (D 56 ) and (D53) will give and directly 

(n odd). With these known, equations (b6o) will furnish the values of 
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s n f S n J7 S m ’ S m e( I lia ^ ions (®57) to (B59) and (B6l) the values of 
j . Equations (Bl6) and (B19) to (B25) give the stiffener and plate 
stresses . 


Limiting Cases of Large Stiffener Areas 
Five different sets of limiting conditions will now be considered, 
all for the case in which the structure and loading are symmetrical 
about the lines x = a/2 and y = b/2. In terms of the notations defined 
by equations (C80), these limiting conditions are; (l) X^— *-0, X^ finite; 

(2) X^-^O, followed by X — M); (3) X.p*-0, X g finite; ty) X:p*-0, 
followed by X^ — >-0; (5) X^ = X^ = X, followed by X— *-0. The reduction 
of the general equations to these limiting cases results in considerable 
simplification; in particular it is no longer necessary to solve simultaneous 
equations, except for condition (3)* Conditions (2), (4) and (5) which 
are physically identical also turn out to be mathematically identical. 

The details of the reduction follow. 


Condition (l): X^ — >-0 while X ^ maintains a finite value . - Use 

equations (D53) to eliminate the in equations (D55)j and then divide 
equations (D55) by A^E,., and rearrange them to obtain 


M V' Z" 

+ it!-) = --?- + 


I 


V" 


c' 

n 


qp- 4, mp A,E, A_E, mn A,E, 

5 5 D 3 2 J J 

P-1,3... n-1,3... 


N 

I 


where 


r- N 


V' = 6 
mp mp 


4 -(nit/b) 

b E 

L n=l, 3 . • • 


N 

£ 


H 


H ' 

pn 


ran (1) (2) 

J n=l, 3. . . 'n 7 n 


(m=l,3...,M) 

(D58) 

(D58a) 
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I 


z" 

m 


= S (T) 


in 


N 

♦■I 

n=l,3- 


H 


mn 


6 n l} 

7 ^ - 7 ^ 2) 


(D58b ) 


V" 

mn 


N 

2H ran (1) (2) X! V pn 

7 n 'n p=l,3... 



(mjt/a) 2 

E 

mn 


(D58c) 


Substituting from equation (D 56 ) into equation ^> 58 ) to eliminate V , 
and then examining the coefficients of the unknown g^ in equation (D 58 ), 
it is observed that some of these coefficients are of the order of 1, 
while the others are of the order of l^a^E-^A^E^) . Retaining only- 
terms of the order of 1 in these coefficients, one reduces equations 
(D 58 ) to the following system; 



m 


A 3 E 3 



n-1,3. • • 


n 


A 3 E 3 


(m = 1,3, 

(D59) 


Substituting from equation (D59) into equation (D 56 ), one obtains the 
following equation for c^ 


c ' 
n 



M 

«„ + I V 

p=l,3. . . 


■ z" 

E_ 

VS 


Vj 


N 

I 


V" 

jga 


T, 


-]) 


q=i,3--- 


(d6o) 


(n=l, 3, ■ • - ,N) 

Thus a solution (eqs. (D59) and (D 60 )) is obtained which gives the 
unknowns and c^ explicitly without the necessity of solving simultane- 
ous equations. This solution is seen to be correct to terms of the first 
in l/(a^E^A^E^) . With c^ and g^ known, then equation (D53) furnishes 
the values of B^, and the procedure for computing stresses is the same 
as described earlier for the special case of symmetry about x = a/2, 


y = b/ 2 . 
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Condition (2): *~Q with finite, followed by »-0 . - Equations 

(D59) and (d 6 o) can be used as a starting point for this case. It is 
first noted, from equations J^T>^6a) , (B 67 ), (C22), and (C23) that 


K = A 1 E 1 [1 " E n + A^; ] 


(D 6 l) 


where 



(D 6 la) 


,(D _ it 


n 


M 

I 


(nut/a)* 


mn 


M M 

TTTTI) 1 1 f!Q V 

n n m=l, 3- • • p=l,3- • • 


m=l, 3 - • • 

(D 6 lb) 

Substituting from equation (D 6 l) into equations (D59) to eliminate 0 , 
examining the coefficients of the unknown g^ in these equations, and 
retaining only terms of the order of 1 in these coefficients, one reduces 
equations (D59) to the following system: 


Z" 

- , m 

g m A E, 
3 3 


{m = 1, 3, M) (D62) 


A similar reduction of equation (D 60 ) can be effected with the aid of 
the fo 
(C23): 


the following form of | , obtained from equations (D56c), (C22), and 


e _ a E f?-- c II ' — 7 + i 1 

Sn A 1 E 1 L b C 3 V (1) (2) A E pn J 

7 n 7 n 


(D 63) 


where 
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a 


M 


/\(2) ^ r \ 1 2 f TUT \ n *T . ^ f liTt \ 1 

pn " Jl) (2 ) C Z- E a ( b } C l ] V " b^ b } E 


7 ' - 7 ' ' “mn 

n n m=l, 3 


pn 

(D63a) 


Substituting from equations (D6l) and (D 63 ) into equation (D60), examining 
the coefficients of the unknown c^, and retaining only terms of the order 
of 1 in these coefficients, one reduces equation (D60) to the following 
equation: 

Z’ 


c = 


n 


n A 1 E 1 


(n = .1, N) 


(D64) 


where 


3 >. 


n 


71 _ f rf . f f HI r 2 \ ,, 1 7” 1 

Z n L5 n A 3 E 3 b c 3 ,(!)_, (2) L ^rm V 


M 


7 ' -7' 

n 'n m=l, 3. 


(D64a) 


with 


< 1 >. 


,(D = 1 - E" 

n n 


(D64b ) 


Equations (D 62 ) and (D64) are the pertinent results for this limiting 
condition. 

Condition ( 3 ): while 7.^ maintains a finite value . - This limit- 

ing case (A^E^ and A^E^ approaching infinity) can be studied with equations 
(D57) as the starting point. These equations are written in the following 
form, in order to be more readily suitable for solution by the Gauss-Seidel 
iteration method: 


■ >re> + 
mm 
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CD 
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i * n 

n=l, 3* • • 


fj /_•> I — ' 26 ^ ' r £ 

mn _ CT) , V w n V In 
mn m 2_ mn -,CD_.y(2) Z_ ^ mn 


7 '• ' _ y • / n 

n=l, 3- • • n 'n n=l,3... 


M 
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Z w pn - Z 

n=l, 3- • ♦ 


mn 

y (D~ (27 V 

P=l,3- • • n=l, 3 --n 'n 
(p / m) 

(m = 1, 3, M) (D 63 ) 
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Substituting from equations (D 6 l) and (D 63 ) into (D 65 ), then examining 
the coefficients of the unknown g^ in equations (D 65 ), and retaining 
only terms of the order of 1 in these coefficients as A^E-^ and A^E^ 
approach infinity, one reduces equations (D 65 ) to the following 
system: 


K (0 in + 2 


N 

Z 

n=l, 3 ■ 


(nn/b) C 0 ^ 

[1 _ E"][ 7 (l) - 


n n 
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- s « ♦ 7 

D1 


26 
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H 
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mn (!) J 2 T ' Z A E [1 - Z " ] 
n=l, 3* • • n n n=l,3... 


i 8" 
n mn 
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N 


I K (k I 


H pi. ’ 

mn pn ^ 


-P A— 7 (D _J2) 

p=l, 3 ... n- 1 , 3 ...n 'n 

(p f m) 


N 

I 


(nit/b)C^ ^ 6^ 

[1 - H"][7 ( 1 ) -7< 2) ] 


n= 1 , 3 • • • n n ' n 
(m = 1, 3, M) (D66) 


The solution of equations (d 66 ) can be effected by the Gauss-Seidel 
iterative procedure (see appendix C - Procedure for numerical solution.)* 
Using equations (D53) to eliminate in equations (D5^ ) , and then 
dividing equations (D 5 ^) by A^E^, and rearranging terms, one obtains 


m 

(D 67 ) 

in which the are defined by equations (D 66 ). As >-0 (A^E !— ), 
equations (D 67 ) are simplified to the following: 
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A 1 E 1 


nn 
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■{ay Z 

n m= 1 , 3 < 
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mn 
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mn - 
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Z A E 

= 1 , 3 ... 


J m 
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- , _ 1 r n . nn n 2 

C n ~ 7 ( 1 ) L A 1 E 1 b C 3 1102 ) 


' 11 
n 


M 

I 


m=l, 5* • • 


LL ' g ' ] 

K mn & m 


(D68) 


where is defined by equation (f 64 b) 


With g^ and known, through equations (D 66) and (D68), equations 
(D 53 ) will give the values of B^. From that point on, the procedure for 
computing stresses is the same as described earlier for the special case 
of symmetry about x = a/2 and y = b/2. 


Condition ( 4 ): — >-Q with \ finite, followed by — >- 0 . - Dividing 

equations (D 66) by A^E^ and examining the coefficients of the unknowns 
g^, it is observed that these coefficients contain terms of the order of 

z 

1 and terms of the order of l/ ( a^E^A^E^ ) . Retaining only the terms of 
the order of 1, one reduces these equations to the following system: 



Z" 

m 

A 3 E 5 


(D69) 


where Z^ is defined by equation (D58b). Substituting from equation (D69) 
into equation (D68), one obtains 


c ' 
n 


Z' 

n 

A 1 E 1 


(D 70 ) 


where Z^ is defined by equation (D6Ua). 

With g^ and c^ known from equations (D69) and (D70), equation (D 53 ) 
will then furnish the values of B^. The subsequent procedure for comput- 
ing stresses is the same as described earlier for the special case of 
symmetry about x = a/2 and y = b/2. 
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Condition ($): = A, A— »-0 . - This case can be thought of 

as one in which all stiffener cross-sectional areas are approaching 
infinity simultaneously while maintaining fixed ratios with respect to 
each other. For the study of this case we return to equations (D54) 
and (D55), but rewritten as equations (D 67 ) and respectively. 

Examining the coefficients of c^ and g^ in equations (D 67 ) and (D 58 ), 
it is observed that some of these coefficients are of the order of 1 , 
while the others are of the order of l/(a^E^A^E^) . Retaining only 
terms of the order of 1 in these coefficients, one reduces equation 
(D 67 ) and (D 58 ) to the following system: 


c ' 
n 


n 





M 

I * 

m=l,3- 


mn 



(D71) 


g' 


Z" 
m 

'm ” A_,E, 


(D72) 


Using equation (D72) to simplify equation (D7l).> the latter becomes 


Z’ 

c’ = 

n A i E x 


(D73) 


where Z^ is defined by equation (D64a). 

Thus a solution (eqs. (D 72 ) and (D73)) is obtained which gives the 
basic unknowns c^ and g^ explicitly without the necessity of solving 
simultaneous equations. This solution is seen to be identical to those 
obtained for conditions ( 2 ) and (4). 
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Illustrative Thermal-Stress Problem 


A particular example will be presented to illustrate the details 
involved in the application of the foregoing analytical results, this 
example has the following characteristics: 

a) Two opposite edges of x = 0 and x = a kept straight; therefore 
the and in equations (8) and ( 9 ) are all zero. 

b) Plate isotropic; therefore elastic constants are given by 
equations (3) • 

c) Plate and stiffeners have the same Young's modulus E. 

d) = A^, A^ = A^ . 

e) No force loading. 

f) Stiffener temperature constant at the value T . 

g) Plate temperatures T(x, y) symmetrical about both centerlines 
(x = a/2, y = b/2) and varying sinusoidally in accordance with the 
following equation 

*d r, (0 < x < a) 

T(x,y) = T +9 sin (— — ) sin (— jp) 

° a b (0<y<b) 

where P and Q are odd integers. 

h) Plate and stiffeners have the same coefficient of expansion a. 

These are the only specializations to be employed at present. However, 
in the numerical example, to be presented later, the problem will be 
further specialized to the case of a square plate (b=a) with all stiffeners 
identical (A^ = A^ = A^ = A^ = A), subjected to a "pillow- shaped” temper- 
ature distribution (P=Qr=l), and having v = 0.3* 

R educ t ion of general equations to special case . - From the given 
temperature distribution, one obtains the following equations for the 
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known coefficients in terms of temperature distribution and the 
coefficients of expansion a (see appendix C): 


mi _ qiH _ rp t »» _ rp'ltl _ Q 

n n m m 


(D7M 


T 

mn 


8 mP & nQ 


CtG 



[P 2 + Q 2 B 2 ] 


(D75) 


V’ = -V" = a© — 5 _ (D 76 ) 

n n a nQ, 

Due to the absence of prescribed forces, the following quantities are 
all zero: 

P 1 (0), P x (b), P 2 (0), P 2 (b) (fig. 4b) 

T x , T 2 , M 2 (fig. 4b) 

B;", B”" (see eqs. (4)) 

<%" (see eqs. (6)) 


(it should be noted that P^(0), P^(a), P^(0), P^(a), and B^ do not 
necessarily vanish. ) 

Because in this example the structure and loading are symmetrical 

about both centerlines, x = a/2 and y = b/2, the simplified system of 

equations, namely equations (D 5 6 ) and (D57) will be used for the 

determination of the c 1 and g 1 with odd m and n. (it should be noted 

that c’ = c" for n odd, g’ = g" for m odd, c’ = c" = 0 for n even, 

n n m m n n 7 

and g’ = g" = 0 for m even). With c' and g' known, equation (D53) 
will furnish the values of odd subscripted B^. (it should be noted 
that B^ = for n odd, and B^ = B^ = 0 for n even). Equations (D57) 

must first be solved simultaneously for g. The quantities needed in 
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order to use these equations will now be evaluated. Substituting from 
equations (3) and (C8l) into equation (D 6), one obtains 



-22/ V-* C- 1 ) 1 ^ 2 + (2+v)n 2 B 2 ](2-6 ) 

(n^B / aEh) ) ^ 

i=0 [i 2 + n 2 B 2 ] 


(D6 1 ) 


Prom equations (CIO') and (d 6')> one obtains 


n n 


aEh 


M 

z 


(D77) 


[m g + (2+y) n 2 B 2 ] 

r 2 . 2^2 
Lm + n B ] 


m=l, 3- • 


Substituting equations (3), (C8l), (D75) and (D76) into equation (D?), 
one obtains 


s u> . os si rh ? tJazMgfl . x , 

n a 1 r _2 . _2_2, ‘ 


[P 2 + Q 2 B 2 ] 


nQ 


(D7') 


Substituting into equations (Dl8) and (D20) from equations (D7^) and 
(D75), one obtains 



6 . a0Eh 
nQ 


PQB 

[P 2 + q 2 b 2 ] 


(Dl8' ) 



5 a9Eh 
mP 


PQB 

[P 2 + Q 2 B 2 ] 


(D20 ' ) 


Substituting from equations (B67), (D77), (C22'), (C23') and (C13'), 
into equation (D56a), one obtains 


0 = A,E 0 
r n 1 r n 

where 


(D56'a) 


133 



M 


\ - 1 + \ I 


m 


m=l, 5? • 


r 2 . 2^2 ,2 

Lm + n B J 


M 


[v UjllB Y , 2 2 2, 2 1 M 

* nE l f 

nt=l,3. 


M 

I 


[m 2 +(2+v)n 2 B 2 ] 
Cm 2 + n 2 B 2 ] 2 


[n^ 2 - vm 2 ] 
[m 2 + n 2 B 2 ] 2 

(D^6 M a) 


Substituting from equations (D77), (Dl8'), (C22')> j[C23') and (D7' ) 
into equation (D 56 b), one obtains 


z 


n 


CtGEh t, 6 . 

n nQ, 


(D 56 ’b) 


where C n are known quantities defined as follows 


M 


Fl b + ^ Y 




FOB 


p O O 

[P +Q B ] 




1 l r [P 2 +(2+v)Q 2 B 2 ] _ ? 

(i 2 +Q 2 B 2 ) 2 [P 2 + Q 2 B 2 ] 


M 

SB £ 




[i 2 + (2+v) Q, 2 B 2 ] 
[i 2 + q 2 b 2 ] 2 


(D 56 "b) 


Substituting from equations (C8l), (C12'), (C22’), (C23') and (D77) 
into equation (D 56 c), one obtains 


| = A E£ 

pn 1 pn 


(D56'c) 


where £^ n are known quantities as expressed by the following equations 
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1 


M 


B[p 2 +(2+v) n 2 B 2 ] [v+X^B 2 ^ 




/ . 2 . 2 . 2 x 2 ' 

(l + n B ) 


X x nV 


M o P ? 

(p 2 ♦ n 2 B 2 ) 2 f [i ^ | “» ] 
._T% Ci +n B ] 

l— x, p . • . 


/ 2 . 2.2 v2 

(P + n B ) 


(D56"c) 


Substituting from equations (c82), (C24 1 ), (C12*) and (D77) into equation 
(D57a)> one obtains 


N 


e mp - a jE (U * V 4 I (m 2 7 - 

trt — l -S N ' 


n=l, 3* • • 


N 


V 


*■/ 2, 2.2 x 2 " 2.4 ][ P + ^- 2+V ^ n B ^ 


' I 

n=lj 3 • • 


(m +n B ; n B 


M 


5 (D57'a) 


(p2 + nV)2 £ 

1=1,3-.. 


2.2- 


r.2. 2.2 ,2 

Li +n B J 


Substituting from equations (C8l), (C13'), (C24')> and (D77) into 
equation (D57b), one obtains 


0" = A..E 0" 

mn 1 mn 


(D57'b) 


where 0" are known quantities and are defined as follows 
mn u 


nV 


0 " = 




- ^3 


[m +n B T 


2 2.2 v 2 nB 

(m +n B ) 


M 

nE 

1=1,3- 


M 

I 


r . 2 , 2.2 -.2 

Li +n B J 


[n^B 2 


vi 2 ] 


[i 2 +(2+v)n 2 B 2 ] 1=1, 3 . • • 


[i 2 +n 2 B 2 ] 2 


(D57”b) 
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Substituting from equations (C24')> (D7'), (D77), (D20'), (D56'a), (D56’b), 

(D5 6'c), (D^7'a), and. (D57't ) into equation (D57)> and substituting 

g' = g' C.. one obtains 
P B P 1 or 


Nr- N 

y _»s 5 .. + tl t x y y 

77 mn 2 / . 


n=l, 5 . .^n 


n=l,3... 


n 

/ 2 . 2 2 v 2 

(m + n B ) 


] C 


N 

£ 

n=l, 3 • • • 


n 


(m 2 +n 2 B 2 ) 2 


M 


(m 2 +n 2 B 2 ) 2 ^ 

i=lj 3 * 


-|-^][m 2 + (2+v) n 2 B 2 ] 
n B 


[i 2 + (2+v) n 2 B 2 ] 

r . 2 , 2 2 n 2 

Li + n B J 


G(m) 


- 

" & mP ),fr£ rt 2-.2x “ 4m 5 nQ 77 6 mn 


it A. Q B 


4(P + QB) 




n 


nV 


P[- 


[P 2 +(2+v)Q 2 B 2 ] 

3 b 3 / 2^2 D 2,2 Jl r*2 , „2„2- 


QB 


^ _Jl Q^B (m +QyB ) [P" + Q B ] 

4m M 

[i 2 + (2+v) Q 2 B 2 ] 


~ 1 ) 


I 


1=1,3-. 


[i 2 + q 2 b 2 ] 2 


M 


I 

P=l,3« 


N 


^ 6 . [7XT33 - -ixlC-^vJnV] 

[ > ( -ES 8" . x b 4 o ] — g_§ 

• 1-1,3... t s + n 2 B 2,2 V [1 +( 2 +v)n g B I 

1.1,5 [i+nB 1 


)] 


(1-6 ) £ G(p) 

' mp m 


(D57’) 
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where 


aGEh 


and 


c = a 3 /a x 


In the solution of this system of equations the Gauss-Seidel 
iterative procedure can be employed (see appendix C - Procedure for 
numerical solution). 

After equations (D57') have been solved, the G(p) (therefore g^) 
will be known. The c^ can then be determined from equation (D 56 ). 
Substituting from equations (D56'a), (D56'b) and (D56 ’c) into equation 
(D 56 ), and substituting c^ = c^ one obtains 


C(n) = 


jr K 


M 


1 v 


4(? nB 
r n 


p=l, 3- • • 


P I 


pn 


(D56’) 


where 


a© Eh 


With G(n) and G(m) known, substituting from equations (D77), 

(DY' ), (C12')> and (C 13 1 ) into equation (D 53 ), and substituting 

c’ = c' C and g' = g' C. — , one obtains 
n n 2 b to m ^ 1 a' 

2 , 2^2 


M 


B(n) = 


p - 1 ) 
[p 2 +q 2 b 2 3 


4q 2 B" 


M ? 2 nQ M 

2V [i 2 +(2+v)Q. 2 B 2 ] t"* 

^ [ i 2 +Q 2 B 2 ] 2 L 


V ^ * 2 (2 7l\ B 1 G(m) 

.7, lm 2 +n 2 S 2 f 

6 + SSl i li 


i=l, 3* • • 


M 


. m=l, 3 - • • 


[n B - vm 3 
" 2 b 2 3 2 


V in B - 
“ , [m 2 + h 


i=l,3- •• 


C(n) 


[i 2 +(2+v)n 2 B 2 ] 

r .2 . 2_2 - i 2 

Li + n B J 


M 

I 

i=l, 3* • • 


[ i 2 + (2+v ) n 2 B 2 ] 

[i 2 + q 2 b 2 3 2 


(D53’) 
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where 


B(n) 


B’ 

n 

a© Eh 


With B’ as known, one obtains 
n 3 


N 


P 5 (0) = P 4 (0) = -O0A 3 EX 2 £ ^ B(n) 

n=l, 3. • . 


(D78) 


from equations (D13) and (Dlk). 

With C(n), G(m) and B(n) now known, equations (b6o) will give the 
values of the odd- subscripted s^, s^, s^", s^", and equation (B6l) the 
values of j (m and n odd). Equations (Bl6) and (B19) to (B25) will 
give the stiffeners and plate stresses. One thus obtains the following 
results: 


N 

P^y) = P 2 (y) = ^ Cc(n) - vB(n)] sin (^) 

n=l, 3* • • 


(0 < y < b) (D79) 

M 

P^(x) = P 4 (x) = QA^Ea ^ G(m) sin (2~) (0 < x < a) 

m= 1, 3 • • • (D80) 


M 


N 


N = a0Eh 
x 


Z Z sin sin c^) 

10=1, 3 . • • n=l, 3- • • 

(0 < x < a) 

(0 < y < b) 


(D8l) 


M 


(N^o = (N x } y-b = a0Eh Z G(m) Sin (I5 ? ) ^ < x < a) 


m=l,3. . . 


(D82) 
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M 


N 


N = o9Eh 

y 


Y c ( m > n ) sin X 2 —) sin (^) 

m=l, 5 . . . n=l, 3 . . . 

(0 < x < a) 

(0 < y < b) 


(D83) 


N 


^V *=0 = ^y^a = Y C ( n ) Sin ^b^ (0 < y < b) 


n=l,3. . . 


XD 84) 


M 


W 


N = -OJ0Eh 
xy 


Y Y j ( m > n ) c ° s (^in) c ° s (^) 

m=l, 3- • • n=l,3- • • 

(0 < x < a) 

(0 < y < b) 


(D85) 


where G(m,n), C(m,n), and j(m,n) are known quantities, and are defined 
by the following equations: 


kn 2 -R 2 o 

G(m,n) = "p - 3 " 5 [m C(n) + nB 2 G(m)] 

*[m 2 +n 2 B 2 ] 2 


4m [m 2 + 2n 2 B 2 ]B ( n ) 

r 2 . 2 2 -,2 

irLm + n B J 


2 2 
QB 


5 „8 


(P 2 + q 2 b 2 ) mP nQ 


(d86) 


C(m,„) - * *(f) - 1 ] - 

n 2 B 2 ir (m 2 + n 2 B 2 ) 2 (P 2 +Q 2 B 2 ) ^ Q 


it [m 2 +n 2 B 2 ] 2 


[m C(n) + nB G(.m) ] 


(D87) 


J(m,n) = ^p-p - p fn 2 B 2 B(n) - m 2 C(n) - mnB 2 G(m) ) 

n[m +n B ] 


PQB 


8 „8 


(p 2 +q, 2 b 2 ) 


(d 88 ) 
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And from the first two of equation (ii), one obtains the following 
running tensions between the rigid stiffeners and plate edges at 
x = 0 and x = a. 

N N 

N x (0,y) = N x (a,y) = £ sin (^) = a©Eh £ B(n) sin (— ^) 

n=l, 3. . . n=l, 3- • • 

(0 < y < b) (D89) 

Numerical results for X and X^ 4 - 0 » _ The numerical procedure 
and equations described above were applied to the special case of a 
square plate (b=a), will all stiffeners identical (A^ = = A^), 

and Poisson's ratio v equal to 0.3. The above assumption implies B = 1 
and X = X^ , and the common symbol X will therefore be used for 
both X and X^. 

The results obtained for the stiffener tensions and the plate 
stresses are presented in dimensionless form in figure 10 for X = 2.0 
and figure 11 for X = l.o. The values of M and N employed in the 
calculation are indicated on the figures. In general, stresses are 
computed at x/a and y/a interval of 0.1. 

Limiting case of large stiffener areas, condition (l) . - Five 
different limiting conditions were considered in the earlier analysis. 
However numerical calculations were made for only two of these, conditions 
(l) and (2). The equations needed for condition (l) will now be presented. 

Substituting from equations (D7' ), (D20'), (D77) and (024') into 
(D58b), one obtains 
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II 


z" = a0Eh 

m 


r SB. 


« _ (m 2 +Q 2 B 2 ) 2 




(p 2 +q, 2 b 2 ) 


M 




[p c + q 2 b 2 ] 


^ £if •*• (2+v) Q 2 B 2 ] 
i=l,3... 


[i 2 + 


(D58'b) 


Substituting from equations (Cl}'), (C2V), (c8l) and (D?7) into (D58 c), 

one obtains the following equations for V" : 

mn 


V" = 
mn 


4aEh 


{[■ 


2 2 
n B 


(m 2 +n 2 B 2 ) 2 


M 

I 


[n 2 B 2 - vp 2 ] 


r 2 2_2-,2 

i j * nB ] 


2„2 

n B 


m 


M 


Y [p g tf^v)n a B a 1 ^ 2 B S f 
Z- r 2 , 2 P-,2 

p=l,3... [p + n B ] 

(D58'c) 


Then substituting from equations (P58'b), (D58 'c), (D 56 ’a)and (D56’b) into 
(D59), one obtains 


i; - oe S- G(m) 


(D59 1 


where the G-(m) are completely defined by the following equation: 


[ SB 1 ]f [P"+(2+v)Q 2 B 2 J _ . 

G( m ) = ^—6 - ( m5+ Q 2 B 2 ) 2 Q? E? [P 2 + Q, 2 B 2 1 

(pAq 2 B 2 ) mP M p T7 

Y [i 2 + (2+v ) Q 2 B 2 I 

Z_ r . 2 2 2 ,2 

i«l,3... [l +Q B ] 


M 


y ^ 2 b 2 - y p^3 

p p rA^ 2 i 2 

- M [ — S-Bf i ] ] 


1 VWW qV 


m 


V [p^+WqVI [m 2 +Q, 2 B 2 J a % 
Z^ r 2 . _2_2 ,2 

(D 59 'a) 


£ 0 

} 


P=l,3-«. 


[p 2 + Q 2 B 2 ] 2 


Q 
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with ^ and Cq are defined by equations (D56"a) and (D56"b). 

With g^ as known, substituting from equations (D59')> (D 56 ' a), 
(D56'b) and (D56 'c) into equation (D5 6 ), one obtains 


c * 
n 



M 


n 

— 5 n + ae 

nQ 


* 


n 


^ I 5(m) f 


mn 


(D56") 


Substituting from equations (D59' (D 56 " ), ^D7'), (D77)j (C12') and (C13') 
into (D53)> one then obtains 


B ’ = -Q0Eh V B 1 5 + O0Eh X. r~- B* 

n 1 4QB n nQ 2 4n n 


(D53” ) 


where 


M 


B' = 
n 


rf [p 2 +(2+v)Q 2 b 2 ] _ 1} 
[p 2 + q 2 b 2 ] 

M 


4 I 


m=l, 3 , • • • 


[Q, 2 B 2 - vm 2 ] 
[m 2 + Q 2 B 2 ] 2 


M M 

K cm V " [m 2 +(2+v)Q 2 B 2 ] * V [m 2 +(2+v)Q 2 B 2 ] 

1^ r 2 . „2^2 n r 2 , A2 i 


m=l,3. . 


[m + QB ] 


m=l,3* • 


[m + QB] 


(D53"a) 


B 


M 

I 


[m 2 +(2+v)n 2 B 2 ]G(m) 


M 


M 


m=l,3... 

n M 


r 2 2 2 -,2 

Lm +n B J 


[ Z G(p)l pn ] Z [n2B2 ~ vm2 ^ 




, , r 2^ 2 2,2 

m=l,3- . Lm +n B ] 


I 


m=l, 3 • • • 


[m 2 + (2+v)n 2 B 2 ] 

r 2 ^ 2 Q 2 n 2 
[m + n B J 


M 

m=l,3- 


[m 2 + (2+v) n 2 B 2 ] 

r 2 . 2 2,2 

Lm + n B ] 


(D53"b) 


With and B^ known through equations (D 56 "), (D59')j and 

(D53" ), the Fourier coefficients of the stresses can now be evaluated. 
From equations (b6o), one obtains 
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* r^n 


(D90) 


s 1 = s" — a0A,E 5 T-^ [— + VB’] 

n n 1 1 nQ 4QB 77 n 

^n 

s'" = s"" = O0A,E X Q I 21 G(m) 
m m 5 2 4m 

From equation (B6l), with c^, g^, and B^ as known, one obtains 

j = aGEh & j (D91) 

mn nQ, mn 

(m/0) 

Wo) 


where 


PQB 


m 2 \„ 


3 n 


J mn 


(P 2 + Q 2 B 2 ) ^ [m 2 + Q 2 B 2 ] 2 


M B B A 

[m 2 +Q 2 B 2 ] 2 


(D90a) 


With c ' , 
n 

Kt andB^ 

as known, 

obtains 

®mn 

6 nQ 

aGEh 

^irai 


c 

= 5 . 

aGEh 

cmn 

where 

mn 

nQ 

o o 



- 


^ B6 mP , 

X 1 B 


®mn 


4 1 b n 1“ 


n 


(P 2 +Q 2 B 2 ) $ 


2 tr ZJ2 ,2 


r 2 , „2 2 .2 

[m + Q B J 


(»9i) 


(D91a) 


- _ X l m r-, m 2 (m 2 +2Q 2 B 2 ) 1 P 2 

C mn " q 3 b 3 ' (P 2 + Q 2 B 2 ) 






n 


p p p p _ 

QB[m + Q B ] 0 


n 


(D91b) 


Mow with g^, s;", V W and too™,, equations 

(Bl6) and {B19) to (B25) will give the following stiffener and plate 
stress equations: 
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N 


P l^ y) s p 2^ y) = ^1^1 Y 5 nQ 2^ ^ + vB ^ sin ( B f L '> (0 < y < b) 


n=l, 3- 
M 


n 


(D92) 


P 3 (x) = P 4 (x) = 0© A ? E *, £ G(m) sin {£2*) (0 < x < a) 


m=l,3. . 


^93) 


M N 


N = QSEh 
x 


z 1 ^ «*“ <=? > ■*» <t > ;: 

rn=l,3... n=l,3... 


< y < b) 
(D9*0 


M 


(N ) „ = (N ) ^ = a©Eh / jt G(m) . ,mjtXs „ ^ ^ N 

x y=0 x y=b 2 i—> -r A sin ( ) = 0 (0 < x < a) 

-y ^ful Q. 


m=l,3- 


M K 


^ Y 6 c sin sin (2£Z) 

L-> t_j nQ mn ' a ' b ' 


N = O0Eh . , ~ rs 

y l—> nQ, 

®=1, 3 • • • n=l,3- • • 


N 




(N ) = (N ) = o©Eh 7 5 . r^- - sin <~Z) 

y x=0 ' y x=a nQ, 4QB ^ v b 


n=l,3-- • 


0 


n 


(D95) 

(.0 < x < a) 

(0 < y < b) 
(D96) 

(0 < y < b) 

(D97) 


N = -a©Eh 
xy 


M IT 

y 5 3 cos cos (^r^) 

L-, nQ u mn x a ' ^ b ' 

m=l, 3- • • n= 1,3... 


X'X 

I 


(0 < x < a) 
(0 < y < b) 

(D9 8) 


Substituting from equation (D53" ) into equations (D13) to (Dl6) and 


neglecting terms involving K^, one obtains 


\ »t 


P 3 (0) . P 3 (a) , p 4 (0) = P t (a) = aeA 5 E V, 5^-^- 


(D99) 
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With as known, from the first two equations (4), one obtains 


N.. jc __ _ 

N x (0,y) = N x (a,y) = -aGEhb^ ^ B^ sin (“) (0 < y < b) 

(D100) 


Numerical results for condition (l) . - The numerical result for 
= 1, — >"0 for square plate (B=l) with all stiffeners identical, 

subjected to a pillow- shaped temperature distribution j(p=Q=l) is shown 
in dimensionless form in figure 12 for Poisson's ratio v equals to 0.5* 


Limiting case of large stiffener areas, condition (2) . - In this 
case, equations (d62), (d 64), and (D53) will be the governing equations. 
The additional quantities needed in order to use these equations will 
now be evaluated. 

Substituting from equations (3), (C13'), and (D77) into equation 

(D64b), and with aid of equation (D6la), one obtains 

M 

Y [n 2 B 2 - vm 2 ] 

*(!) = ]_ + v m=13...[m 2 + A if ( D 64'b) 

n M P 2 o 

V [m + (2+v) nB ] 

" r 2 . 2^2,2 

[m + n B ] 

m=l, 3. . . 

Substituting from equation (D58'b) into (d 62), one obtains 




where 


(D62 ' ) 


2J2- 


G'(m) = 


PQB 


(/+® 2 b 2 ) mP 


r §g_ 

L O r 


f [P + (2+v) Q B g ] •, 


, 2. r ,2J2y2. 

(m +Q B ) 


[P 2 + q 2 b 2 ] 


qV 


M 


I 

i-1,3. 


[i 2 + (2+v) q 2 b 2 ] 
[i 2 + q 2 b 2 ] 2 

(d 62 ' a ) 
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Substituting from equations ( 3 ), (C12’), (D7'), (Dl 8 '), (D 77 ), (D56’b), 
(D58’b) and (D64'b) into equation (D64), one obtains 


2 c w 

c' = -ryfl ^ r /-yo ^ ) 

n 4a $ (l) S nQ 09 A ± ~(Tj C n 

n n 


(d64') 


where 


VP ( a 2 + frv) pti . x , 

.(l) _ [P 2 + q 2 b 2 ] 


M 

® Z 

m=l,3- • • 


[m 2 + (2+v) Q 2 B 2 ] 
[m 2 + Q 2 B 2 ] 2 


0d64") 


M 


-( 2 ) _ r PgB 
n “ L2 . „2^2, 


(P + Q B ) 


pqb [ i£+istvi^!i . l3 Z 

(P + Oil m=l. 3. . . , R 

M _ j & nQ 


2,^2 ,2 


V U + C2+v) qV] 

Jr 3 ... [» 2 ^ 2 b 2 ] 2 


M 

+ vB( £ 

”i=l, 3 



+ (2+v) n 2 B 2 ] 
[m 2 + n 2 B 2 ] 2 


M 

G(m)} /{ C ^ 

m=l, 3 ■ • • 


[m 2 + (2+v) n 2 B 2 1 

[■2 2 2 -,2 J 

Lm + n B J 


{D64’") 


With and g^ (m, n odd) known, and substituting from equations 
(D 63 ’ ), (D62 1 ) , (D7* ), (D77), (C 12 ' ) and (C13’) into equations (D53), 
there follows 


B' = G©Eh B^ 1 ^ 5 „ + aGEh K B^ 2 '* 
n n nQ In 


(D53'") 


where 
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M 


B 


(1) 


n 


4q 2 b 2 


c' 11 1 r l p 2 1 if " F[ [^v) . 13 

, m=1.5... [m + <* B ] . fp 2 + a. g B g l 

l w J 


M M 

[m 2 +(2+v)Q 2 B 2 3 V U 2 +(2 +v)Q 2 B 2 ] 


®(1> V U^(2 + v)q^3 y 

*n r 2 2 2 -,2 Z _ r 2 _2 2 -,2 

-12 Lra + Q B J m , [m +Q B J 
m=l,3... m=l, 3... 


M 


I [ f^t ~ G 'W *<£° I 7% 


M 


Bit 


B 


( 2 ) 


(D53'"a) 


2 2 2 n 
B -vm J 


[m c +n"B c ] c u [m"+ n 2 B 2 ] 2 

m=l, 3 « • . m=l, 3 • • • 


knC 


M 

I 


[m 2 + (2+v) n 2 B 2 j 


M 


m=l,3- • • 


r 2 . 2 2 -,2 

Lm + n B J 




m=l, 3, . . . 


[m 2 +(2+v)n 2 B 2 ] 

r 2 , 2 2 -.2 

Lm + n B J 


(D53"'b) 


With the basic quantities c^, and. B^ known (equations (D64 ' ) 

(D62* ), and (D53'"))j the additional quantities needed in order to deter- 
mine the stresses will now be evaluated. Equation (b6o), (b6i), (B3^), 

and (B35) will give the values of the odd- subscripted s' s", s'", s'",' 

% ° n’ n’ m J nr 


j . g , c as follows: 
mn & rwr mn 


s’ = s" = OSA [£- 
n n 1 1 4n 


< 6 . 


.( 2 ) 

Tiy VB n 2>: 


(D101) 


m 


= s 


m 


O0A^EA 2 


sS G '« 


(D102) 


J mn 

(V°) 

(n^O) 


O0Eh & „ J 

nQ mn 


(D103) 
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g = OC0Eh 5 _ G 
°mn nQ mn 


c = O© Eh 8 - C 
mn nQ mn 


(d 105 ) 


where 


m 2 C^ ^ 


4 B^ 


j _ — ±J3E 5 + £ — _ + £ 

11111 (P 2 +Q 2 B 2 ) B[m 2 +Q 2 B 2 ] 2 0^ *[m 2 + Q 2 B 2 ] 2 


(D103a) 


q 2 b 2 mQ, 4m[m 2 +2n 2 B 2 ]B^ 1 ^ 

= ' (pVb 2 ) 8mP ' ^[b^ n 2 B 2 ] 2 + „[„ 2 + n 2 B 2 ] 2 


(DIOUa) 




W B^ iy 2,2 22, m 

r n_ r m Cm + 2Q B J n _ 

mn ” .2 2 L , 2 . .2^2,2 XJ 2 2 r 2 7 A 2 B 2 n 2 _ (l) 

it QB (m+QB) Q B Lm + Q B J 


(D 105 a) 


(P 2 + q 2 b 2 ) ^ 


Now with c', g', B ' , s', s' , c , g , and j known, equations 
n m n' n m mn °mn mn 7 

(Bl6), and (B19) to (B25) will give the following stiffener tensions and 
the plate stresses: 

N c (2 ) 

P l(y> - P 2 W * ^1^1 L Zk)- vB n ] sin ^ 


n= 1, 3 < 


(0 < y < b) (D 106 ) 


P 5 (x) = P^(x) = Q©A 5 E\ 2 Y. fS G '( m ) sin C— 2 ) (0 < x < a) 


m=l, 3* • 


(D107) 
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N ** a© Eh 
x 


I I G ® sin ^ sin (S f J 8 n«. 


m=l, 3 * • • n=l, 3 ... 


(0 < x < a) 
(0 < y < h) 
(D108) 


(Vy=o - CVy=b = ^ \ Z T? sin (S f > <° < * < a) 


m=l, 3 • 


(D 109 ) 


N = oQEh 

y 


n 

V Y C sin (H 5 S) sin (SH) 6 nQ 

£_ Z-. mn v a b nQ 


m=l, 3 - • • n=l, 3 ■ 


(0 < x < a) 
(0 < y < b) 


(DUO) 


N c (l) 

(N ) = (N ) = -a© Eh V — /y\ sin (~^)5 

y x=0 y'x=a 4Q 2g2 $ (1) b nQ 

n=l, 3- • • n 

(0 < y < b) 


(Dill) 


N xy = “ Q0Eh Z Z J mn cos ^ C0S } 6 r 


nt=l, 3* • • n= 1, 3 • • 


(0 < x < a) 
(0 < y < b) 
(D112 ) 


With the known, equations (D13) to (Dl6) yield the following 

tensions at the ends of the x-wise stiffeners where they join the rigid 


vertical stiffeners: 


1 *B (1) 


Pj(p) = r 3 (») - P 4 (0) = P 4 (a) - <*»A ? EV, y -jS— & nQ 

n-1,3... 

(D113) 


And from the first two of equations (4), one obtains the following running 
tensions between the rigid stiffeners and the plate edges at x=0 and x=a: 


I 


(o < y < b) 


N 

N x (0,y) = N x (a,y) = 0 ©Eh £ B^ l} sin (^) b^ 

n=l, 3. • • 


(Dll4) 


Numerical results for condition (2) . - The numerical result for 
— *"0 and for square plate (B = l) with all stiffeners identical, 

subjected to a "pillow- shaped" temperature (P = Q = l) is shown in 
dimensionless form in figure 13 for Poisson's ratio v equals to 0.3* 


Illustrative Prescribed-Force Problem 
As another illustration of the application of the general theory 
of this appendix, the case will now be considered in which the stiffeners 
at x = 0 and x = a are rigid and straight and the plate is stretched 
by means of forces applied perpendicular to these stiffeners. The loading 
in this case is that of figure 4b with T^ = T^ = T ^ 0 and all other 
loads vanishing. As in the previous example, the plate is assumed to be 
isotropic and to have the same Young’s modulus as the stiffeners, while 
the stiffeners are assumed to be symmetric about the plate centerlines 
(that is A^ = Ag, A^ = A^). The temperature is assumed to be uniform 
and will therefore produce no stresses. Later on, for the sake of 
numerical calculations, the problem will be further specialized to the 
case of a square plate (b=a) with all stiffeners identical (A^ = A^ - 



Reduction of general equations to special case . - Due to the absence 
of prescribed normal forces, shearing flows, and prescribed moments, the 
following quantities are all zero: 
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(fig. 4b) 


P L (0), P L (b), P 2 (0), P 2 (b) 

M x , M 2 (fig. 4a) 

B-', B-’ (see eqs. (4)) 

<£, V' <4" (see eqs. (6)) 

(it should be noted that P^(0), P^(a), P^(0), P^(a), B^ and B^ do not 
necessarily vanish. ) And due to absence of thermal loading, the 
following quantities are also all zero: 


n' 


mit 

V 


m 


rpii ir 

m 


T 

mn 


V’, 

n' 


V” 


n 


(see eqs.(l4)) 
(see eq. (l6)) 
(see eqs. ( 23)) 


Because in this example the structure and loading are symmetrical 


about both centerlines (x = a/2 and y = b/2) it follows that c^ = c^, 
g j = g^, and B^ = B^, and the simplified system of equations, namely 
equations (D53), (D56), and (D57)> may be used for the determination of 

the b a- s a> “ d 

From the given loading conditions the quantities 5^^, defined by 


equation (D7), defined by equation (Dl8), and S' >,; , defined by 


(7) 


n 


m 


equation (D20), can be reduced to the following 

5^ = 0 
n 


= 0 
n 


s (7) = gT 
m a 


(M = 1, 3, M) 


(D7" ) 
(di8") 

(D20”) 
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Substituting from equation (D7"), (Dl 8 "), (C22'), (C23’) and (D7l) into 
equation (D 56 b), one obtains 


L = 0 


(D56"'b) 


Equations (D57)> which define the quantities g^, can be reduced to the 
following form if use is made of equations (C24'), (D7")> (D7l), (D 20 ' t ), 
(D56’a), (D56’b), (D56'c), (D57'a), and (D57'"b): 


N j N 

f y ^ y 

-rj ran 2 Z_ 


n 


n= 1, 3 • • • G ^ n 


n= 1, 3 • • • 


{ 2 . 2 2 v2 

(m + n B ) 


] 


N [■ 


■^b’ y 


2 '"'2 2 2 " “r"4 ]Cm2 + ( 2 +v )n 2 B 2 ] 


/ ^ , 2 2 >2 2 
(m + n B ) n B 


2 • - (m 2 * a 2 ) 2 y Ci " : ( 2 + ;>ff ] 

M... + 


•} G(m) 


M 


N 


2_ 

ran 


ie 0 


■ y y 

C0 

p=l, 3 * • • n=l, 3 " • • n 


mn 


V 3 


C -g-g - g 2 - ~ T ]U 2 +(2 + v)n 2 B 2 ] 

4 (m 2 +n 2 B 2 ) 2 n 2 B 4 

M 


(M 2 + n £ B 2 ) 2 y 

, Ci +n B ] 

1“ -Lj v* ■ • 


(1 ‘ 6 mP ) P 

-G^P) (m = 1,3, (D57('l)) 


m 


where 


G(p) 


g^ 

T h 


3 


With g^ known from the solution of the above equations, the c^ can be 
determined from equation (D 56 ). If use is made of equations (D56'a), 
(D56 'c), (D56'"b), and (b 66 ), then equation (D 56 ) can be put in the 
following form: 
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C(n) = 


n BC ^ 


M 

I 


m | G(m) 
ran ' ' 


n m=l, 3 . . . 


(n = 1, 3, .... M) 

(D56(l)) 


where 


C(n) = 


c 'A n 
n 1 

Th 


With the and g^ known, equations (D53), (D77)> (D7 f ) ^ (Cl 2'), (C13')> 
and (B 66 ) lead to the following expression for the determination of 


B' : 
n 


B(n) = 


M M 

f G(m) c(n) £ *4 

Cm 2 + n 2 B 2 ] 2 _ [in' 


m=l,3- . 


m=l, 3- 


2 B 2 - vm 2 ] 
2 + n 2 B 2 ] 2 


M M 

_ V [ni + ( 2 +v)n j V 

r 2 ^ 2 P n 2 r 2 . 2 TB 2 1 2 

[m+nB] , [m+nB] 

m= 1, 3 • • • m= 1, 3 • • • 


[ra 2 + (2+v) n 2 B 2 ] 


(n = 1, 3, . .., N) 


(D53(1)) 


where 


B(n) = 


B'A, 
n 1 

Th 


With the known, from equations (D13) and (Dl4) one obtains the follow- 
ing tensions at the ends of x-wise stiffeners where they join the rigid 
vertical stiffeners: 

N K 

P 5 (0) = P 3 (a) = P 4 (0) = P 4 (a) = T[i - £ ^ B(n)] (D115) 

n=l, 3* • • 

And from the first two of equations (4), one obtains the following running 
tensions between the rigid stiffeners and the plate edges at x = 0 and 
x = a: 
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N 

N x (0,y) = N x (a,y) ^ B(n) sin (^£) (0 < y < b) 

n=1,5 ’“ (Dll6) 

With S', g^, as known, equations (B6o) will give the values of 
the odd subscripted s^, s^, s^", s'”', and equations (B6l), (B3*0> and 
(B35) values of j , g^, C mn m an< ^ n odd. The expressions 

obtained are: 


s 


r 

n 


s 

j 


m 


mn 


g mn 


c 

mn 


s^ = T[C(n) - vB(n)] 

(D117) 

s"" = TG(m) 
m 

(Dll8) 

IT J(m,n) 

(D119) 

‘*1 

r 1 G(m,n) 

(D120) 

1 

p C(m,n) 

(D121) 


1 


where 


J(m,n) = 


4nB 


r 2 , 2-2 n 2 
itLm + n B J 


[n 2 B 2 B(n) - m 2 C(n) 


mn B G(m)/C] 


(D119a) 


G(m,n) = 


n[m 2 +n 2 B 2 ] 2 


[mn 2 B 2 C(n) + n^B^G(m)/C + m(m 2 + 2n 2 B 2 ) B(n)] 

(D120a) 


C ( m n ) _ CsJ + 4m 2 n B 2 G(m) 

v\,m,n; - 0 0 0 . o o o. 


r 2 2 2 -i2 r2 2—2 -i2 ~ 

rt[m+nBj nLm + nBj C 


♦ 4m B(n) r a (m 2 ^nV) . x] (D12Ja) 
n B it (m + n B ) 
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Equations (Bl6) and (B19) through (B25) then give the following 
stiffener tensions and plate stresses: 


P 1 (y) = P 2 (y) = T [C(n) - vB(n)] sin C 2 ^) (0 < y < h) 

n=1,5 **' (D122) 

M 

P 5 (x) = P 4 (x) = T £ G(m) sin (5~) (0 < x < a) (D123) 

m=l,3- • • 


N x “ r Z Z G(m ' n) Sln O Sln (S ? ) (0 < y < b) 

m=l, 3* • • n=l,3--- 

(DI24) 


n 

<Vv=0 * <Vy=b * P Z sin (=f ) (0 < x < a) 


m=l, 3* • • 


(D125 ) 


» y = p z Z c(Di ' n) sin (E ?> sin (S f > (0 < x < a) 


m=l, 3- • • n=l, 3* • • 


<Vx=0 ’ < B y>x=a * r 


^ C(n) sij 


n=l, 3* • • 


(0 < y < b) 
(D126) 


(0 < y < b) 

(DI27) 


N xy - ' P Z Z J(m ’ n) cos O 008 C T> (0 < y < b) 


m=l, 3- • • n=l, 3. . 


(D128) 


Numerical results for and \ 0 . - The equations described 

above were applied to the special case of a square plate (B = l), with 
all stiffener areas equal (A^ = A 2 = = A), and Poisson's ratio 
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v equals to 0.3* The assumption that B = 1 and that all stiffener areas 
are equal implies that X^ = X^ (see eqs . (C8l)), and the common symbol 
X will therefore be used for both X^ and X^. The Gauss-Seidel iteration 
procedure was used in solving equations (D57(l)). 

The results obtained from the stiffener tensions and plate stresses 
are shown in dimensionless form in figure l4 for X = 1.0. The values 
of M and N employed in the calculations are indicated on the figure. In 
general, stresses were computed at x/a and y/a interval of 0.1. 


Limiting case of large stiffener areas . - For the case in which all 
the stiffener cross-sectional areas are very large compared to the plate 
cross-sectional area, equations (D53).> (D7l); and (D 72 ) may be employed 
as approximations which become more and more accurate as the ratio of 
stiffener to plate cross-sectional areas approach infinity. 

Equation (D72) can be reduced to the following equation by substituting 
equations (D7"), (D20"), (D77), and (C2V) into equation (D72) with Z" 
defined by equation (D 58 b): 


G(m) 


2 _ 

nut 


where 



(D72') 


Substituting from equations (3); (CI2 1 ), (D7"), (Dl8"), (D77), and (D64'b) 


into equation (D7l)> one obtains 

M 

Z m[m 2 + (2+v) n 2 B 2 ] - 

r> r' \ r\ ^ 


C(n) - 


nC c£> 


( 1 ) 


. [ 2 , 2 2 -i2 

m=l, 3 - • • Lm + n B J 


G(m) 


M 


(D71’) 


n 


I 

m= 1, 3 • • • 


[m 2 + (2+v) n 2 B 2 ] 

r 2 , 2 2 ,2 

Lm + n B J 
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where 


C(n) = 


c'A n 
n 1 

Th 


and the G(m) are known from equations (D72 ' ) . 

With the and g^ known, equation (D53) can be expressed in 
the following form if use is made of equations (D77)* (D7' ), (C12').> 
(C13'), and (B66 ) : 


M 


B(n) = 


m[m 2 +(2+v)n 2 B 2 ] . 

r 2 2 2 -,2 ^ 

m=l, 3. . ■ [m + n B ] 


M 


I 


C(n) y 

m=l,3- 


r 2_2 2 -i 

[n B - vm J 

r 2 2.^2 -i2 

Lm + n B J 


M 

nC f. 

m= 1, 3 • • • 


[m 2 + (2+v)n 2 B 2 ] 

r 2 . 2 2 -.2 

Lm + n B J 


M 


I 

m=l,3- • 


[m 2 + (2+v) n 2 B 2 ] 
[m 2 + n 2 B 2 ] 2 


(r»53(2)) 


where 


B(n) = 


B'A, 
n 1 

Th 


With the B^ thus known, from equations >(D13) and (Dl4) one obtains the 
following tensions at the ends of x-wise stiffeners where they join the 
rigid vertical stiffeners: 


N 

P 3 (0) = P 3 (a) = P 4 (0) = P 4 (a) - T[| - (*\/4) £ B(n)/n] 

n==1,5 ‘"“ (D129) 


And from the first two of equations (4), one obtains the following running 
tensions between the rigid stiffeners and the plate edges at x = 0 and 


x = a: 


N x (0,y) = N x (a,y) 


N 


Th 

A, 


y B(n) sin (nny/b ) 


n= 1, 3 • • 


(0 < y < b) 

(D130) 
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With c^, g^, known, equations (B6o) will give the odd- subscripted 
s n' S n' S m"> S m" , and - e< l u ation (B6l), (Bjk), and (B35) the odd- subscripted 

W c mn> aS foilows: 


s^ = s" = TtC(n) - vB(n)] 


(D131) 


s’" = s"" = TG(m) 
m m K ' 


°mn 

(V°) 

(n/0) 


T 1 J(m,n) 
A 1 


(D132) 

(D133) 


Th nt \ 
8 = a(m ’ n) 


Th - 


'inn A, 


C (m, n ) 


(D13M 

(D135) 


where 


J(m,n) = q ' ^Hq o 0 [n^B^B(n) - m 2 C(m) - ran B 2 G(m)/C] (D133a) 


r 2 2 2 -,2 

jrLm +n B J 


G(m,n) - 


r 2 2^2 -|2 
jrLm +n B J 


[mn 2 B 2 C(n) + n^B G(m)/C + m(m 2 + 2n 2 B 2 )B(n)J 

(D13ta) 


F(m,n) . + B g S(°-) 


r 2 , 2 2 -i2 r 2 2 2 

jtLm+nBj rtLm+nBjC 


+ 


2 2 

4roB(n) r tn (m 

2 2 L f 2 
n B k (m 


+ 2n 2 B 2 ) 


+ n 2 B 2 ) 2 


1] 


{D135a) 


Finally, equations (Bl6), and (B19) through (B25) will give the following 
stiffener tensions and plate stresses: 
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N 


P i ( y ) = P 2 (y) = T £ [C(n) - vB(n)] sin (mty/b) (0 < y < b) 


ri— ' 1 j 3 • • • 


(D 156 ) 


M 


P^(x) = P^(x) = T ^ G(m) sin (nutx/a) 


M 


N 


Th 

N = 
x A l 


Z Z G(m, n) sin (mitx/a) sin (nxy/b) 


'm=l, 3- • • n=l, 3- • 


(0 < x < a) 

(D137) 

(0 < x < a) 
(0 < y < b) 

(D138) 


M 


(■^ x )y S _0 = (N x } y=b = f I [G(m)/C] sin (mnx/a) (0 < x < a) 


m=l, 3- • 


M 


N 


N 


= I C(m,n) sin (nutx/a) sin (nny/b) 


m=l, 3. • • n=l, 3 . . 


N 


O* y ) x ,0 = <V~a " £ Z C(n) Sin ^ /h) 


M 


n= 1,3. • . 

N 


(D139) 

(0 < x < a) 
(0 < y < b) 

(Dl40) 

(0 < y < b) 

(DlUl) 


n 


xy 


Th 

~ " A /_ 


_ (0 < x < a) 

2 ^ J(m,n) cos (nutx/a) (mry/b) (0 < y < b) 


m=l,3. • • n=l, 3* • • 


(Dl42) 


Numerical results for the case of large stiffener areas . - The numerical 
results for \ - ^ — »- 0 , obtained from the above equations, are shown in 

figure 15 for square plate (B = l) with all stiffeners identical and 
Poisson's ratio v equals to 0.3- 
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APPENDIX E 


ANALYSIS FOR ALL FOUR STIFFENERS 
WITH PRESCRIBED DISPLACEMENT CONDITIONS 


In the present appendix the case of figure 4c is considered. In 
this case all the edges of the plate are assumed to be forced into the 
prescribed shapes by means of attached rigid stiffeners (shown shaded 
in fig. 4c), which shapes are defined by known values of K^, K^, K^", 
and K^" in equations (8) through (ll). Correspondingly, the Fourier 
coefficients B*. B". B'" and B"", which describe the running tensions 
between the stiffeners and the plate edges, are now unknowns. In 


addition the loading resultants T^, M^, T^, M^, T y My T^ and 


constitute eight new knowns, supplanting P (0), P (b), P (0), ,P (0), 

1 ’ 1 <=• /< j 


P-^(a), P^(0), P^(a), which are now unknowns . 




Formulation of Boundary Condition of Prescribed Curvature 
Equations (d4) apply also in the present case. They express the 
requirement that at x=0 and x-a the curvatures of the stiffeners and 
the curvatures of the adjacent plate edges to which they are joined 
must be equal. These equations are rewritten in slightly rearranged 
form as follows: 


160 
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2 2-6 


K' = 
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C 1 I * C ^~J 

„ mn 


m=0 
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2 2-6 
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4 — • mn 
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I t - 1 )” 1 rz ‘ C 2 ( S S ' 2 
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mn 
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n „ n / mo , 


(El) 


2 rnr 
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-C (=) ](Sf) [cA-(-D‘V)(~) 


+ |(K)(M)c 1 [(c j -c 4 )(^) -c a (=S) ][g;-(-i) n <]) 

I (- i> m f - “ w - <- i > X " ][0 1 (B § )2 - c 3 ( S f )2] 

4 — 1 mn 
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Two more equations, analogous to these, will now he obtained from 
the conditions of prescribed curvature along the edges y = 0 and y = b. 
Differentiating the last of strain-displacement relations, equations (Bl), 
with respect to x, one obtains 



<$7 be 

*y _ x 

bx by 


Eliminating the strains in terms of the stresses by means of equations 
(2), and then the stresses in terms of the stress function through 
equations (b 4), this becomes 




b^by 


ci 5 F 


1 V 


de 

> 

dy 


(E2) 


2 2 

The curvatures of 5 v/bx of the edges y = 0 and y = b are therefore 


^2 

rO V - 


CI^F 




de 


(^) = (c 3 -c 4 )(~V-) - c (^) - (^) 

cbc y=0 ^ Bx c>y y*=0 By y=0 ^ y^=0 


rb 2 V 


B^F 


B 5 F. 


- (c -c 4 )(-Y~) - C (^) - 

dx y=b Bx By y-h by y-h 


be 

(sr) 


(23) 


y^b 


The terms on the right-hand side of the first of equations (E3) can be 
expressed in series form with the aid of equations (21), (b 42), and 
(BU 3 ). The result is 

M N N 


( SU = £ t(c *-v ^ 1 - c i I 


mirx 


mn 


V"'} sin 
m a 


m=l 


n=0 


n=0 


Comparing this equation with equation (10), one obtains the following 
M equations representing the boundary condition of prescribed curvature 
along the edge y = 0: 
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N 


N 


K' " 
m 


= (c.-c,) —Vo - c n y h 

'4 y a °mn 1 £_ i 


_ V'" 
mn m 


n=0 


n=0 


(m = 1, 2, ..., M) 
(E4) 


Similarly, the condition of prescribed curvature along the edge y=b is 


N 


N 


K"" = (C.-Cj — y (-l) n d - C. y (-l) n h - V"" 
m 4 y a v ' °mn 1 mn m 

n=0 n=0 

(m = 1, 2, ..., M) 


(E5) 


The unknowns and h in these equations can be expressed in 
terms of basic unknowns c^, c^, e ^ c * To accomplish this, we 
first observe that equations (B 58 ) and (B6l) can both be represented 
by the following single equation, in which any undefined quantities are 
to be considered zero: 

2 


mn 


mnn 

abE 


(T + £ “Hi [ c (-1) c" ]C + £ ^ [g’-(-irg"]C 1 } 
L mn aa n v 'n 2 bb B ra ' 'in 1 


(m^O) 


mn 


2 5 3 3 

+ IT t-TT* <“> ♦ f(2f) 


mn 


(e6) 


Furthermore, from equation (B47), (Bl4), (B50), and (B51), one obtains 

- r ~" - 1 (E7) 


*v> * b [ <; - & 


Eh = (22) T * - — [o'-(-l)“ C"]C 

(m f 0) 

(n f 0) 


2 mn ,nn . 


mn a a ' b 


.m 


n 


'n 2 


- fC^,-C-3-) n CKC 2 {2f) + «y2C 3 )(2J) (22) ] 


+ £(H2)(^)t B ..(.ir B ;][(=) c 2 + (if) (c 4 -ac 3 )] 


a v a' v b n 


n 


- |(=) 2 (1!2) 2 [b;"-(-i)V"]c, 


m 


m J 2 


(E8) 
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and these two equations can he represented by the following single 
equation in which, once again, undefined quantities are to be regarded 
as zero: 


E h 
mn mn 

(m / 0) 


(2£) 3 t + 2 2* (S2) J [ C . . (- 1 )" c"]c„ 

v b' mn a a v b' n n 2 


2-6 


no 


2 2 
mn n / nn , 




2 2 
r nn- 


+ c 2 * (Sf) (Cj, - 2 Cj)] 

(E9) 


- 1 (“M”) 2 [b;" - (-i) n b;"]o : 


m J 2 


Substituting from equations (e 6) and (E9) into equations (E4) and (E5), 
one obtains the following equations, which are analogous of equations <El) 

K i" - - V m’ + I IT T mn T [ < C 3 ' ' C lO ] 


n=0 


mn 


N 


n-0 


mn 


I FI (|(=)(f)[(c 3 -o 4 )(!H) - Cl (Sf) ][=.-C-i)V]c 2 

^V 2 (=i) 2 [ s;-(-i)\]fc 2 (”) 2 - c 3 (Sf) 2 ]c 1 3 


N 


!= I r-(“>tBA-(-i)\ ][c 2 0 - c 3 (2 e> ]0 i 

, mn 


n=l 


N 


) y 

m a Ls 


n=0 

N 


^ ^ Uc 5 - c 4 )< E ) 2 . Cl < s §) 2 ] c 2 

mn 


2-5 


I <-«" F~ -E^UC,- 0 ,,)^ - ]c 2 

~ mn 


n=0 


(E10) 
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N 


K”" = 
m 


- C l (2 5> 1 

„ „ mn 


n=0 


N 


- 0lt )(SJ) -C l( H) ][c;-(-l) m c"]C 

_ ^ mn 


n=0 


2-6 


- “ E 32 <“> -° 5 <“ ) V 

- ! = I (- 1 )" t^)[B A -(-l)V][C 2 (2|) 2 . 0 p&)\ 

* 1 TWV1 


N 


n=l 


mn 


N 


- B ' "(— ) y (- i) n 

m v a' Z— 


E - ^TT 2 [(c 5 -c,)(M) 2 - Cl (25) 2 ]c 2 


n=0 

N 


mn 


+B 7^> Y. I - 

~ mn 


2-6 


no 


tcy VO - c i't) * a 


nrt, 


n=0 


(E10) 


Formulation of Boundary Conditions of Equilibrium 
Equilibrium considerations for the stiffeners at x=0 and x=a lead 
to the following relations (see appendix D): 


N 


,(°) - - £ 


. T M 

— B ' + — - — 
nn n 2 b 


n= 1 
N 


.(a) - - £ 


S_ B" + — - — 

nit n 2 b 


n=l 


N 


T, M, 


P, (0) = y (-l) n — B' + + -r 

4 V L nir n 2 b 


n=l 

N 


T„ M 


p, (a) = y (-l) n — B” + -§ + 
4 Z_ N ' nit ti 2 


n=l 


(Ell) 
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Additional equations of the same type result from considering the 
equilibrium of the other two stiffeners. The equilibrium equations for 
the stiffener at y = 0 sure 


a • 


a 

■P-jjCO) + P 2 (0) + J N^(x,0) dx = T, 
0 

a 

/ T a 

x N^(x,0) dx = + M, 


M 


By substituting 

V X ’ 0) ’ E B m" sin (“) 

m=l 

and solving these equations, one obtains 

M 

r- a T, 

p (0) = - ) B"' + - -2 

1 4_. mn m 2 a 

m=l 


M 


T„ M, 


Z -L -j L'l- 

(-l) m — B ' " + — ^ 

mn m 2 a 


(E12) 


m=l 


Thus, in effect, the unknown P^Co) and P^O) have been expressed in ter m s 
of the known T^ and My Similarly, from the equilibrium equations for 
the stiffener at y=b, one obtains 

Ti M, 

p (b) = - y b" , f + ~4 - — 

1 L nut m 2 a 


m=l 


M 

r— , T, M, 

P '(b) = y (-l) m — B'”' + -A + Ji 
2 v ' mn: m 2 a 


(E13) 


m=l 
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Separating B', B", B’", and B"" Terms in R', H", R"\ R ,,M 

— — rr — n nr^ m — 4 n ; 

Substituting from equations (Ell) to (E13) and (B 69 ) into equation 

(B 68 ), and isolating the B', B", B"' and B"" terms, one obtains 

n n m m 


M 


R' = S' + > H [B ' " - (-l) 1 ^""] + 7’B ’ - y"B" 
n n mn m v ' m 'n n n n 


m=l 

M 


K = - V (“l) m n [B"’ - (-l)V"] - 7"B 1 + 7' "B" 

n n mnm ' m nn nn 


m=l 


N 


(El4) 


R’" = S’” + ) H [B ’ 1) 1 ^" ] + r-B’" - r"B'"' 

m m mn n N ' n mm mm 


n=l 

N 


R" " = S"" -X' (-l) n H [B' - (-l) m B"] - r"B’" + r”’B"" 
m m mn n n mm mm 


n=l 


where S^", S^" are known loading terms defined as follows: 

M 


S A - % + A ’ E 


11 


m t . _ y jl 

n n / . E 

1 mn 


„ 2 T 3 2M 3 T, 2M, 

— T + -2. - — 2 - (- l) n [— -] 

ab mn b ab ' ' L b ab J 


m=l 

M 


(E15) 


S n ‘ -K + A 2 E 2 T T n + Z ( - X) 


T, 2M, 


m JL_ 5 mn_ T + ^ + ^ , } n . 

E ab mn b ab ' ' 


m=l 


mn 


T. 2M, 

[— + — ] 
L b ab J 


(E16) 


M . 2 

g 1 tt _ A E rp t u _ 1. ^ -L mnjr 


m 


1 


3 3am ) E ab mn a ab 


mn 


T t 2M_ 2M„ 

T + ^ -(-i) m [_2 2 ] 

' a ab 


m=l 


(E17) 
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N 


T n 2M_ 


n=l 


mn 


S"" = -Q"" + A, E, — T »" + y (_p) n .^L 225- T + _J: + — 1 -(_]_) 
m m 4 4am ' ' E ab mn a ab ' ' 

(E18) 


T ? 2M 

[— + — ] 
L a ab J 


7 ^ and 7 ^ have been defined in equations (CIO) and (d 6 ) . 

r m throu S h and ^ are also known quantities and are defined by 
the following equations: 

w 

3 


r' = 


m bnur 3 3 


N 


+ A 3 E 3 ” C 3 + I iH 

, mn 


n=l 


r rr , V ( \ n 1 2 r mn \ n 

r m - bm * + L (_1) E b C ; 

„ , mn 


n=l 


N 


Pi " — — + A F p + \ ^ / HUT \ _ 

m " brrut 4 E 4 a °3 L E b C a } °Z 

mn 


and 


n=l 


n_ - I =? t=i- {=?) c. 


mn b a E A a / 2 2 2' 

mn m -n 


•<E19) 


(E20) 


(E21) 


(E22) 


System of Simultaneous Equations 

Using 7 ^ and as defined in equations (CIO) and (d 6 ), equations 

(El) can be written in the following form: 

M M 

B , 7 ( 1 ) _ B "y( 2 ) = 8 ' - y A F [B ,M -(-ljV’”] + V n' [g'-(-l) n ?"] 
n n n n n £_ mn m ^ 1 m Z_ mn Lfe m v ' b m 


m=l 


m=l 


M 

I 

m=0 


v’ [o’ 
mn n 


(-1)" c"J 

n 


(n = 1 , 2, 


,, N) 

(E23) 
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M 


M 


B' 

n n n n 




m=l 

M 


m=l 


■z<-« 


m 


v ' [c'-(-l)“ c"] (n 
mn n v 7 n v 


= 1, 2, 


m=0 


■■, N) 

(E24) 


where 6^, A^, and 5^ are known quantities, given as follows: 

M 2 2 

S' = -K’ - V’ + / — - T (— )[(C -C, )(— ) - C (— ) ] 
n n n t — • E mn' a 7 ' 3 4 7V b 7 2' a 7 


m=l 


mn 


A 


2 2 

, _ 2 mr _JL_ mn f ,mt. /®2\ 

'mn b b E a 1 b ; ~ °3^ a' J 2 
mn 

M 


(E23a) 

(E23b) 


5" = -K" - V" + Y (-l) m T (— )[(C -G,)(SJ) - C 0 (5S) ] 

n n n L-, E mn ' a 3 4 7 ' b 7 2 ' a 

, mn 

111-1 (E24a) 

Similarly, equations (E10) can be written in the following form: 


B i "r^ ^ = 6"' -V A" [B* - (-i) m B"] 

mm mm m mn n n 


N 


n=l 


N 


N 


y q n [c'-(-i) m c"j + y V-- [i’-(-DV] 

A_. mn n ' 7 n Z__ mn & m 7 to n 


n=l 


n=0 


B ' " p(2 ) 

m m 




(m = 1, 

2, . 

M) 


N 




5""- 

m 

I> 

l) n A" [B 1 
mn n 

£ 

i 

i 

ni 


n=l 




N 



N 


d- 

-l)V [c’-(-l) m c 

mn n ' 

+ 
i — i 

= c 

(-D 

n=l 



n=0 




*\ 
i — 1 

ii 

2, 

M) 
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(E25) 


(E26) 


II 



where 


2 2-6 


= (^) T - c-(^) K-^Jc, (ES5a) 

T MM ' 


r m 2) * (S f )£ I! ( ' 1)n j^V'V^ ' VT^K-E 5 ^ 0 ; 


2 2 

6’" = -K"' - V’" + Y ~~ T (~c) [(C,-Cj ) (— ) - Ch(^) ] 

m m m E mn' b' ' 3 4 /N a 1' b' 


„ _ 2 mn , 1 /Hu c (—) 2 - c (— ^lc 
A mn " a a E ^ b H V J V b' JC 1 


— re f— ' 

E lC 2 ( ' a' 
mn 


- 

1 mn 
J a 

■ n «v 2 -. mn. 

' b' J a 

2-8 

no 

b 

D” / T 

mn 

(— ) 
mn' b ' 


(E25b ) 

(E25c) 

(E25d) 

(E25e) 

(E25f) 


n_i (E26a) 

Substituting from equations (El4) into equations (b 62) to (B 65 ) 
yields the following equations: 

M 

c , df.(n)-c ,, P 1 (n) = S’ + V n CB’" - (-l) 1 ^""] + 7 ’B’ - 7 "B" 
n 1 ' n 1 ' n/^mnm ' m nn nn 


I 


2 A i; - (-D n i: 


(n = 1, 2, .... N) 
(E27) 



M 


■c'e..(n) + c"a_(n) = S" - V (-l) m n [B ’ " - (-l)^""] - 7"B ' + 7"’B" 
n l x n 2 V n mnm v ' m n n nn 


m=l 


, 2 f n ,, 2 y . -Cl- 1 )" K, , . 

b ( ~> L'- 1 ' E (n = X > 2 > H > 

. mn 

Dt=1 (E28) 

N 

ila_(m)-g"3 0 (m) = S'" + V H [B’-C-DV] + r'B"’ - r"B"" 
nn 3 2 m mn n n mm mm 


n=l 


N 


c;-(-D c 


m -n 


2,mjt\ ” n v n t , wN 

- a ( “i ) L E (m “ X ’ 2 ’ •••’ M) 

, ran 


n=l 


(E29) 


N 


-g^P_(m)+g"a. (m) = S"" - V (-1)^ J - r"B ' " + r"'B"" 

V2' / & m4 v/ m /_' / mnn' ,/ n mm mm 


n=l 


N 


c’-(-l) m c" 


2 /imt / , \» w n ' ' n / w \ 

(— ) 2^ C” 1 ) E ( m = 2, . . . , M) 


n=l 


mn 


(E30) 


Equations (E23) through (E30) constitute 4 n + 4m equations that have 
to he solved, simultaneously for the 4 n + l+M unknowns c^, c||, g^, g^, B^, 


B", B"’, and B"". 
n’ m. ’ m 


Procedure for Use of Equations 

All the pertinent equations for this case have been presented above. 

The procedure for using them will now be summarized: Equations (E23) to 

(E30) are first solved simultaneously for the c^, c^j, g^, g^, B^, B^, 

B^", and B^". With these known, equations (b6o) then give the s^, s^, 

S m"' S m"’ e< 3. ua "ki° ns (B57) to (B59) and (B6l) the j . Finally, equations 

(Bl6) and (B19) to (B25) give the stiffener and plate stresses. 
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Special Case: Square Plate, Structure and Loading Symmetri cal 

About Centerlines and Diagonals 

If the plate is square (b=a) and if the structure, loading, and 
thermal strains are symmetrical about the centerlines (x= a/2, y=b/2) and 
about the diagonals, considerable simplification of the simultaneous 
equations can be effected. Symmetry about centerlines and diagonals 
implies that A ] _ = = A, P 1 (0) - P^b) = P 2 (0) = .P 2 (b) - 

P 3 (0) = Pj(a) = P 4 (0) = P^(a), q x (y) = -q^y) = q^(x) = -q^(x), e^y) = 
e 2 (y), e^(x) = e^(x), = 0, for m even or n even, T ± = = T, 

M x = M 2 = = U k = M q , 


K' = K" = 0 
n n 


V' = V" = 0 
n n 


for n even 


K'" = K"" = 0 


m 

m 

V " = 

V"" 

m 

m 

K' = 

-K" 

n 

n 

V' = 

-V" 

n 

n 

K ,M = 

-K" 

m 

m 

V " = 

-V"' 

m 

m 


for m even 


6' = -5" for n odd 

n n 


ft’" = -6"" 
m m 


for m odd 


By the property of symmetry about the centerlines, it follows that 


c ' -c" = g’=g" = B' = B" = B"' = B"" = 0 (m, n even) 

nnTnTnnnm m 


(E31) 
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B* = B" 


n n 


B • " = b"" 

m m 


c ’ 
n 


n 


_ t "“.n 

^ ~ Ki 


1 


(m,n odd) 7 


(E32) 


In order to insure that the physical symmetry about the diagonals is mani- 
fested in the mathematical solution, M and N will be restricted to be 
equal. It can then be expected that 


B ' " 
m 


= B’ 
m 




(E33) 


Consequently, the original simultaneous equations system, equations (E23) 
to (E 30 ), can be reduced to the following: 


M 


B . [( 7 ( 1 )_ 7 ( 2 )) +2 a' ] 

n w n n ' nn 


= 6 ' + 2 


E 

P=l,3- 


P' c' 
pn p 


+ 2 c 


M 

7 v 

n pn 

P=l, 3- . . 


M 

-2 ) A* B’ (1 - 6 ) (n = 1, 3, . . ., M) 

L. pn p np' v ' 

p=1 > 5 "” (E3U) 

2 M 

{[a.(n)-p_(n)] + = S' + B'(7'-7 U ) + 2 V n B' 

n l v ' l v ' b x b' E n n' n n y pn p 

P-1,3 - • • 


nn 


M 


- !<5g> I 

p-i, 3- • • 


c'(l - 5 ) 

p P n 

E 

pn 


(n = 1, 3, •••, M) 


(E35) 
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The procedure for the special case can he summarized as follows: 

Solve equations (E3M and (E35) simultaneously for the c^ and B^ 

(n = 1, 3> •"> N). (This can he done by the Gauss-Seidel iterative 

method, which will he described in detail in a later section dealing 

with an illustrative thermal- stress problem.) With these known and 

using the relations c" = c', g* = g" = c', B" = B', B'" = B"" = B' 

n n* Tn mu nr n n m m n 

(m, n odd), equations (b 6 o) will furnish the values of s^, s^, s^", s^" 
(m,n odd), and (B 6 l) the values of the 3 ^ (m, n odd). Equation (Bl 6 ) 
and (B19) through (B25), in which only the odd values of n and m are 
included, will then give the stiffener and plate stresses. 


Limiting Case of Large Stiffener Areas 

Still referring to the special case of symmetry with respect to 

plate centerlines and diagonals, the limiting condition in which the 

stiffener cross-sectional areas A are large compared with the plate 

cross-sectional area will now he given brief consideration. This 

limiting condition can he analyzed by means of a first-order perturbation 

applied to equations (E3^0 and (E35), as follows: Assume c^ = c + kc^ 1 ^, 

B^ = + substitute these expansions into equations (E3^) and 

(E35), and equate separately terms of the zeroth and first degree in 

This will lead to four systems of simultaneous equations, two for c^'* 

and B^ and two more for c^ and B^. The latter two equations will 
n n n 

involve the c^^ and Therefore the c' ,< ” > ^ and equations must 

n n n n ^ 

be solved first, after which the c^^ and equations can be solved. 

n n 

This technique will be described in more detail in connection with 
a particular application in the next section. 
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Illustrative Thermal- Stress Problem 


In order to illustrate the details involved in the application of 
the foregoing analytical results, a particular example will be considered 
which has the following characteristics: 

a) All four edges kept straight; therefore the etc. in equations 
(8) through (ll) are all zero. 

b) Plate isotropic; therefore the elastic constants are given by- 
equations (3)- 

c) Plate and stiffeners have the same Young's modulus E. 



e) Plate with the same dimension in x and y directions; therefore 
B in equation (C79) equals to 1 (a = b). 

f) No force loading. 

g) Stiffener temperature constant at the value T . 

h) Plate temperature T(x, y) symmetrical about both centerlines 
(x = a/2, y = b/2) and varying sinusoidally in accordance with the 
following equation: 


T(x,y) = T q + 0 sin (^) 



(0 < x < a) 
(0 < y < b) 


where P and Q are odd integers. 

i) P = Q 

j ) Plate and stiffeners have the same coefficient of expansion a. 
It will be recognized that in this problem, as described above, there 
exists symmetry about the plate centerlines and diagonals. When 
numerical calculations are considered, the problem will be further 
specialized to the case of a "pillow- shaped" temperature distribution, 
namely P = Q = 1. 
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Reduction of general equations to special case . - From the given 
temperature distribution, one obtains the following equations for the 
known coefficients in terms of temperature distribution and the coefficient 
of expansion a (see appendix C): 

T' = T" = T’" = T"" = 0 (E3 6) 

n n m m v 

2 

t = -6 B 5 n ae(-) [P 2 + q 2 b 2 ] (E37) 

mn mP nQ v a ' v 

V' = -V" = a© — 8 . (EJ8) 

n n a nQ 

Si mi larly, 

V”’ = - V"" = a© ^ 5 ^ (E39) 

m m a mP 

Due to the absence of prescribed forces, the following quantities are 
all zero: 

T i> M i" T 2' M 2’ T y M y T v m 4 ( fi g- 
Q n^ Sii"' Sm" ( see eqs ’ 

(it should be noted that P^(0), P^(b), P,-,(0), P,_,(b), P^(o), P^(a), 

P, (0), P, (a), B', B", B'", B"" do not necessarily vanish.) 

Because in this example the structure and loading are symmetrical 
about both centerlines and both diagonals, the simplified system of 
equations, namely equations (E3*0 and (E35), are the governing equations. 
From equations (E 36 ) and (E37 )j (C8l), and (E15): 


S' = 8 . a© Eh 
n nQ 


PQB 

(P 2 + q 2 b 2 ) 


(EI5 ’ ) 
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Substituting from equations ( 3 ) and. (C8l) into equation (E22), one 
obtains 


n 

mn 


2mB 

n 



m 





(E22 ' ) 


Equation (E23a) can be reduced, in the following expression by substituting 
from equations {3), (C 8 l), (E37) and (E 38 ): 


6 ' - ofi (Si) r C^ * (g+v) j£ l _ D 8 
a [P^ + Q ,V] ^ 


n 


(E23'a) 


Substituting from equations (C8l) and ( 3 ) into equation (E23b), one 
obtains 

2 2 2 

_ 2mn[n B - vm ] B (E23’b) 

mn . r 2 . 2-2,2 - 

b Lm + n B J Eh 


Substituting from equation (C8l) into the first and fifth of equation 
(B 67 ), one obtains 


M 


a^(n) = AE (1 + 


± y - 

2 L c 

m=l L 


m 


2 x 2 t > 2 ! 2 

m + n B J 


M 


Pl(n) = AE -i X 


/ 2 
(-n m 




r 2 . 2_2 -,2 
Lm + n f J 


Therefore, 


M 


a L (n) - P x (n) = AE [l + ^ £ 


m 


nt=l,3- • • 


r 2 . 2„2 -,2 

Lm + n B J 


} (e4o) 


Equation (E3^ ) can be reduced to the following expression by substituting 
from equations (C8l),- (C12'), (C13 1 ), (D77), (E23*a) and (E23*b): 
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M 


, 2(B -v)B .^2 V 

[i + b 2 r _ ft 

p=l,3- • 


(g*y . ) nV] , B(n) 

_ O O 0.0 •* 


r 2 . 2 2-i2 

Lp + n B ] 


_ g ?! f [P g +(g-nv)Q 2 B 2 ] _ j _ 2 y p[p 2 +(2+v)n 2 B g ] C(p) 

2 [P 2 +Q 2 B 2 ] ft tp 2 + n 2 B 2 ] 2 

p=l,3. . . 

M M 

+2n 2 B 2 C(n) Y [ | B - g ~g P g ] -2nB 2 V - 

p,XJ... tP 2 *n 2 B 2 ] 2 p.jf,... t(p 2 + n 2 B 2 ] 2 


where 


B’ 

= o©ih' 


(n =1, 3, M) 


Q0Eh 


(E34 ' ) 


Substituting from equations ( C8l ), ( C22 ' ), (C23'), (E15'), ( E22'), and 

(E40) into equation (E35), one obtains 

2 


M 


(n[X + ^ £ 


(p 2 + n 2 B 2 ) 2 


V 

] + i — ) C{n) 

n( 1 + B 2 ) 2 


P=l,3..i 

=6 no T + Cvn + I 


M 


[P +Q B ] 


p=l,3- • • 


(p 2 + n 2 B 2 ) 2 


]B(n) 


4 * \ 


M 

I [■ 


, M 
72^2 - 5 lB < p) - V Sb2 I 


P ( 1_ ^ n p ) c ( P ) 

[p 2 + n 2 B 2 ] 2 


P=i,3... ~ ' P=l,3... 

(n = 1, 3, M) (E35 ’ ) 

Procedure for numerical solution. - In the solution of this system 

of equations (E3^ ' ) and (E35 ' ) the Gauss-Seidel iterative procedure of 

reference 6 was employed. Briefly this involved an initial assumption 

that all the B(p)'sand C(p)'s in equations (E35 ' ) were equal to zero 

except C(l). This allowed an approximate value of C(l) to be obtained. 
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Substituting this approximate value of C(l) into equations ( Ejl ' ) and 
setting B(3), B(5 )j etc. and all other C(p)'s equal to zero gave an 
approximate value of B(l). Substituting these approximate values of 
B(l) and C(l) into equation (E35 ' ) and setting all other C(p)'s except 
C(3) and all other B(p)'s equal to zero gave an approximate value of 
C(3) • Substituting these approximate values of C(l), C(3)j and B(l) 
into equations (E3^') and setting all other C(p)'s and all other B(p)'s 
except B(3) equal to zero gave an approximate value of B(3)* Continuing 
in such a fashion it was possible to obtain one set of approximate 
values for the B(p)'s and C(p)'s. This set is called the first iteration 
solution to the system of equations (E3^') and (E35’)- Additional 
iterations were obtained in the same manner as the first iteration 
except that initial values of unknowns are not assumed to be zero are 
taken from the results of the previous iteration. The procedure was 
stopped when no B(p) or C(p) changed more than 0.000001 from one 
iteration to the next. Alternatively, one can reverse the order in 
which equations (E3^ ' ) and (E35'), solving for a B(p) in equations 
(E3U ' ) first, then for the corresponding C(p) in equations (E35 ' ) • The 
final results are the same as in the above procedure. 

With the B' known, and T, = 0, = 0, from the first two of 

equations (Ell), one obtains the end tensions in x-direction of stiffener 
at y-0: 

N 

P 5 (0) = P 5 (a) = - Q9AE \ fe B(n) (Ell 1 ) 

n=l, 3- ■ • 
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Due to the properties of double symmetry and symmetry about the 
diagonals, with B’ and c’ as known, then B" = B'” = B" M = B' and 

= c^. With these known, equations- (b 6 o) will furnish the 
values of s', s". s'", s"": (B 6 l) the values of j : (B34) the values 
of g^; and (B35) the values of (only the odd- subscripted quantities 
are needed). Equations (Bl 6 ) and (B19) to (B25), with the only odd values 
of n and m included, will then give the following stiffener and plate 
stresses: 


N 


Pl(y) = Pg(y) = 0A m £ [C(n)-vB(n)] sin (^) (0 < y < a) 


n=l, 3 - • • 


(E4l) 


M 


P 5 (x) = P^(x) = 9ABCC ^ [C(p)-vB(p) ] sin (— ^ -) (0 < x < a) 

P=b3--- (E42) 


N = CtQEh 
x 


M N (0 < x < a) 

Z Z G(m ’ n) Sin (S ? ) Sin (2 f > (0 < y < a) 

m=l, 3- - • n=l, 3 . • • 

(E43) 


M 


( N x)y=0 " ( N X ) y=a ’ 00,511 Z SiD O (0 < X < a) 


P-l,3- ■ ■ 


(e44) 


N = Q©Eh 

y 

m=l, 3 . . 


M N 

Z Z 


n=l,3- - • 


„ (0 < x < a) 

C(m,n) sin (-^ )sin(2S) (o < y < b) 

(E45) 


N 


(N y>x -0 = <Vx.a = ^ Z C(n) Sin (2 f > (0 < y < a) 

(e 46) 



M N 

-Q©Eh I I J(m,n) cos (^) cos (^) 

m=l,3. . . n=l,3- • • 


(0 < x < a) 
(0 < y < a) 

(E47) 


where G(m,n), C(m,n), and j(m,n) are known quantities, and given by 
the following equations; 


G(m,n) = 


^mn 

o©Eh 


>in 2 E 2 r . „2„, M . lnn[m 2 + 2n 2 B 2 ]B(n) 

T 2 “ “2,2T2 tmC(n) + "= C < m)] + “ r 2“ 2„2^ 

rtLm+nBj jrLm+nBJ 


4nB 2 B(m) , n 2 B 2 [n 2 B 2 + 2m 2 ] 
rtm 2 [m 2 +n 2 B 2 ] 2 


1) - 


2 2 
Qja 


6 „6 


(P 2 + Q 2 B 2 ) ^ nC * 

(e48) 


C(m,n) = 


mn 

Q!0Eh 


4mB(n) r m 2 (m 2 +2n 2 B 2 ) 4nB 2 B(m) (n 2 B 2 + 2m 2 ) 

22 1 , 2^ 2 q 2v 2 " J n r 2 M 2^2,2 

n B it (ra+nB) Lm+nBj 


bnf 


rt[m 2 +n 2 B 2 ] 2 


[raC(n) + nB 2 C(m)] 


8 „8 


(p 2 + q 2 b 2 ) 


(E49) 


j(m, n) = 


v mn 

CtQEh 


PQB 


-8 ^8 


4m nBC(n) 


mP n.Q r 2 . 2^2,2 r 2 . 2^2,2 
(P +Q B ) jtLm + n B J irLm + n B J 


4m^B 


r 2 , 2 2 -i2 

itLm + n B J 


B(m) + 


4n^B 5 


r 2 2 2 -i2 

nLni +n B J 


B(n) 


(E50) 
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With B' = B" = B'" = B"" (n odd) known, from the first of equations 
n n n n v ^ 

(4), one obtains the following running tension between stiffeners and 

-i . N 

plate: i 

N x (0,y) = oGEh 2 ^ B (“) s in (-^) (E5l) 

n=l,3. . . 


Numerical results for X = 1 . - The numerical procedure and equations 
described above were applied to square plate subjected to a "pillow- 
shaped" temperature distribution (P = Q = l) and having stiffener and 
plate cross sections corresponding to X = 1.0. The results obtained 
for the stiffener tensions and plate stresses are shown in dimensionless 
form in figure l6 for v = 0.3* 


Limiting case of large stiffener areas . - Still referring to the 
special case of symmetry with respect to plate centerlines and diagonals, 
the first-order perturbation solution corresponding to stiffener 
cross-sectional area large in comparison with plate cross-sectional 
area (X — >0) will now be given. This solution is obtained by assuming 

that C(n) = C^(n) + XC^(n) 

(E52) 


B(n) = B^(n) + ^^(n) 


where 


a© Eh* 


B(n) = 


B' 

n 

a©Eh 


Substituting from equations (E52) into equation (E34')> the following 
two equations are obtained from those terms which do not involve X and 
those which do involve X, respectively: 
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IF 


- 2» 2 b 2 V Cp^a+vK^] , B (o) 

[i + B d r ^ r - 2 ~^ 2 ‘ |2 


p=l,3. • • 


[V + n B ] c 


N 

= Pit r iP g +(2+v)Q 2 B 2 3 _ ■, _ 2 V" p[p 2 + (2+v)n 2 B 2 ] C (0) (p) 

nQ 2 L r2^ _2?,2 


[P +Q B ] 


P=l,3- • 


[p c + n B^r 


+2n 2 B 2 


M 


(n) £ 
P=l,3.. 


y p[n 2 g 2 -^p 2 ]B ( 0 ) (p) ^ 

r 2 4 . 2,2,2 ' U %' 


r 2^2 2-i _ 

[n 2 n i - gnB 

[p 2 *n 2 B 2 ] 2 p= ^ 3 ... [p + n B ] 


(E53) 


p p ^ P p o 

(^^5 - 2 „¥ f [ P | g ] ) B^) (n) 

[1 + B 2 ] 2 f~t [ p 2 + n 2 B 2 ] 2 


P=l,3. • • 


= -2nB 


2 £ ptp^v^V^p) t 2 n 2 B 2 c (l) (n) y fn 2 B 2 -vp^] 
p-tls... tP 2+ °B 2 ] 2 „_f, [p 2 + n 2 B 2 ] 2 


P=i>3. • 


~2nB 


M 

? z 

P=l>3» • • 


P£ngg 2 -_v p 2 ] B (l) ( nl ( } 

[p 2 + n 2 B 2 ] 2 U V 


(E54) 


Similarly, substituting from equations (E52) into equation (E35'), one 
obtains 

C^ 0) (n) = vB^(n) (n = 1, 3, •••, N) (E55) 

and 
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L 


C (l) (n) = 5 £ g F 2 - 2 ~ + B (0) (n) n 2 B 2 £ g ' \ g -g~ 

nQ 4 [P 2 +Q 2 B 2 ] _ [p 2 + nB 2 ] 2 


y [ £_ - I] b( - ( £) - nB 2 

^ (p 2 + n 2 B 2 ) 2 P n 


p=l, 3 . 


P=l* 3 - • • 


r 2 2 2 -i2 

Lp +n B J 


(n) y — P g g : 2 - + >-B (l) (n) 

f-, P + n¥f 

P=l, 3 . • • 


The C^(n) and B^(n) are needed only in the combination C^(n) - vB^(n) 
for later stress calculations; from the last equation, one obtains the 
following expression for the combination: 


ch)( n )-vB< 1 )(n) - 6 g F g g - + B«”(n) n 2 B 2 £ ■ 2 * g g g 

^ [p +Q B ] r n + n R I 


(n) n 2 B 2 V - 

-■ L [p + n B J 


p=l, 3 . . . 


■ I t 

p=l, 3 . . . 


p 3 V 0) ' 

(p 2 + n 2 B 2 ) 2 2 " 


p=l> 3 * 


pC (0) (p) 
[p 2 +n 2 B 2 ] 2 


'<-) I 


p=i, 3 . . . 


r 2 f 2 2 -|3 ("- 1, 5, M) 

[p + n B ] 

(E56) 


Equations (E 53 ) through (E56) are to be solved for B^(n), c£°^(n), 
B^O), and C^(n). Equations (E53) and (E55) must be solved simultaneously 
for the B^^(n) and C^^(n). The procedure for solving simultaneously 
for B^(n) and C^(n) is the same as described previously for the case 
of any With B^(n) and C^(n) known, equations (E56) will furnish 
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the values of (n)- vB^ l) (n) directly. By substituting from equations 

(E52) into equations (Bl6) and (B19) "to (B25).> and neglecting terms in 
\ for the plate stresses, one obtains the following equations for 
calculating stresses: 

N 

P^y) = P^(y) = eAECt*. Y [C^ 1 ^(n)-vB^ 1 ^(n) ] sin (pp) (0 < y < a) 

n=l, 5- . . ( E57 ) 


M 

P 3 (x) = P 4 (x) - 9AmX ^ tC (l) (p) “ vB^(p)] sin (^) 

p=l,3> * • 


(0 < x < a) 
(E58) 


N 

x 


M 

O0Eh Y^ 
m=l, 3- • 


N 


I 

n=l, 3 ■ 


vnn 


sxn 


(I^X) 


sxn 




(0 < x < a) 
(0 < y < b) 

(E59) 


M 


(N ) = (N ) = Q0Eh Y C (0) (p) sin (—p) (0 < x < a) 

v x y=0 v x 7 y=a Z-* a 


p=l,3- 


(e6o) 


M N 

V y -c sin (=£) sin (Sf ) 

Z— mn a a 


N y = o9Eh ^ ^ -nn 

m=l, 3- • • n=l, 3- • • 


N 


(0 < x < a) 
(0 < y < a) 

(E6l) 


(N ) = (N ) = Q9Eh V C'°'(n) sin (^^p) (0 < y < a) 

v y'x=0 x y x= a 4-. a 

n=l, 3* • • (E62) 

M N (0 < x < a) 

\y ■ ' aeEh Z Z 3 nrn 003 O COS (B ? L) (0 < y £ a) 

111=1, 3- n=l, 3 - (e65) 

where e , c , and 3' are known quantities, and given by the 
w mn' mn 

following equations: 


185 



_ 

g mn ~ oSEh 


4n 2 B 2 


4m[m 2 +2n 2 B 2 ]B^°^(n) 


r 2 2 2 2 -ta (0l W*^ (0) W]* a 2 2,2 

sLffl + n B J irLm + n B ] 


4nB 2 B^(m) f n 2 B 2 [n 2 B 2 + 2m 2 ] _ 


3tm 


r 2 . 2 2 -,2 

[m + n B J 


n q 2 b 2 _ _ 

(P 2 + q 2 b 2 ) 

(E 64) 


- _ mn 

"mn _ o©Eh 


_ 4mB^(n) r m 2 jm 2 +2n 2 B 2 ) , 4nB 2 B^°^(m) (n 2 B 2 + 2m^) 

= n 2 B 2 « (m 2 + n 2 B 2 ) 2 31 ° ° ° " 


[m 2 + n 2 B 2 ] 2 


4'm c 


n[m 2 + n 2 B 2 ] 2 


[mC^(n) + nB 2 C^(m)] - 


5 „6 


(p 2 +q 2 b 2 ) nQ 
(e65) 



£mn_ 

O0Eh 


PgB 4m 2 nBC^°^(n) 

(P 2 + Q 2 B 2 ) mPn ^ _ ,[m 2 + n 2 B 2 ] 2 


4mn 2 B^ G^°^(m) 
jt[m 2 + n 2 B 2 ] 2 


4m^B 

r 2 . 2 P -,2 

irLm + n B J 


B^(m) + 


4n 5 B 5 


i- 2 . 2 2 v 2 

jtLm + n B ) 


B^°) (n) 


(e66) 


With the B^ known, from the first of equations (4), one obtains in 
a similar manner the following equation for the running tension between 
the rigid stiffener 'and plate edge at x = 0: 
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N 

N^(0,y) = Q)9Eh ^ B^^(n) sin (mry/b) 

n=l, 3- . . 

And from first two of equations (Ell), one obtains the end tensions 
in x-wise direction of stiffener at y-0 as follows: 

N 

P 3 (0) = P 3 (a) = -O0AEX £ jjj B (0) (n) (Ell") 

n=l, 3 - • • 

Numerical results for large stiffener areas . - The result for 
A. — >-0 for a square plate (B=l) with all stiffeners identical subjected 
to a "pillow- shaped" temperature distribution is shown in dimensionless 
form in figure 17 for Poisson's ratio v equals to 0-3 • In general, 
stresses were computed at x/a and y/b interval of 0.1. 


(0 < y < b) 

(E 67 ) 
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Figure 6. Plate stresses and stiffener tensions for the case of 
one edge held straight, pillow- shaped temperature 
distribution, \ « 2.0, and v m O.J. 

(M - 50, N - 59) 
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Figure 6. (continued) 
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Plate stresses and stiffener tensions for the case 
of one edge held straight, pillow- shaped temperature 
distribution, \ = 1.0 and v = 0.3. (Dashed curves, 
from fig. 5c of ref. 1, are for the case of all edge 
stiffeners perfectly flexible.) 

(M - 30, N = 59) 
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one edge held straight, pillow- shaped temperature 
distribution, *- 0 , and v = 0. (M = 79, N = 79). 
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Figure 9. (continued) 
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Figure 10. Plate stresses and stiffener tensions of the case of two 
opposite edges held straight, pillow-shaped temperature 
distribution, \ = 2.0, and v = 0.3. n(M =59, N = 59) 
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Figure 12. Plate stresses and stiffener tensions for the case of 

two opposite edges held straight, pillow- shaped tempera- 
ture distribution, X = 1, X — *-0, and v = 0.5* 

(M = 59 , N = 59) 
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Figure 13. Plate stresses and stiffener tensions for the case of 

two opposite edges held straight, pillow- shaped tempera- 
ture distribution, — ► 0, X — *-0, and v = 0.3* (Dashed 

curves, from fig. 5a of ref.l, axe for the case of all 
edge stiffeners perfectly flexible. ) (M = 79: N = 79) 
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Figure 13* (continued) 
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Figure 14* Plate stresses and stiffener tensions for the case 

of two opposite edges held straight, stretching forces 
applied normal to the straight edges, \ = 1.0, and 

v = 0.3. (M = 59, N = 59) 







MM! 

m sniSRi 


1(91; 




iHIHISliiSlii lt 

|1M| 

r 


in 


N (0,0)A 


■ » ■ m 

! IfSi 




■iii 


■ a a 

* a I 

:(! 


in 
IIM 

£«< 

iip 

■ a i 

M I ■■ • 

i liTpi ape 
a IbbbB aha 
■ ■■■ ■■■ 

■■■■■■■■■ ■■■■■■! 
■■■■■■■■a ■■■■•■ 

■ ■■■iiibi an 


II 


III Ml 
■llllllll 


I Hi!! 


Illllllll II 


m 

iii 


ill IIIIIIHIlii 


iiiiiiiiiiiiinniii! 

■iiifiiimiiMM 

ills 

MMiiMHHBS 


Hill lllllllllllillliiillliilllllllslllllllll!ll 


Hllllllillllljlill 


i 


m 

iiiE 


i iiiiniiiifiiiiim r 


III Hi 


llil 

IIHI 


!lll 


1::::!:::::!::!:::!:!::: 
■ ■■■■ (BIBaaMlI MBBIBIIBI 
■■■■■•■■•■■■■BBBBBaaiaisi 


0.4 0.5 


Figure 15. Plate stresses and stiffener tensions for the case of two 
opposite edges held straight, stretching forced applied 
normal to the straight edges, and v = 0.5. 

(M = 59, N = 59) 
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Figure 15* 

(Continued) 
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(f) Demonstration of singular corner shear stress behavior. 

Figure 15- (continued) 
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Figure l6. Plate stresses and stiffener tensions for the case of all 
four edges held straight, pillow- shaped temperature 
distribution, A. = 1.0, and v = 0.3* (Dashed curves, from 
fig. 5c of ref. 1, are for the case of all edge stiffeners 
perfectly flexible.) (M = 59# N = 59) 
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Figure 16. (continued) 
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Figure 16. (continued) 
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Figure 19- Results from figure 11 compared with those obtained when 
overall expansion in x- direction is completely prevented 
by the external constraining forces T-^T^-O . 4439 ©CtehE. 
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